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ABSTRACT 

Foss,  John  F .  ,  Ph.  D.  ,  Purdue  University, 

January  1965.  A  Study  of  Incompressible  Bounded 
Turbulent  Jets.  Major  Professor;  Dr.  J.  B.  Jones 

The  intent  of  this  investigation  was  twofold. 

A  portion  of  the  study  was  directed  t  the  establish¬ 
ment  and  investigation  of  the  three-dimensional  effects  of  a 
single  bounded  jet  flow  between  parallel  plates.  From 
velocity  traverses  taken  between  the  bounding  plates,  the 
existence  of  an  important  three-dimensional  effect  was 
observed.  A  physical  flow  model  was  postulated,  the  vortex 
stretching  hypothesis,  to  explain  this  observation.  A  series 
of  experiments  were  conducted  to  examine  predictions  made 
on  the  basis  of  the  model;  the  agreement  between  the  observed 
and  predicted  mean  flow  characteristics  and  distributions 
provides  substantial  support  for  the  vortex  stretching  hypothe¬ 
sis.  A  capsule  description  of  the  vortex  stretching  hypothesis 
follows . 

Vortex  filaments  will  be  generated  in  the  velocity  gradient 


of  the  bell- shaped  jet  velocity  profile.  The  vortex  filaments 
which  will  end  on  the  bounding  plates,  are  always  formed 
from  the  same  fluid  particles  until  they  are  dissipated  by 
viscous  action.  Because  of  the  no- slip  condition  at  the  wall, 
the  filament  is  both  stretched  and  reoriented  as  the  central 
portion  is  swept  downstream.  This  reorientation  induces 
an  inflow  to  the  jet  near  the  midplane  between  the  bounding 
plates  and  an  outflow  near  the  bounding  plates.  The  stretch¬ 
ing  increases  the  vorticity  and  thus  intensifies  the  effect. 

The  second  portion  of  the  experimental  program  was 
an  investigation  of  the  jet  interaction  problem  which  is  the 
governing  factor  in  all  proportional  fluid  jet  amplifier  device 
and  is  an  imp>ortant  factor  in  the  operation  of  bistable  ampli¬ 
fier  elements.  Because  few  data  were  available  in  the  litera 
ture  concerning  detailed  flow  phenomena,  the  investigation 
of  the  intersecting  jet  problem  was  also  an  exploratory  one. 
Four  studies  were  conducted  with  the  intent  of  establishing 
"information  bench  marks"  from  which  further  information 
could  be  inferred.  The  tests  were; 

i)  A  jet  deflection  study.  The  measured  quantities 
were;  the  deflection  angle,  the  apparent  pivot  point  and  the 
static  pressures  assisting  in  the  deflection  for  momentum 


xiv 


flux  ratios  of  0.  05.  0.10,  0.15  and  0.  20  and  for  three  geo¬ 
metries. 

ii)  A  detailed  mean  flow  survey  in  which  the  velocity 
profiles,  flow  angles,  static  pressure  distributions,  and  mass, 
momentum,  and  energy  flux  ratios  were  determined  for  seven 
planes  between  the  bounding  plates  at  three  longitudinal  loca¬ 
tions. 

iii)  and  iv)  Two  series  of  tests  to  establish  the  effect  of 
the  nozzle  geometry  and  the  effect  of  the  momentum  flux 
ratios  were  conducted. 

From  the  data  of  study  (i)  it  was  determined  that  the 
geometry  used  caused  a  significant  increase  in  the  deflection 
angle  with  respect  to  the  momentum  flux  action  alone.  Also, 
the  analytical  model  which  was  proposed  provides  a  satisfactory 
prediction  method  for  the  geometry  with  the  highest  gain. 

From  ii),  iii)  and  iv)  it  was  determined  that  the  inter¬ 
section  of  the  jets  leads  to  an  increase  in  the  vortex,.stretching 
effects  and  that  the  nozzle  geometry  and  the  momentum  flux 
ratio  have  significant  effects  on  the  resultant  jet. 


INTRODUCTION 


Since  1960  considerable  effort  has  been  expended  on  purely 
pneumatic  elements  for  control  and  logic  systems.  The  oroportional 
fluid  jet  amplifier,  which  has  been  successfully  used  in  certain  con¬ 
trol  systems,  is  such  an  element.  Figure  1  shows  a  schematic  dia¬ 
gram  of  a  typical  proportional  amplifier,  and  by  use  of  terms  defined 
in  the  figure,  a  qualitative  description  of  the  operation  of  the  device 
can  be  given  as  follows. 

For  a  zero  difference  in  the  control  jet  flows,  the  power 
jet  flow  is  evenly  divided  by  the  splitter  and  the  flows  from  the  two 
output  ports  are  equal.  If  the  control  jet  flows  are  unequal,  the 
imbalance  will  deflect  the  power  jet  such  that  the  flows  from  the  out¬ 
put  ports  are  no  longer  equal.  The  difference  in  mass  flux,  momen¬ 
tum  flux,  energy  flux,  or  pressure  of  the  two  outlet  streams  is  related 
to  the  difference  in  the  same  quantities  at  the  control  jet  input  stations. 
From  a  comparison  of  the  flow  conditions  at  the  inlet  and  outlet,  the  ga 
of  the  amplifier  is  defined.  For  example,  the  mass  flux  gain  —  G  — 
is  commonly  defined  as 


(-1)  [  m 


G  = 


left 


m  .  .  ] 

__^i^htJ_out£i^ 


[*1 


eft 


-ml 

right'  control  jets 


where  m  is  the  mass  rate  of  flow. 

Comparing  the  purely  pneumatic  device  with  the  electronic  amplifier 
shows  that  the  power  jet  flow  is  analogous  to  the  plate  current  and 
the  control  jet  serves  the  same  function  as  the  grid  circuit.  If 
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simplifying  assumptions  are  made  for  the  form  of  the  velocity 
profiles  (both  entering  and  leaving  the  mixing  region)  and  if 
the  static  pressure  distribution  is  similarly  dealt  with  by 
assumption,  it  is  possible  to  predict  quantitatively  the  per¬ 
formance  of  the  proportional  fluid  jet  amplifier.  *  Much  of 
the  effort  to  analyze  this  particular  device  has  been  directed 
toward  this  approach. 

From  a  designer's  viewpoint,  this  method  of  pre¬ 
dicting  the  performance  would  be  greatly  improved  if  the  mean 
velocity  distribution  leaving  the  interaction  region  were  known. 
In  addition,  if  other  details  of  the  flow  leaving  the  interaction 
region  (e.g.,  static  pressures  and  turbulence  velocities) 
were  known,  this  information  could  provide  valuable  insight 
into  the  actual  interaction  phenomena.  Finally,  if  detailed 
measurements  were  taken  in  the  interaction  region  itself, 
then  the  character  of  the  momentum  interchange  and  the  role 
of  the  static  pressure  might  become  known.  This  would 
greatly  facilitate  the  design  of  such  elements. 

Figure  1  shows  one  of  the  many  possible  geometries 
which  can  be  employed  for  the  proportional  fluid  jet  amplifier; 
however,  other  configurations  may  also  be  employed  for  this 
device.  The  interaction  of  two  jets  is  the  governing  phenomenon 
which  is  common  to  all  proportional  amplifiers  and  Is  an  im¬ 
portant  factor  for  the  bi-stable  elements  (oscillators,  and 
switches);  therefore,  any  insight  into  this  interaction  phenom¬ 
enon  can  lead  to  a  better  understanding  of  the  operation  of 
all  such  devices. 


*  This  method  would  employ  the  conventional  control 
volume  approach  (see  e.g.  Reference  1,  Shames) 
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The  intent  of  the  present  work  is  to  provide  infor¬ 
mation  that  will  be  important  to  the  understanding  of  the  jet 
interaction  problem  from  both  a  fundamental  and  a  design 
viewpoint.  In  order  to  separate  the  chavacteristics  of  the 
flow  which  result  from  the  interaction  of  the  two  jets  from 
those  characteristics  induced  by  the  downstream  geometry, 
the  experimentol  configuration  was  that  of  Figure  1  with 
the  section  above  the  line  A  -  A'  completely  open.  An 
aspect  ratio  of  6  was  used. 

A  characteristic  of  all  three-dimensional  tur¬ 
bulent  shear  flows  is  that  the  system  of  equations  describing 
the  flow  is  incomplete.  That  is,  for  the  four  basic  equations 
of  continuity  and  momentum,  there  are  ten  unknowns,  viz, 
three  mean  velocity  components,  the  static  pressure,  and 
six  independent  components  of  the  Reynolds  stress  tensor, 

*  P  .  Thus  to  render  the  system  of  equations  mathe¬ 

matically  tractable,  one  must  make  certain  simplifying 
assumptions.  A  plausible  assumption  for  this  geometry 
would  be  that  the  flow  has  a  two-dimensional  character  over 
a  significant  region  between  the  bounding  plates.  If  the 
two-dimensional  assumption  is  applicable  in  this  region, 
then  the  phenomenological  theories  (e.  g. ,  Prandtl' s  mixing 
length  theory)  may  be  useful  in  gaining  an  analytical  solution 
for  the  mean  velocity.  An  early  portion  of  the  experimental 
program  was  designed  to  check  this  assumption  of  two-dimensionality 

*  Since  typical  aspect  ratios  (height  to  width  ratio  of 
the  power  jet  nozzle)  for  the  fluid  jet  amplifiers  are  from 
6  to  8,  the  flow  is  three-dimensional. 


by  means  of  velocity  traverses  between  the  plates,  viz.  u(z) 
for  a  fixed  x  value  and  various  y  values  (see  Figure  1  for 
coordinates) .  These  data  were  taken  for  the  single  bounded 
jet  case  in  the  absence  of  side  jets. 

The  remaining  experimental  investigation  was  directed 
to  the  following: 

i)  mean  velocity  and  static  pressure  profiles  of  the 
jet  which  emerge  from  the  interaction  region  (resultant  jet), 

ii)  deflection  angles  of  the  resultant  jet  as  a  function 
of  (a)  the  ratio  of  the  control  and  power  jet  momentum  fluxes 
and  (b)  the  nozzle  geometry. 

iii)  the  effects  on  the  character  of  the  resultant  jet 
which  are  caused  by  different  nozzle  geometries  and  different 
momentum  flux  ratios. 
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NOMENCLATURE 


app 

B 

b 

cj 

C_ 


E 

hs 

k 

k' 

K 

1 

m 

M 

mfr 

N„ 


nozzle  exit  width 

apparent  pivot  point  of  resultant  jet 

barometric  pressure  in  inches  of  mercury 

/•  • 

momentum  flux  thickness,  —  \  (u/u_)^  d  y 

2 

control  jet 

^  ■  ^ATM 

pressure  coefficient,  ,  (except  as  noted) 

2 

energy  flux  thickness,  I  (q/qj.)(u/uj.)  d  y 

energy  flux,  e  qc  u^ 

total  pressure  head 

static  pressure  head 

ratio  of  geometric  lengths  Y/X 

ratio  of  geometric  lengths,  Y'/X' 

proportionality  factor  for  velocity  computation,  15,  89'^T/B) 
mixing  length 

mass  flux  width  \  u/  u  dy 
mass  flux,  p  m  u^, 

2  2 

momentum  flux  ratio,  a,  ut  /  a,  u, 

*  t  in 


Reynolds  number 


static  pressure 


power  jet 

dynamic  head,  (h^  -  h^);  except  for  section  VII-A 
where  q  c  momentum  flux  per  unit  depth 
radius  of  curvature  of  jet 

dial  reading  for  the  x  -  y  -  G  traverse  device 
X  -  component  of  mean  velocity  vector 
t  component  of  mean  velocity  vector 
y  -  component  of  mean  velocity  vector 
magnitude  of  mean  velocity 
z  -  component  of  mean  velocity  vector 
longitudinal  distance  (see  figure) 
apparent  origin  of  the  jet 
length  of  standoff  wall  (see  figure) 

X  +  a/  2 

transverse  distance  (see  figure) 
length  of  setback  wall  (see  figure) 

Y  +  a/  2 

transverse  displacement  of  the  centerline  for  the 
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vertical  distance;  x  -  y  -  z  form  a  right  hand  ortho- 
ganol  co-ordinate  system  (see  figure) 

X 


standoff  wall 
setback  wa 
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T 

X 


“i - 

1 

^  ^  y  A  ^ 

z 

♦  ,  « 

1 

1 

9  9  9  9  19 

i  '  ^ 


a 

y 

6 

( 

f 

ij 


Greek  Symbols 

deflection  angle  of  the  resulUnt  jet  (see  Figure  29) 

angle  used  in  analytical  study  (see  Figure  29) 

angle  used  in  analytical  study  (see  Figure  29) 

length  ratio  x  /x' 
app 

eddy  viscosity 

longitudinal  co-ordinate  measured  along  the  axis  of 
the  resulUnt  jet  (origin  at  the  power  jet  nozzle  exit) 
non-dimensional  co-ordinate,  y/ b  or  y  /b 
«/b. 


P 


density  of  fluid  (room  air) 
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O’ 

T 

e 

Q* 

{ 


app 

c 

L 

m 

0 

R 


1 

2 


* 


Standard  deviation 
shear  stress 

flow  angle,  arc  tan  (v/ u) 

flow  angle  with  respect  to  resultant  jet,  0  -  Gc 
transverse  co-ordinate  measured  with  respect  to 
the  resulUnt  jet  ( t,  4  .  *  form  an  orthoganol  right 
hand  co-ordinate  system. ) 


Subscripts 

quantity  at  the  apparent  pivot  point 

quantity  at  the  maximum  velocity  point  of  the  profile 

left  pocket  of  the  defined  region  (for  the  present  study 
this  is  the  near  pocket) 

quantity  at  the  maximum  velocity  point  of  the  profile 
and  for  z/  a  =  0 

quantity  at  the  power  jet  nozzle  exit  (except  for  x  ) 

o 

right  pocket  of  the  defined  region  (for  the  present 
study  this  is  the  far  pocket) 
power  jet  nozzle  flow 
control  jet  nozzle  flow 


quantity  referred  to  the  t,  4.  z  co-ordinate 


system 


Superscripts 

turbulent  fluctuation  of  (  ) 
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(  ) 
r-) 


time  average  of  (  ) 

vector  quantity  of  (  ) 
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SURVEY  OF  THE  LITERATURE 

Introduction 

Although  purely  pneumatic  devices  are  a  recent  develop¬ 
ment,  a  sizeable  literature  dealing  with  them  is  already  available. 
Commercial  applications  account  for  most  of  the  interest,  and 
this  product  orientation  is  primarily  responsible  for  the  litera¬ 
ture  emphasis  on  the  gross  performance  or  system  characteristics 
of  the  devices.  References  2  and  3  are  representative  of  the 
majority  of  the  current  literature.  Although  these  papers  provide 
valuable  information  for  the  systems  designer,  they  provide  little 
insight  into  the  details  of  the  interaction  phenomena. 

Significant  simplifications  in  the  analytical  model  for  the 
intersecting  jets  problem  can  be  effected  if  the  two-dimensional 
approximation  is  valid  for  the  central  region  of  the  flow  field. 

Since  previous  investigations  (17  -  21)  of  jet  interaction  between 
parallel  plates  have  made  use  of  the  two-dimensional  jet  studies, 
a  review  of  the  two-dimensional  jet  flow  literature  is  given  below. 

The  Two-Dimensional  Jet 

In  1926,  Tollmien  (4)  applied  the  then  recently  developed 
Prandtl  mixing  length  theory  to  the  two-dimensional  jet  flow  problem. 
To  achieve  an  analytical  solution  for  the  mean  velocity  distribution 
he  introduced  the  following  assumptions: 

i)  The  mixing  length  is  constant  in  any  plane  normal  to  the 
mean  flow. 
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ii)  The  static  pressure  gradient  in  the  direction  of  flow 
is  negligible. 

iii)  Viscous  shear  is  negligible  in  comparison  with  the 
Reynolds  (or  turbulent)  stress. 

iv)  The  instantaneous  velocity  fluctuation  in  the  direction 
of  flow  is  very  small  in  comparison  with  the  mean  velocity  com¬ 
ponent. 

v)  The  flow  is  in  a  self-preserving  state. 

With  the  above  assumptions  and  the  defining  equation  for 
mixing  length,  the  system  of  equations  (continuity  and  momentum) 
is  determinate. 

Starting  with  the  continuity  and  momentum  equations  and 
using  u(x,  ij),  where  jj  is  the  similarity  variable,  y/ g(x) ,  Tollmien 
was  able  to  form  the  solution  for  the  velocity,  u(x,  r>) .  g(x)  is 
a  characteristic  piarameter  which  describes  the  jet  as  a  function 
of  X  (the  velocity  on  the  jet  axis  --  u^(x)  --  is  such  a  parameter) . 

The  solution,  which  is  not  readily  presented  analytically, 
appears  in  tabular  form  in  Reference  4, 

From  an  analysis  of  the  transverse  equation  of  motion, 
Tollmien  predicted  that  the  static  pressures  in  the  jet  would  be 
greater  than  atmospheric.  He  also  applied  the  conservation  of 
momentum  criterion  to  show  that  u  (x)  is  proportional  to  x  t 
downstream  of  the  potential  flow  core. 

Two  other  phenomenological  theories  were  applied  to  the 
two-dimensional  jet  problem.  Gortler  (5)  postulated  a  constant 
eddy  viscosity,  c  ,  across  the  mixing  layer  and  after  making 
assumptions  as  those  of  Tollmien,  ii) ,  iv) ,  and  v) .  he  obUined 
the  following  expression  for  the  mean  velocity: 

—  e  sech^  A 
u  2 
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where  ri  •  yl  x  and  A»(i/x)c  • 

Reichardt  (6)  started  with  the  fact  that  the  Gaussian  error 
curve  provides  an  accurate  representation  of  the  experimentally  de¬ 
termined  non-dimensional  velocity  profile.  The  one-dimensional 
diffusion  equation  yields  a  solution  of  this  form  and  this  led  Reichardt 
to  make  the  following  analysis.  Starting  with  the  x-direction  mo¬ 
mentum  equation  (neglecting  viscous  effects  and  the  x-direction 
static  pressure  gradient), 

~  (u^  )  +  (uv)  •  0 

ox  8  y 

he  made  the  substitution 


uv  •  -  X(x)  8(u*  )  /  8  y 

and  thus  obtained  the  one-dimensional  diffusion  equation 


8  u^  ...  8*  u 


T  =  "W 


Note  in  the  above  equations  that  u  is  the  instantaneous 
turbulent  velocity  and  that  X(x)  is  independent  of  y.  The 
assumption  implied  by  the  substitution  is  that  the  transport  of 
the  x-momentum  by  the  velocity  transverse  to  the  man  flow  is 
propoxtional  to  the  gradient  of  the  x-direction  momentum  flux 
in  the  transverse  direction.  This  assumption  serves  as  the  basis 
of  the  Reichardt  Inductive  Theory.  The  non-dimensional  velocity 
profile  for  the  Reichardt  theory  is: 


u/  Uj.  e  exp  (-ir  /  8)  Tj*  * 


*  The  constant  {-ir/  8)  was  experimentally  determined  by 
Miller  and  Comings  (9). 
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where:  tj 


y/  b(x)  and  b(x) 


(u/  u^)  *  dy. 


The  above  phenomenological  theories  have  all  been  shown 
to  be  lacking  in  complete  physical  rigor.  Hinze  (7)  provides  an 
excellent  summary  of  the  phenomenological  theories  and  their 
shortcomings  (seepages  277-293). 

In  1934,  Forthman  (8)  conducted  an  experimental 
program  motivated  by  the  analysis  of  Tollmien.  The  accuracy 
of  the  results  of  this  investigation  are  uncertain  inasmuch  as 
the  static  pressure  probe  used  by  Forthman  gave  a  reading 
which  was  greater  than  atmospheric.  It  is  interesting  to  note 
that  Forthman*  s  first  probe  (an  L-head,  closed  end  type) 
gave  negative  gage  static  pressure  readings  in  the  jet.  Be¬ 
cause  Tollmien*  s  analysis  had  incorrectly  predicted  positive 
8^8*  pressures,  a  new  probe  which  did  give  positive  press¬ 
ures,  was  built  and  used. 

Miller  and  Comings  (9)  have  experimentally  verified 
the  relation  proposed  by  Townsend  (10),  viz,,  P(y)  +  p  v*  m 
^(atmospheric)*  This  formula  is  based  on  an  order  of  mag¬ 
nitude  analysis  of  the  transverse  equation  of  motion.  The  re¬ 
lation  and  its  verification  show  that  the  static  pressures  in  the 
jet  are  less  than  atmospheric. 

The  work  of  Miller  and  Comings  also  shows  that  the  experi¬ 
mental  data  for  the  mean  velocity  profile  are  well  matched  by 
the  Gaussian  error  curve  and  that  the  prediction  of  Tollmien  for 
the  center  line  velocity  decay  is  accurate  to  within  three  percent. 
The  velocity  profile  of  Gortler  is  also  presented  by  Miller  and 
Comings  along  with  their  experimental  data.  For  the  central 
portion  of  the  jet  (u/  Uc  >  0.3)  the  agreement  is  quite  good; 
however,  for  the  outer  portion  the  percentage  error  of  the 
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"t  '>'>1 


C 


1/2 


corrected  for  u*  . 


An  analytical  contribution  of  (11)  is  that  of  calculating  the  eddy 
viscosity  as  a  function  of  y  for  an  assumed  Gaussian  distribution 
without  recourse  to  turbulent  stress  measurements. 

Albertson,  et  al.  ( 12)  have  also  investigated  the  two-dimen¬ 
sional  jet.  Their  data,  in  agreement  with  the  results  of  (9)  and 
(11),  indicate  that  the  Gaussian  error  curve  provides  a  satisfactory 
fit  to  the  experimental  mean  velocity  profile.  However,  it  should 
be  noted  that  the  profiles  of  Albertson,  et  al.  show  a  much  larger 
scatter  in  the  data  than  the  profiles  of  either  (9)  or  -(11),  The 
authors  have  shown  that  the  velocity  profiles  in  the  mixing  zone 
adjacent  to  the  potential  core  region  (i.  e.  the  zone  where  Ug  »  u^) 
may  be  well  approximated  as  a  single  curve  by  the  relation 

log,o  Z  *  +  y  ] 

10  Uq  X  ' 

for  the  majority  of  the  jet  profile  (i.  e.  u/  uq  »  0,  1).  The 

lateral  component  of  velocity  was  computed  from  the  measured 

magnitude  and  direction  of  the  mean* velocity  vector.  Neglecting 

dp 

and  all  components  of  the  Reynolds  stress  tensor  except 
-  p  W,  they  were  able  to  compute  the  t  (y)  distribution.  From 
this  and  the  known  velocity  profile  they  computed  the  mixing 
length  and  eddy  viscosity  as  functions  of  y. 

Three  important  characteristics  of  the  two-dimensional 
jet  are  presented  by  Albertson  et  al. .  Plots  of  the  following 
ratios: 
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M 

—  s  (mass  rate  of  flow  at  x  =  x)  /  (mass  rate  of  flow  at  x  =  0) , 
0 

2bq 


(momentum  flux  at  x  =  x)  /  (momentum  flux  at  x  s  0) 


—  c  (kinetic  energy  flux  at  x  =  x)  /  (kinetic  energy  flux  at  x  =  0) 

0 

show  that  downstream  from  the  potential  core,  M/  M  increases 

0 

exponentially  with  distance  (viz. ,  (M/ M^)  =(x/a)  '  ),2bq^/aqQ 
remains  constant,  and  E/  Eq  decreases  exponentially  with 
distance  (viz. ,  E/  E^  =  (x/  a)  ) .  The  second  impor¬ 

tant  characteristic  is  the  angular  deflection  of  the  flow 
from  the  axis  y  =  0.  The  velocity  vector  was  found  to  have 
a  maximum  outward  angle  of  10  degrees  (away  from  the  x- 
axis)  in  the  neighborhood  of  y/  x  =  0.  15;  angles  approaching 
90  degrees  inward  were  found  for  the  extremities  of  the  jet. 

A  third  characteristic  of  the  two-dimensional  jet  is  the  length 
of  the  potential  core,  i.e.  the  region  where  u  =  u„. 

Albertson  et  al.  indicate  that  this  potential  core  closes  at 
approximately  fi /e  nozzle  widths;  Miller  and  Comings'  data 
indicate  that  this  length  is  about  six  nozzle  widths  and  the 
data  of  V.  d.  H.  Zijnen  indicates  that  the  core  closes  at  4  to 
4.  5  nozzle  widths.  Although  these  differences  may  appear  in¬ 
significant,  the  prediction  of  the  downstream  value  of  u  is 

c 

strongly  dependent  on  the  x  location  of  the  core  termination. 

Since  the  physical  flow  property  u^  is  not  exactly  constant 
in  this  core  region,  as  indicated  by  the  data  in  (9),  (11), 
and  ( 12) ,  it  is  difficult  to  establish  an  exact  value  for  the 
length  of  the  region.  The  measurements  of  Albertson,  et  al. 
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covered  a  Reynolds  number  range  of  1,  600  to  13,  000.  Miller  and 
Comings'  data  are  for  a  Reynolds  number  of  19,  000  and  v.  d.  H. 
Zijnen*  s  data  are  for  a  Reynolds  number  of  13,  300. 

In  a  separate  study  using  the  same  flow  system  as  in  (11)  V.  d. 
H.  Zijnen  has  also  made  turbulence  measurements  in  a  plane  jet  of 
air  (13).  Measurements  of  11^,  as  a  function  of  y  and 

,  ^  ,  as  a  function  of  x  for  y  *  0  are  presented.  The  trans¬ 
verse  measurements  were  made  for  six  x/a  values  from  17.5  to 
35.  The  u^ ,  vlv  measurements  indicate  reasonably  similar  dis¬ 
tributions,  whereas  7^  shows  a  large  scatter.  This  may  be  a 
direct  result  of  the  nozzle  geometry.  In  (13),  v.  d.  H.  Zijnen 
shows  that  turbulence  measurements,  made  with  the  hot  wire 
and  with  the  diffusion  technique  are  in  agreement  when  the 
(^)‘  data  are  non-dimensionalized  using  (  u  +  ("u^  )‘/^  ]  in¬ 
stead  of  the  mean  velocity  u. 

Heskestad  (14)  has  made  mean  flow  and  turbulence  measure¬ 
ments  in  a  two-dimensional  jet.  The  geometric  configuration  of 
the  flow  system  is  different  from  that  of  any  of  the  previous 
investigators,  A  channel,  which  was  60  inches  by  1  inch  in 
cross-section  and  60  inches  long,  was  placed  between  the 
plenum  chamber  and  the  nozzle.  The  nozzle,  with  an  outlet 
cross-section  of  60  inches  by  inch,  had  straight  walls 
making  an  angle  of  45  degrees  with  the  larger  walls  of  the  approach 
channel.  No  reason  was  given  for  the  long  approach  channel  up¬ 
stream  of  the  nozzle. 

The  major  contribution  of  this  paper  is  turbulence  measure¬ 
ments  for  very  large  x/a  values.  The  rms  fluctuations,  d^,  a 
turbulence  energy  balance,  the  triple  moment,  intermittency 
and  the  value  (a[u^]‘/^  /  9  t)  ^  were  measured.  A  brief 
comparison  between  the  turbulence  data  available  from  (9), 

(13),  and  (14)  is  presented  in  the  following  table. 
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Tab)e  I 

Turbulence  data  from  the  literature 


Investigator 

u^  max/  u  * 
c 

v^  max/  u  * 
c 

uv  max/ 

Reference  (9) 

0.  068 

0.055 

0.  025 

(13) 

0,  032 

0.  023 

0.  018 

(14) 

0.070 

0.035 

0.  017 

Note:  The  values  are  all  approximate. 


A  more  complete  comparison  for  u^  is  given  in  the  section 
"Single  Bounded  Jet  Study." 

The  important  differences  among  the  data  sources  are:  (i) 

the  values  from  (9)  and  ( 13)  were  taken  at  40  and  35  nozzle  widths 

whereas  (14)  is  for  10 1  nozzle  widths,  (ii)  thev^/u*  and 
TT  c 

uv/  u^^  values  from  (9)  were  caluculated  from  the  two-dimensional 
equations  of  motion  and  the  continuity  equation  whereas  the 
data  of  ( 13)  and  ( 14)  are  taken  from  direct  measurements. 

Angular  deflections  of  the  mean  flow  and  the  mean  velocity 
profiles  as  a  function  of  x/  a  v/ere  also  investigated  by  Heskestad, 

If  the  mean  velocity  distributions  are  compared  with  those  of 
previous  investigators  the  effect  of  the  unique  upstream  geometry 
may  be  inferred.  Most  significantly,  the  mean  velocity  profiles 
do  not  indicate  a  self-preserving  form  before  x/ a  =  65,  whereas 
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Reference  9  shows  that  self-preservation  is  obtained  by  x/  a  =  10. 
Heskestad  compares  his  self-preserving  mean  velocity  profile 
with  only  that  of  Gbrtler;  the  data  from  (14)  are  in  good  agree¬ 
ment  with  Gortler' s  curve  up  to  a  (y/  x-x  )  value  of  0.  17 

o 

(a  u/  u  value  of  0.  2) ;  for  larger  (y/  x-x  )  values  the 
c  o 

analytical  curve  is  significantly  higher  than  the  experimental 
points.  From  this  agreement  it  appears  that  for  the  fully 
developed  region  there  is  little  memory  effect  of  the  upstream 
geometry.  This  is  certainly  to  be  expected  because  of  the  high 
entrainment  rates  of  a  jet  flow  and  the  physical  distance  from 
the  nozzle. 

Positive  values  of  the  angle  of  the  mean  flow  with  respect 
to  the  x-axis  show  a  maximum  of  about  3  degrees  at  x/  a  =  100. 
When  the  angular  distribution  is  compared  with  that  calculated 
from  the  continuity  equation,  the  agreement  is  quite  good. 
Albertson,  et  al.  (12)  also  made  flow  angle  measurements 
and  compared  them  with  computed  values  from  the  continuity 
equation.  The  measurements  from  (12)  show  a  maximum 
angle  of  about  10  degrees;  consequently  this  angular  distri¬ 
bution  differs  from  the  measurements  of  (14)  as  well  as  the 
computed  values  of  ( 12)  and  ( 14)  .  It  is  presumed  that  the 
experimental  technique  of  (12)  is  in  error.  A  dual  toUl 
probe  arrangement  with  a  splitter  plate  separating  the  tubes 
was  used;  in  order  to  obtain  a  null  reading  of  this  probe 
in  a  large  mean  velocity  gradient,  the  probe  would  have  to  be 
shifted  to  a  larger-than-true  angle.  This  could  readily  account 
for  the  presumed  experimental  error  of  about  7  degrees  in  (12). 

Liepmann  and  Laufer  ( 15)  have  made  an  extensive  ex¬ 
perimental  study  of  the  half  jet  which  is  formed  by  extending 
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one  of  the  nozzle  walls  past  the  nozzle  exit  (parallel  to  the  mean 
flow  velocity  at  the  exit)  and  by  terminating  the  other  wall  in  the 
usual  manner  for  a  two-dimensional  jet.  An  aspect  ratio  of  eight 
(60  X  7.  5  inches)  was  employed  for  a  nozzle  exit.  The  downstream 
section  was  unbounded,  normal  to  the  60  inch  dimension,  and  all 
measurements  were  taken  within  an  x/  a  range  of  0  to  4.  92 
nozzle  widths.  They  treated  the  flow  as  two-dimensional  for  the 
region  bounded  by  the  plane  of  maximum  velocity  and  the  atmos¬ 
phere. 

Liepmannand  Laufer  measured  velocity  profiles  and  pre¬ 
sented  them  with  the  analytical  predictions  of  Tollmien,  Gortler, 
and  Reichardt.  For  the  presentation  of  the  mean  velocity  data, 
Liepmann  and  Laufer  chose  the  non-dimensional  variable  tj  =  tr  y/  x 
where  (T  is  a  free  parameter  which  is  chosen  to  provide  the  best 
fit  between  the  data  and  the  analytical  curve  under  consideration. 
This  free  parameter  was  used  to  force  the  analytical  profiles  to 

fit  the  experimental  ones  at  the  location  where  u/  u  =0.  5.  The 

c 

physical  location  of  u/  u^  =  0.  5  also  served  as  the  y  =  0  point 
for  the  co-orUinate  selection.  In  order  to  match  the  experimental 
data,  the  solution  of  Tollmien  requires  that  a  be  equal  to  12.0; 
the  solution  of  Gortler  requires  that  a  be  11.0,  and  the  error 
curve  of  Reichardt  requires  u  to  be  13.  5.  Liepmann  and  Laufer 
conclude  that  the  Gortler  solution  provides  the  closest  approxi¬ 
mation  of  their  experimenul  data.  The  mean  velocity  profiles 
of  Liepmann  and  Laufer  are  not  considered  valid  for  the  two- 
dimensional  jet  case  for  the  following  reasons; 

i)  Their  geometry  is  not  truly  two-dimensional, 
ii)  The  co-ordinate  selection  is  arbitrary. 
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iii)  References  (9),  (11),  and  (12)  indicate  that 
the  profile  is  best  approximated  by  the  error 
integral  curve. 

Liepmann  and  Laufer's  measurements  of  the  turbulent 
fluctuation  velocity  ,  the  correlation  iI7,  the  microscale, 

and  the  turbulence  energy  terms  are  believed  to  be  unique  for  this 
type  of  two-dimensional  free  turbulent  mixing.  From  measure¬ 
ments  of  the  mean  velocity  profiles  and  the  Reynolds  stress, 

-p  ^  they  were  able  to  calculate  the  values  of  the  mixing 
length  and  eddy  viscosity  across  the  mixing  layer.  Even  though 
the  geometry  is  not  exactly  the  same,  these  measurements  and 
computed  quantities  are  of  interest  in  connection  with  the  two- 
dimensional  jet  problem. 

A  comparison  of  the  <  and  1  curves  between  Albertson, 
et  al. ,  and  Liepmann  and  Laufer  reveals  gross  dissimilarities  in 
the  shapes  of  both  sets  of  curves.  These  differences  arise  from 
the  different  geometries  or  from  the  fact  that  the  latter  curves 
(15)  are  from  a  direct  calculation,  whereas  the  former  curves 
(12)  were  obtained  from  a  graphical  integration  of  experimental 
data.  Although  references  ( 15)  and  ( 12)  do  not  agree  on  the 
distributions  for  «  and  1  ,  the  two  sets  of  curves  indicate  that 
€  deviates  less  from  a  constant  value  than  does  l  . 

The  relation  between  the  half  jet  and  the  two-dimensional 
jet  is  not  well  understood.  All  the  measurements  of  Liepmann  and 
Laufer  were  taken  within  a  length  (4.92  =  x/a)  that  would  not  allow 
for  the  closing  of  the  potential  core  in  a  two-dimensional  jet. 
Although  the  self-preserving  velocity  profiles  are  somewhat 
similar  to  those  of  the  two-dimensional  jet,  it  is  difficult  to 
compare  the  non  self-preserving  turbulence  data  because  of  the 
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dissimilar  downstream  establishment  rates.  Another  factor,  with 
an  influence  that  is  hard  to  estimate,  is  the  effect  of  the  high  pro¬ 
duction  of  turbulence  energy  near  the  extended  plate  on  the  turbu¬ 
lence  and  mean  flow  properties  of  the  region  which  was  analyti¬ 
cally  treated  as  a  two-dimensional  jet. 

The  nozzle  used  by  Liepmann  and  Laufer  was  (apparently) 
a  beam  deflection  nozzle.  This  type  of  nozzle  design  provides  for 
a  smooth  acceleration  of  the  fluid  through  the  nozzle  (see  Appen¬ 
dix  A.  ) .  A  circular  nozzle  curvature  was  used  by  Miller  and 
Comings  (9)  and  by  Albertson,  et  al.  ( 12) ;  this  is  known  to 
produce  a  non-uniform  velocity  profile  at  the  location  of  the  abrupt 
change  in  the  radius  of  curvature.  The  ratio  of  radius  of  curva¬ 
ture  to  nozzle  width  was  greater  than  6:1  for  (9),  whereas  for 
(12)  this  ratio  was  approximately  1.  The  scatter  in  the  ex¬ 
perimental  data  (12)  for  the  mean  velocity  distribution  as  op¬ 
posed  to  the  smooth  profiles  of  (9)  may  be  a  result  of  this 
difference  in  the  nozzle  construction.  One  of  the  nozzle  designs 
used  in  (11)  and  that  employed  in  (14)  seriously  violate  this 
requirement  for  a  smooth  acceleration  by  using  a  straight  -  wall 
nozzle  design . 

Olson  and  Miller  ( 16)  present  a  rather  clever  analytical 
development  leading  to  an  expression  which  allows  for  the  pre¬ 
diction  of  the  length  of  the  potential  core  region,  the  center  line 
velocity  decay,  and  the  axial  variation  of  the  y  coordinate  at 
the  u/u^  value  of  0.  5.  The  latter  can  be  used  as  a  measure 
of  the  jet  spreading.  As  presented,  the  analysis  is  for  compres¬ 
sible  flow  and  involves  the  assumption  of  constant  eddy  viscosity 
and  a  Gaussian  mean  velocity  distribution.  These  two  assumptions 
would  be  incompatible  for  incompressible  flow. 
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In  summary: 

i)  The  mean  velocity  profile  for  the  two-dimensional 
jet  can  be  accurately  approximated  by  the  Gaussian 
error  curve  (9).  (11),  and  ( 12) ;  the  Gortler  so¬ 
lution  provides  a  satisfactory  approximation  to  this 
profile  (9),  (11),  and  (14)  . 

ii)  Characteristics  for  the  two-dimensional  jet,  such 
as  mass,  momentum,  and  energy  fluxes  (12), 
center  line  velocity  decay  (4),  (9),  and  (12), 
velocity  vector  angles  (11)  and  (14),  and  the 
static  pressure  distribution  (9)  and  (10)  are 
known, 

iii)  Turbulence  data  near  the  nozzle  exit  (9)  and  (11) 
and  a  more  complete  turbulence  survey  for  large 
x/a  (14)  are  available.  However,  certain  im¬ 
portant  discrepancies  in  the  various  results 

exist  as  indicated  by  the  previous  discussion, 

iv)  Details  of  the  turbulence  structure  as  well  as  the 
eddy  viscosity,  mixing  length,  and  mean  velocity 
distributions  are  known  for  the  half  jet  ( 15) . 

Intersecting  Jets 

Several  references  in  the  general  area  of  fluid  jet  am¬ 
plification  are  of  direct  interest  to  the  present  study.  A  general 
description  of  these  is  given  below.  Specific  items  to  be  refer¬ 
enced  from  each  will  be  treated  in  later  sections  wherein  results 
of  the  present  study  are  discussed. 
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Manion  and  Goto  ( 17)  have  analytically  and  experi¬ 
mentally  investigated  the  jet  deflection  problem  for  a  "defined 
region.  "  The  study  was  initiated  as  a  result  of  indications  from 
another  source  that  the  presence  of  the  "setback"  and  "standoff" 
walls  induced  a  greater  deflection  angle  than  that  predicted  from 
momentum  flux  considerations  alone.  From  the  input  data  of  the 
power  and  control  jet  momentum  fluxes,  the  pressure  in  tho  far 
side  pocket,  and  geometric  constants  of  the  defined  region,  the 
analysis  of  Manion  and  Goto  provides  a  prediction  of  the  an¬ 
gular  deflection  of  the  resultant  jet.  A  key  step  in  the  analysis 
is  the  prediction  of  the  near  side  pocket  pressure  by  assuming 
an  inelastic  collision  between  the  power  and  control  jets  and  an 
absence  of  any  shear  -  induced  entrainment  effects  of  the  control 
jet  by  the  power  jet.  The  experimental  program  was  conducted 
on  an  air  flow  system  of  small  physical  dimensions  (a  *  0.  032", 
aspect  ratio  *  8).  The  flow  regime  was  compressible;  power 
jet  plenum  pressures  were  5,  10,  and  20  psig  for  the  tests 

reported.  The  deflection  angle  o  was  determined  by  a  null 
reading  from  a  dual  total  pressure  probe.  The  agreement  be¬ 
tween  the  analysis  and  the  measured  values  of  tan  a  is  sat¬ 
isfactory  if  the  correct  value  of  the  arbitrary  far  side  pocket 
pressure  is  used. 

An  analytical  investigation  of  the  jet  interaction  prob¬ 
lem  in  a  defined  region  was  also  conducted  by  Goto  ( 18) .  Two 
simultaneous  control  jet  flows  are  considered  and  as  in  ( 17) 
a  prediction  equation  for  the  angle  a  was  the  main  objective 
of  the  study.  The  computation  equation  for  tangent  a  involves 
the  geometric  constants,  input  momentum  fluxes,  and  the  pres¬ 
sure  in  one  of  the  pockets.  This  pressure  is,  in  turn,  given 


as  an  independent  function  of  the  input  momentum  fluxes,  the 
geometric  constants,  and  the  deflection  angle  a.  Consequently 
a  trial  and  error  solution  is  implied. 

Reilly  and  Kallevig  (19)  and  Reilly  and  Moynihan 
(20)  have  presented  some  of  the  work  on  the  more  common 
type  of  fluid  jet  amplifier  (i.  e.  X  =  0  and  Y  ^  0)  which  has 
been  carried  on  at  the  Military  Products  Research  Department 
of  Honeywell  Incorporated.  These  presentations  deal  with  the 
gain  of  the  amplifier,  a  prediction  of  the  gain  by  different  flow 
models  and  corresponding  experimental  investigations.  Of 
particular  interest  to  the  present  study  are  the  velocity  pro¬ 
files  taken  in  the  z  =  0  plane  of  a  bounded,  incompressible, 
three  -  jet  flow  system  (aspect  ratio  =  6).  The  profiles, 
shown  in  ( 19),  indicate  that  the  resultartt  profile  rapidly 
approaches  a  symmetric  -(but  not  self  preserving)  shape 
(about  the  x  axis) .  They  also  Indicate  that  a  smooth, 
two-dimensional  like,  velocity  profile  is  not  attained  by  10 
nozzle  widths. 

Olson  (22)  has  taken  jet  deflection  data  for  a  power 
and  control  jet  geometry  similar  to  that  of  the  "defined  region" 
with  the  exception  that  the  "standoff"  and  "setback"  walls 
are  absent.  An  interesting  result  of  this  investigation  is  that 
the  apparent  impingement  point  (apparent  pivot  point)  has  been 
displaced  downstream  of  the  geometric  intersection  point. 

This  displacement  appears  to  be  correlated  with  the  momentum 
flux  ratio  (mfr),  having  a  value  of  8a  for  an  mfr  of  0.016 
and  decreasing  in  an  exponential  fashion  to  a  value  of  0.  4a  at 
an  mfr  of  0.  1.  Olson  attributes  this  displaced  effective 
origin  to  the  elapsed  time  between  the  physical  impingement 
and  the  effective  momentum  interchange.  The  velocity  profiles. 
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which  are  apparently  taken  from  total  pressure  data  only,  show 
good  symmetry  and  smoothness  by  eight  nozzle  widths  downstream. 
These  profiles,  as  well  as  those  from  (19),  indicate  that  if  a 
serious  velocity  discontinuity  is  present  at  the  intersection  point, 
it  is  rapidly  smoothed  by  turbulent  mixing.  The  data  from  (22) 
was  taken  for  a  pow  er  jet  exit  Mach  number  of  0.  66. 

The  majority  of  the  paper  by  Pepperone,  et  al.  (21) 
is  not  of  direct  interest  to  the  present  study;  however  several 
items  from  the  second  section,  "Jet  Streams,  "  are  of  im¬ 
portance.  With  reference  to  an  unpublished  investigation  of 
velocity  profiles  in  the  z  =  0  plane  of  a  bounded  (aspect  ratio  e  8) 
jet  flow,  they  note  that  the  momentum  flux  indicated  by  these 
profiles  is  not  constant  but  falls  off  in  the  downstream  direction 
as  a  result  of  the  shear  effects  on  the  j<  by  the  bounding  plates. 
They  further  note  that  this  would  make  the  jet  appear  to  emanate 
from  a  point  farther  upstream  than  for  the  regular  two-dimensional 
jet  type  flow.  Although  no  indication  is  given  as  to  how  it  was 
defined,  they  state  that  the  apparent  origin  is  four  nozzle  widths 
upstream  of  the  physical  nozzle  exit.  These  data  were  apparently 
collected  with  a  total  pressure  probe  only  and  for  the  compressible 
flow  case  of  a  power  jet  plenum  pressure  of  10  psig. 


OBJECTIVES 


General  Objectives 

The  intent  of  the  present  study  was  to  establish  certain 
important  features  of  jet  flows  between  parallel  plates,  a  com¬ 
plex  flow  phenomenon.  A  complete  and  fully  dcUiled  investi¬ 
gation  which  considered  all  the  parameters  would  require  an 
inordinate  amount  of  time.  Therefore,  specific  topics  were 
chosen  and  tests  (ii)  to  (v)  were  conducted  to  establish 
'information  bench  marks"  from  which  a  much  wider  infor¬ 
mation  area  could  be  surmised.  The  five  specific  subject 
areas  considered  include 

i)  the  single  bounded  jet; 

ii)  the  deflection  angle  (a)  of  the  resultant  jet 

as  a  function  of  nozzle  geometry  and  the  momen¬ 
tum  flux  ratio; 

iii)  the  mean  flow  characteristics  of  the  resultant 
jet  as  a  function  of  the  three  co-ordinates  (x,  y,  z) 
including  flow  angles,  mean  velocity  and  static 
pressure  distribution,  and  mass,  energy,  and 
momentum  flux  ratios; 

iv)  the  effect  of  nozzle  geometry  on  the  resultant 
jet; 

v)  the  effect  of  the  mfr  on  the  resultant  jet. 

Although  the  tests  will  be  discussed  separately  their 
common  intent  is  to  obtain  information  that  will  allow  the  results 
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for  a  specific  set  of  conditions  to  be  extrapolated  to  other  geometric 
and  flow  conditions. 

To  obtain  a  complete  solution  of  the  jet  interaction  problem 
for  the  geometry  chosen  would  require  that  the  turbulent  velocity 
V(x,  y,  z,  t)  be  known  as  a  function  of  the  variables  X,  Y,  and  the 
momentum  flux  ratio  (mfr).  Since  the  problem  is  inherently  one 
of  a  three-dimensional  turbulent  flow,  a  complete  analytical  form¬ 
ulation  involves  ten  unknowns  and  four  equations,  rendering  a 
complete  analytical  solution  impossible.  The  mode  of  attack  on 
the  problem  is  then  governed  by  what  information  is  desired;  this 
information  is,  of  necessity,  short  of  the  complete  solution.  The 
five  specific  investigations  and  their  function  in  the  overall  scheme 
of  the  present  study  are  discussed  below. 

Objectives  of  Several  Experiments 

Although  a  desirable  (and  necessary)  analytical  simplifi¬ 
cation  can  be  effected  by  the  assumption  of  two-dimensional  flow, 
the  flow  regime  of  the  majority  of  fluid  jet  amplifiers  will  be  three- 
dimensional.  This  is  true  because  of  the  aspect  ratios  o'  6  to  8 
which  are  commonly  used.  The  relative  importance  of  the  three- 
dimensional  effects  is  dependent  on  the  characteristic  being 
investigated  and  the  downstream  distance  for  any  given  geometry. 

The  objectives  of  the  first  series  of  tests  was  to  establish  infor¬ 
mation  concerning  the  three-dimensional  character  of  a  single 
bounded  jet. 

The  operation  of  a  proportional  amplifier  is  strongly  de¬ 
pendent  on  the  physical  displacement  of  the  resultant  jet  following 
the  interaction.  If  the  resultant  jet  follows  a  linear  path  from  the 
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interaction,  then  the  deflection  angle  a  and  the  apparent  pivot 
point  are  sufficient  to  specify  the  physical  position  of  the  resul¬ 
tant  jet.  As  noted  in  (17)  and  (18),  the  near  and  far  side  pock¬ 
et  pressures  will  have  a  significant  effect  on  the  deflection  of 
the  jet.  A  series  of  tests  was  conducted  to  determine  the  angle 
a  and  the  apparent  pivot  point  as  a  function  of  the  geometric 
variables,  X  and  Y,  and  the  flow  variable,  the  mfr.  Standoff 
wall  static  pressures  were  also  recorded  for  the  near  and  far 
pockets. 

If  the  angle  a  and  the  apparent  pivot  point  are  known 
for  known  input  momentum  fluxes,  static  pressure  distribu¬ 
tions,  and  the  geometry,  the  momentum  flux  of  the  resultant 
jet  can  be  deduced  by  a  control  volume  analysis  without  know¬ 
ledge  of  the  profile  shape.  However,  the  mass  or  energy  flux 
values  of  the  resultant  jet  cannot  be  so  determined  inasmuch 
as  they  are  dependent  upon  the  shape  of  the  velocity  profile 
of  the  resultant  jet.  For  this  reason,  and  because  they  would 
perhaps  provide  information  on  the  interaction  process,  a 
rather  complete  aeries  of  velocity  profiles  were  taken  at 
x/a  values  of  5,  15.  and  30;  z/a  =  0,  i  1,  *  2.  and 

t  2  —  ;  and  for  X  =  Y  *  2a  with  an  mfr  value  of  0.  1. 

There  are  two  obvious  restrictions  on  the  generality  of  the 
results  from  this  detailed  velocity  profile  survey:  these 
are  the  single  value  of  the  geometric  conditions  and  the 
single  value  of  the  mfr. 

The  fourth  specific  teat  conducted  was  to  determine 
the  effect  of  nozzle  geometry  on  a  single  bounded  jet  and  to 
proWde  information  that  could  be  used  to  extrapolate  the  re¬ 
sults  of  the  detailed  velocity  traverses  to  other  geometric 
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conditions.  The  purpose  of  the  fifth  test,  which  was  similar 
in  intent  to  the  fourth,  was  to  provide  information  on  the  re¬ 
lationship  between  the  profile  characteristics  and  the  value 
of  the  Rifr. 
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EQUIPMENT,  MEASURING  DEVICES, 

AND  INSTRUMENTATION 

Preliminary  Flow  Study 

An  early  decision,  in  the  overall  planning  for  this  study, 
was  to  conduct  a  preliminary  experimental  program  to  determine 
the  magnitude  of  the  three-dimensional  effects  present  in  a  single 
jet  flow  between  parallel  plates.  A  general  description  of  the  ex¬ 
perimental  flow  system  for  this  preliminary  study  is  given  below 
with  special  emphasis  on  the  construction  characteristics  which 
had  an  apparent  effect  on  the  data. 

The  output  of  a  small  centrifugal  blower  was  passed 
through  a  diffuser  to  a  plenum  chamber;  the  chamber,  which 
measured  3x12  inches  in  cross-section,  contained  a  honey¬ 
comb  straightener  and  screens,  A  nozzle  section,  of  the  beam 
deflection  design,  ♦  measuring  3  x  0.  5  inches  at  the  exit, 
was  placed  downstream  of  the  plenum  chamber.  The  nozzle 
blbcks  were  cut  from  a  solid. piece  of  wood  and  then  sanded 
smooth.  The  flow  at  the  exit  of  the  nozzle  was  bounded  by 
parallel  plates  normal  to  the  3  inch  dimension.  With  this  con¬ 
figuration,  the  flow  region  was  open  to  the  atmosphere  on  three 
sideit.  AUhough  the  plexiglass  faced  upper  plate  was  one  contin¬ 
uous  piece  from  the  plenum  chamber  forward,  the  sheet  metal 


*  Appendix  A  provides  an  explanation  of  why  this  de¬ 
sign  was  chosen. 
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covered  lower  plate  was,  for  construction  purposes,  continuous 
only  from  the  nozzle  exit. 

The  construction  of  the  travers#*  device,  which  provided 
three  translational  degrees  of  freedom,  allowed  the  total  and 
static  probes  to  survey  the  upper  two-thirds  of  the  3-inch  space 
between  the  plates  for  any  x,  y  location. 

The  velocity  profiles,  u(z),  showed  a  high  degree  of 
asymmetry  with  respect  to  the  y  and  z  axes.  These  results 
indicated  that  either  the  flow  for  this  geometry  was  naturally 
asymmetric  or  that  the  flow  was  highly  sensitive  to  irregularities 
in  the  surface  geometry.  In  order  to  separate  these  two  possible 
causes,  cerUin  components  of  the  final  flow  system  were  assem- 
beled  to  allow  for  a  study  of  the  single  bounded  jet. 


Single  Bounded  Jet  Study 
Flow  System 

Figure  2  is  a  schematic  drawing  of  the  flow  system  used 
for  the  single  jet  study.  The  flow  bleed  ports  at  the  exit  of  the 
diffuser  serve  three  purposes.  They  provide 

i)  flow  rate  control  through  the  nozzle, 

ii)  stabilization  of  the  diffuser  flow, 

iii)  attainment  of  stable  operating  conditions  for  the 
centrifugal  blower  by  allowing  a  high  enough  flow 
rate  through  the  blower. 

The  grid  at  the  inlet  of  the  plenum  chamber  increased 
the  mixing  of  the  diffuser  efflux  flow  with  the  fluid  in  the  plenum. 
This  mixing  was  necessitated  by  the  cross-sectional  area  dif¬ 
ference  between  the  plenum  and  diffuser  exit.  The  contraction, 
from  the  plenum  to  the  egg  crate  straightener,  and  the  nozzle 


BLEED  PORTS 
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shape  both  followed  the  beam  deflection  curve.  The  honey¬ 
comb  was  3  inches  deep  with  I  -  inch  square  openings. 

Three  40-mesh  screens,  spaced  1-^  -  inches  apart,  were 
placed  downstream  of  the  egg  crate  straightener. 

The  nozzle  blocks  were  constructed  from  eight  pieces 
of  —  -  inch  plywood  which  were  cut  simulUneously  to  follow 
the  bejm  deflection  curve.  After  sanding  for  smoo^thness,  four 
of  the  pieces  were  used  with  four  spacing  pieces  {-  -  inch)  to 
form  each  nozzle  block.  A  plexiglass  sheet,  0.  0  15  inch 
thickness,  was  then  fitted  and  cemented  to  the  plywood. 

Since  both  nozzle  blocks  came  from  the  same  original 
cutting,  the  nozzle  was  symmetric  in  its  final  form.  The 
contraction  ratio  of  the  nozzle  itself  was  24;  the  contraction 
ratio  from  the  plenum  chamber  to  the  nozzle  ew  v/as  294. 
After  the  flow  system  was  assembled,  it  was  discovered  that 
the  vertical  faces  of  the  nozzle  blocks  at  the  exit  were  not 
parallel:  this  resulted  in  a  3  per  cent  variation  in  the  nozzle 
width.  The  plates  which  formed  the  bounding  surfaces  for 
the  jet  were  5/  8"  flakeboard  with  a  1/  16"  Formica  facing. 
Three  2-inch  angle  pieces  of  0.  25  inch  thickness  supported 
the  upper  plate.  For  the  region  in  the  neighborhood  of  the  jet. 

the  plate  separation  was  equal  to  six  inches  +  1-^  percent. 


Measuring  Probes,  Instrumentation,  and  Traverse  Devic^ 

Number  22  hypodermic  needles  (0.025  inch  O.  D. , 
0.0135  inch  I.  D. ) ,  ground  to  provide  a  30  degree  conical 

tip  (included  angle) ,  were  used  as  total  head  tubes.  This 

conical  tip  provided  a  constant  value  of  the  impact  pressure 
over  a  rather  large  range  of  yaw;  it  also  guaranteed  that  the 
displacement  effect  on  the  toUl  head  reading  would  be  negligible 


36 


even  for  the  high  velocity  gradients  encountered  in  boundary 
layer  and  jet  flows.  The  response  of  this  probe  with  respect 
to  yaw  is  given  in  Appendix  B,  Static  pressures  were  mea* 
sured  with  a  closed  end,  twelve-hole,  cylindrical  probe.  This 
probe,  proposed  by  Fage  (23)  in  1936,  has  a  known  yaw  and 
turbulence  intensity  correction.*  A  special  probe  was  con¬ 
structed  to  provide  an  accurate  location  of  the  maximum  velocity 
position  in  the  flow  field.  This  probe,  which  employs  two  total 
head  tubes  with  a  fixed  spacing  of  3/  16  inches,  is  shown  in 
Figure  3.  The  probe  was  mounted  with  the  total  and  static 
probes  on  the  variable  z  traverse  device  (discussed  below) . 

A  micromanometer  was  used  to  measure  the  total  and  static 
pressures;  the  instrument  has  a  minimum  reading  of  0.001 
inch  of  water. 

The  longitudinal  velocity  fluctuation,  11^  was  measured 
with  a  Flow  Corporation  constant  current  hot-wire  anemometer 
using  a  tungsten  wire  of  0.00015  inch  diameter  and  0.040 
inch  unplated  length.  The  signal  from  the  amplifier  output 
was  observed  on  a  Hewlett  -  Packard  oscilloscope.  Model 
120A,  and  quantitative  measurements  were  taken  from  the 
Model  12A1  random  signal  voltmeter.  High  frequency  noise 
was  elminated  at  the  amplifier  output  by  a  7  kc  low  pass  filter. 


*  The  factor  0.  25  given  by  Fage  for  the  correction 
coefficient  K,  where  K  =  (probe  reading  true-static  pressure) 
/p[v*  +'^  ]  has  been  recalculated  by  the  writer  for  the  known 
turbulence  structure  of  the  flow  between  parallel  plates.  The 
value  for  K,  used  in  the  present  work,  was  0.  147.  The  yaw 
correction  given  by  Fage  was  for  a  rather  limited  velocity 
range  (30  -  40  fps).  A  new  formulation  for  yaw  corrections 
was  deduced  by  the  writer  and  is  presented  in  Appendix  B. 


In  order  to  survey  the  entire  flow  field  between  the  bounding 
plates,  a  traverse  device  with  three  translational  degrees  of  free¬ 
dom  was  designed  and  constructed  (x  -  y  -  z  device).  This  device 
is  shown  in  Figure  3  with  the  total  head,  static  head,  and  max¬ 
imum  velocity  detector  probes  mounted  for  a  data  traverse  run. 
The  vertical  motion  is  accomplished  by  a  gear  train.  A  pair  of 
bevel  gears  change  the  operator-applied  rotation  in  the  x  -  z 
plane  into  a  rotation  of  two  spur  gears  in  the  x  -  y  plane.  A 
threaded  rod,  which  is  fastened  to  the  second  spur  gear,  drives 
the  triangular  drill  rods.  Thus,  the  vertical  position  of  the  probe 
tip  can  be  directly  related  to  the  angular  position  of  the  external 
drive  handle.  Accurate  probe  positioning  was  possible  since  a 
vertical  motion  of  0.05  inches  corresponds  to  an  angular  ro¬ 
tation  of  361  least  count  divisions  on  the  graduated  dial  used 
for  the  exterior  measurement  of  rotation.  The  initial  z  position 
of  each  probe  was  set  an  exact  distance  from  the  plate  by  the 
following  method. 

i)  Leads  from  an  ohmmeter  were  attached  to  a  steel 
scale  (24  •  inch)  and  to  the  base  of  the  probe 
(total  or  static) . 

ii)  The  vertical  traverse  carriage  was  then  moved 
toward  the  plate  until  the  probe  tip  made  firm 
electrical  contact  with  the  steel  scale  which  was 
held  on  the  plate. 

iii)  Since  the  steel  scale  was  of  known  thickness,  the 
vertical  position  of  the  probe  tip  was  accurately 
fixed  by  advancing  the  vertical  carriage  away 
from  the  plate  and  noting  where  electrical  con- 
Uct  was  first  broken.  This  latter  technique 
eliminated  gear  train  backlash  error. 
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Figure  3  Total,  static,  and  dual  probes  mounted  on  the  x  -  y  -  z  traverse  device 


39 


Interaecting  Jeta 

Flow  System 

The  left  control  jet  nozzle  and  the  right  side  dummy 
nozzle  shown  in  Figure  4  were  the  only  additions  which 
were  required  to  change  the  single  jet  to  the  intersecting 
jet  flow  system.  Both  nozzles  were  constructed  with  seven 
static  pressure  taps  (0.04  -  inch  diameter)  set  on  0.  125 
inch  centers  and  located  in  the  z  *  0  plane  of  the  standoff 
wall  face.  The  control  jet  nozzle  widths  were  fixed  at  1 
inch.  Values  of  X  =  2  with  Y  =  1,  2,  3  were  employed 
for  the  present  study  (X  and  Y  are  shown  in  Figure  1). 

A  thin  plexiglas  sheet  was  again  used  for  the  nozzle  surface. 
The  air  supply  for  the  control  jet  flow  was  obtained  from  a 
separate  blower  and  plenum  system. 

Measuring  Probes,  Instrumentation  and  Traverse  Device 

The  total  head,  static  head,  and  maximum  velocity 
detector  probes  used  for  the  single  jet  study  were  also  used 
for  the  intersecting  jets.  Because  of  its  insensitivity  to  yaw 
and  its  limited  upstream  disturbance  effect,  a  plate  static 
probe  was  used  in  the  pockets  formed  by  the  nozzle  geometry 
of  the  interaction  region. 

An  important  aspect  of  the  intersecting  jets  problem 
is  the  magnitude  of  the  deflection  of  the  mean  flow  leaving 
the  interaction  region.  In  addition,  to  gain  accurate  data 
surveys,  the  local  velocity  vector  must  also  be  known.  These 
two  requirements  were  met  with  a  direction-finding  probe. 
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FIG.  4  MODIFIED  FLOW  SYSTEM  FOR  THE 
INTERSECTING  JET  FUDW  STUDY 


A  three  tube  probe  (trident)  shewn  on  the  x*y-6  traverse 
device  in  Figure  5,  was  used  as  the  direction  finder.  A 
null-reading  between  the  upper  and  lower  tubes,  which  were 
ground  to  opposing  45  degree  angles,  gave  the  mean  flow 
direction.  The  center  tube  has  a  circular  opening  and  re¬ 
sponded  to  the  total  pressure.  The  sensitivity  of  this  de¬ 
vice  was  dependent  upon  the  magnitude  of  the  dynamic  pres¬ 
sure  and  varied  from  +  0.05  to  approximately  +  0.  4  degrees 
throughout  the  entire  flow  field  of  interest.  This  probe  was 
also  used  with  a  second  positioning  technique.  For  a  fixed 
geometric  orientation,  a  differential  readings  was  used  to 
indicate  the  angular  deflection  of  the  flow  for  a  known  value 
of  the  dynamic  head.  Appendix  B  describes  the  calibration 
technique  and  the  calibration  formula  that  was  deduced  for 
the  fixed  position  method. 

Since  the  flow  leaving  the  interaction  region  was 
not  directed  along  the  x-axis,  a  third  geometric  variable, 

0(  ©  =  arc  tan  v/  u) ,  was  required  to  specify  the  flow  pro¬ 
perties.  Since  various  deflection  angles  would  be  encount¬ 
ered  in  this  study,  a  traverse  device  was  designed  and 
constructed  that  would  provide  two  translational  and  one 
rotational  degrees  of  freedom  (x-y-©  device) .  The  tra¬ 
verse  device  is  shown  in  Figure  5  with  the  trident  probe 
in  the  operating  position.  An  important  characteristic  of 
this  device  is  that  the  three  variables  x,  y,  and  ©  are  in¬ 
dependent  since  the  rotating  element  of  the  traverse  device 
is  a  circular  arc  with  the  probe  measuring  station  fixed  at 
its  center.  (The  variables  x,  y,  ©  are  not  independent  if 
the  axis  of  rotation  is  at  the  base  of  the  probe. )  As  in  the 
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igure  5  Trident  probe  mounted  on  x  -  y-  O  traverse  device 
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operation  of  the  x,  y,  z  traverse  device,  an  exterior  rota¬ 
tion  by  the  operator  in  the  x  -  z  plane  is  converted  to  an 
X  -  y  rotation  of  two  spur  gears,  the  second  of  which  is 
mounted  on  the  arc  holding  the  probe.  The  exterior  rotation 
is  measured  and  converted  into  a  measurement  of  the 
angular  position  (  0)  of  the  probe. 

Both  the  traverse  devices  were  checked  for  possible 
upstream  disturbance  effects.  Neither  device  caused  a 
detectable  upstream  effect  on  the  total  pressure  as  proved 
by  the  negative  results  for  the  following  test: 

i)  For  the  carriage  in  various  positions  in  the  flow 
field,  an  obstruction  was  placed  at  the  base  of 
the  probe  (a  2  inch  x  2  inch  steel  block  was  used) , 

ii)  The  reading  of  a  static  pressure  probe,  mounted 
on  the  probe  holding  device  in  the  normal  manner, 
was  recorded  with  and  without  the  obstruction  in 
place. 

iii)  The  procedure  of  (ii)  was  repeated  using  the 
total  head  probe. 

The  static  pressure  was  unaffected  by  the  x  -  y  -  z 
carriage;  however,  the  x  -  y  -  ©device  seriously  disturbed 
the  static  pressure  readings.  Consequently,  all  static 
pressure  measurements  were  made  with  the  x  -  y  -  z 
carriage  and  yaw  corrections  were  applied  to  the  static 
pressure  readings  so  obtained.  The  correction  technique 
is  discussed  in  detail  in  a  later  section. 
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SINGLE  BOUNDED  JET  STUDY 
Resulta 

If  the  mean  velocity  distribution  for  a  jet  flow  is  given  in 
a  non-dimensional  form,  two  further  variables  must  be  specified 
in  order  to  fix  the  actual  magnitude  of  the  velocity  at  a  given  point 
in  the  flow  field.  These  two  variables,  the  center  line  velocity 
and  a  quantity  which  can  be  used  as  a  measure  of  the  jet  width 
(e.  g.  b) ,  are  referred  to  as  the  principal  characteristics  of  the 
jet.  The  apparent  origin,  the  length  of  the  potential  core  region, 
and  mass  flux  and  energy  flux  widths  are  additional  jet  characteristics. 

Figure  6*  shows  the  values  of  the  momentum  flux  thick¬ 
ness,  b,  for  z/a  s  0  from  the  present  study  as  compared  with  b 
for  a  two-dimensional  jet  (Reference  9) .  Table  3**  indicates 
the  values  of  m,  b,  and  e  at  z/ a  s  0  and  +2  from  the  present  study. 
Mass,  momentum  and  energy  flux  ratios,  viz.,  (M/M^),  Zbq^/aq^, 
and  (E/  E^) ,  are  presented  in  Figure  7  along  with  daU  from  ( 12) 

for  a  two-dimensional  jet.  Figure  8  presents  values  of  u  /  u  from 

m  0 

References  9  and  12  and  the  values  from  the  present  study.  Also 

shown  on  the  plot  is  the  value  of  u  /  u„  which  would  be  obtained 

if  the  velocity  measurement  (of  the  present  study)  were  taken  from 

a  total  pressure  probe  only.  Values  of  u  /  u.  for  z/a  =  0  and  +2 

c  0 

are  presented  in  Table  2. 


*  All  figures  presenting  single  bounded  jet  study  results 
may  be  found  at  the  end  of  this  section,  see  page  70  ff. 

**  The  tables  are  presented  in  the  section  "Discussion  of 
Retiults"  page  47  ff. 
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The  reduction  of  the  total  and  static  pressure  data  for 

the  computation  of  m,  b  and  e  as  well  as  the  computations 

foru(y),  v(y),  C  (y)  and  V(y)  are  discussed  in  section 
P 

VIII  A.  This  discussion  is  presented  in  the  later  section  to 

provide  simultaneous  discussion  of  similar  measurement  and 

data  reduction  techniques. 

Velocity  profiles  that  span  the  full  o"  between  the 

parallel  plates  are  presented  in  Figure  9.  These  u{z/a)/u 

1 

profiles  were  taken  at  x/a  =  15  and  y/a  values  of  0,  t-, 

1  1  ^ 

^  and  +2  —  .  The  data  used  in  the  computa  tion 

of  these  profiles  were  not  corrected  for  turbulence  intensity 
or  yaw  effects  on  the  total  and  sUtic  probes.  The  magnitude 
of  these  errors  is  estimated  to  be  such  that  the  computed 
velocities  are  approximately  5  per  cent  too  high  at  y/a  s  +2~; 
the  yaw  induced  errors  will  decrease  to  zero  at  y/a  =  0. 

Static  pressure  data,  obtained  from  the  traverses  for 
u(z),  are  presented  as  a  series  of  isobaric  contours  in  Fig¬ 
ure  10.  These  data  are  presented  as  h^(y,  z) ;  the  intent  is 
to  emphasize  the  magnitude  of  the  negative  static  pressures 
present  in  the  jet  core.  To  permit  a  ready  calculation  of  a  non- 

dimensional  pressure  distribution  the  value  of  -ypu^  (in  the  same 

2  m 

units  as  h^)  is  listed  on  the  figure.  A  total  of  420  data  points 
were  used  to  determine  the  contours. 

Velocity  traverses  were  made  at  x/ a  =■  5,  10,  15,  and 
30  for  z/a  s  0  and  x/ a  =  5,  15,  and  30  for  z/a  =  1-2.  The 
n(Tj)  /  u^  data  points  from  the  z/  a  c  0  traverses  are  shown 
in  Figure  11;  for  reference  the  analytical  curve  due  to  Reichardt 
(Gaussian  error  curve)  is  shown  on  the  same  figure;  The  inclusion 
of  the  error  curve  also  provides  for  the  comparison  with  other 
two-dimensional  jet  data.  In  order  to  provide  a  comparison 
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between  z/  a  ~  0  and  +2  which  is  not  clouded  by  non-dimension- 
alization  on  b,  Figures  12  to  14  indicate  the  values  of  u(y/a)/u 
for  z/  a  r  0  and  +2  at  x/  a  values  of  5,  15,  and  30. 


Distributions  of  6 and  C 

r 

as  the  lateral  velocity  profiles. 


are  presented  in  th-?  same  manner 
The  G(y)  distributions  are 


presented  in  Figures  15  to  18.  For  visual  symmetry,  the  values 

of  +  e  are  plotted  for  +y  and  -Gfor  -y.  The  C  (jj)  distributions 

P 

for  x/  a  c  5,  10,  15,  ZO  and  z/  a  =  0  are  compared  with  those 

of  the  two-dimensional  jet  (9)  in  Figures  19  to  22.  C  (y) 

P 

distributions  for  z/a  a  0  and  +2  are  shown  in  Figures  23  to  25. 


A  complete  tabular  listing  of  the  processed  data  from  the 


y  traverses,  as  obtained  from  the  computer  program,  is  pre¬ 
sented  in  Appendix  E, 

Figure  26  shows  ^  data  from  the  three  published 

sources  and  from  the  present  study  for  x/a  =  15  and  z/ a  =  0, 
The  data  from  Heskestad  ( 14)  were  taken  at  an  x/  a  value  of  47; 


this  was  the  closest  data  station  to  x/  a  =  15  presented  in  ( 14) . 

Date  from  (9)  are  presented  for  x/ a  s  10  and  20  (to  bracket  15); 

the  data  from  (13)  were  taken  at  x/a  s  17,  5.  Figure  27  presents 

a  comparison  of  5^  /  u*  for  x/  a  =  15  and  z/  a  =  0  and  +2. 

c 

In  order  to  check  the  accuracy  of  the  hot-wire  equipment, 
a  single  point  check  was  made  prior  to  the  above  noted  measure¬ 
ments  by  measuring  5^  on  the  center  line  of  a  long  smooth  pipe 
(fully-developed  flow) ;  agreement  to  within  2  per  cent  of  pub¬ 
lished  data  (23)  was  found. 
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Discussion  of  Results 

In  the  execution  of  the  experimental  p  -ogram,  the  initial 
effort  was  directed  to  the  establishment  of  detailed  velocity 
traverses  as  a  function  of  z,  at  x/ a  *15,  for  eleven  values  of 
y/ a.  These  profiles,  shown  in  Figure  9,  were  to  provide  indi¬ 
cations  of  the  flow  symmetry  and  the  degree  to  which  the  flow 
appeared  to  possess  a  two-dimensional  character.  The  profiles 
showed  a  definite  bulge  near  the  bounding  plates:  the  maximum 


point  of  the  bulge  moved  toward  the  wall  and  became  more 
pronounced  for  increasing  values  of  y.  The  difference  between 

“z  =  0  “max  large;  for  example,  at 

y/ a  =  2-,  u  /u  .  *  2.  50. 

2  max  z  =  0 


Vortex  Stretching  Hypothesis 
To  explain  this  unexpected  phenomenon,  a  hypothesis 
based  on  the  action  of  stretched  vortex  filaments  in  the  jet 
was  formulated.  This  hypothesis  and  further  tests  which 
support  it  are  considered  to  be  one  of  the  major  contribu¬ 
tions  of  the  present  study. 

A  vortex  filament,  which  is  generated  from  the  u(y) 
velocity  gradient  and  extends  between  the  bounding  plates, 
is  convected  with  the  mean  flow  in  such  a  manner  that  the 
central  portion  is  swept  downstream,  whereas  the  portion 
atUched  to  the  wall  is  stationary.  It  is  hypothesized  that  the 
resultant  stretching  and  re-orientation  are  responsible  for 
the  u(z)  velocity  profiles,  departures  from  two-dimensional 
jet  characteristics,  and  other  significant  effects. 
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The  hypothesis  is  based  on  the  following  physical 
reasoning.  With  reference  to  Fipure  28,  assume  that  a 
vortex  filament,  which  extends  from  plate  to  plate,  is 
generated  at  the  location  A  at  time  t.  Since  the  vortex 
filament  is,  by  definition,  formed  from  a  given  collection 
of  particles,  the  vortex  tube  will  be  convected  downstream 
by  the  mean  flow.  The  particles  adjacent  to  the  bounding 
plates  will  remain  stationary  because  of  the  no  slip  con¬ 
dition  at  the  wall;  the  particles  in  the  neighborhood  of  z  =  0 
will  not  be  constrained  and  will  therefore  be  swept  down¬ 
stream  with  the  mean  flow.  Figure  28  also  shows  the  assumed 
position  of  the  vortex  filament  at  time  t  +  At.  The  convection 
does  two  things;  it  stretches  the  filament  and  reorients  the 
vorticity  vector.  At  time  t  this  vector  has  a  strong  com¬ 
ponent  only  in  the  z-direction;  at  time  t  ;  At  the  vector  will 
have  a  weak  y  component  ar.d  a  strong  x  and  z  component 
near  the  plates.  In  an  inviscid  fluid,  the  product  of  the  vorticity 
and  the  cross-sectional  area  of  a  vortex  filament  is  constant. 
Therefore,  if  the  filament  is  stretched,  the  cross-sectional 
area  is  decreased  and  the  angular  rotation  of  the  fluid  particles 
will  be  increased.  For  a  viscous  flow  the  effect  described  above 
is  changed  only  in  degree;  the  qualitative  nature  of  the  phenomenon 
is  unchanged. 

The  assumed  reorientation  of  the  vortex  filament  performs 
two  functions.  The  filament  is  stretched,  increasing-the  vorti¬ 
city,  and  a  non-random  pattern  is  established  which  produces  an 
inflow  near  the  z  =  0  plane  and  an  outflow  near  the  bounding  plates. 

Two  further  confirmations  (following  the  u(z)  profiles)  of 
the  above  hypothesis  are  available  from  the  literati  re.  During 


(b.)  VORTEX  FILAMENT  AFTER  TIME  At 


FIG.  28  ASSUMED  VORTEX  PRODUCTION  AND  STRETCHING 
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preliminary  studies,  v.  d.  H.  Zijnen  (11)  attached  two  bounding 
plates  at  the  nozzle  exit  (geometrically  the  same  as  the  present 
study  except  that  an  aspect  ratio  of  20  was  used) .  Measurements 
of  the  total  pressure  at  y  s  0  and  x/ a  s  40  as  a  function  of  z 
show  a  definite  velocity  increase  away  from  the  z  =  0  plane. 

The  maximum  velocity  is  approximately  12  per  cent  higher  than 
Similarly,  Curtet  (25)  has  reported  data  which  is  an  aside 
to  his  principal  investigation.  The  geometry  of  (25)  is  a 
rectangular  jet  exhausting  into  a  rectangular  channel.  The 
system  of  (25)  differed  from  the  present  study  in  that  there  were 
walls  at  +  y  »  5a  (for  the  test  of  interest)  and  a  secondary  flow 
could  be  introduced  parallel  to  the  primary  flow  at  x  s  0,  Also, 
the  aspect  ratio  of  the  nozzle  was  10.  Curtet  made  measurements 
of  the  velocity  in  the  y  =  0  plane  at  x/a  =  200.  The  results  of 
these  measurements  show  a  maximum  velocity  at  z/a  s  4 
(upper  and  lower  plates  at  z/  a  *  5)  which  is  approximately  12 
per  cent  higher  than  u^  (at  the  same  x/a  location) .  Although 
the  geometry  is  somewhat  different  from  the  present  study  (viz., 
the  presence  of  the  side  walls),  the  vortex  stretching  hypothesis 
may  also  be  used  to  explain  the  observations  of  Curtet. 

Although  the  farthest  downstream  check  was  at  x/ a  =  25,  data  from 
the  present  study  indicate  the  same  trend  as  (11)  and  (25).  From  Table  2 

Table  2 

u  /  u-  As  a  Function  of  x/a 
c  0 

x/  a 

z/  a 
+  2 
0 

-  2 


5 

10 

15 

25 

0.995 

0.  870 

0.  716 

0.  552 

0.988 

0.  840 

0.  705 

0.531 

0.997 

0.  843 

0.686 

0.  535 
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Uq  at  z/ a  =  +2  is  approximately  4  per  cent  greater  than 

u  /  u.  for  an  x/  a  s  ?5.  For  z/a  s  -2  u  /  is  only 
m  U  c  0 

slightly  greater  than  u  /u.  at  x/a  =  25;  however,  the 

m  u 

ratio  u  y  ■>/“  at  x/  a  =  25  does  show  a  marked  increase 
z/  as  -4/  m 

with  respect  to  x/  a  =  15.  A  conclusive  statement  about  the 

effect  is  prohibited  by  the  unsymmetric  behavior  of  the  u  /  u„ 

c  0 

values  with  respect  to  z  (probably  caused  by  the  3  per  cent 

deviation  in  the  nozzle  width  over  the  height  of  the  nozzle), 

and  the  variations  of  u  /  u.  with  respect  to  u  /  u,^  as  a 

c  0  m  0 

function  of  x.  These  variations  may  be  indicative  of  a 
significant  flow  phenomenon;  however,  the  data  presented 
are  not  sufficient  to  establish  what  this  might  be. 

These  three  observations,  that  the  velocity  is 
higher  away  from  the  z  s  0  plane,  indicate  an  important 
aspect  of  the  vortex  stretching  hypothesis,  viz. ,  that  the 
phenomenon  is  convective.  Since  there  will  be  a  retarda¬ 
tion  of  the  flow  adjacent  to  the  wall,  and  since  this  will 
manifest  itself  as  a  boundary  layer-like  profile  (near  the  wall) 
the  fluid  near  the  plates  will  be  decelerated.  If  this  effect 
were  to  continue  (in  the  absence  of  other  effects)  a  profile 
would  develop  with  a  maximum  velocity  at  z  =  0,  In  order 
to  overcome  this  effet  t  and  establish  a  higher  velocity  near  the 
plates,  a  convective  artion  must  occur.  This  convective 
action  is  allowed  for  in  the  hypothesis  by  the  inflow  at  z  s  0 
and  the  outflow  near  the  plates  with  the  resulting  flow  from 
the  z  s  0  plane  toward  the  bounding  plates  in  the  neighborhood 
of  y  =  0. 

A  natural  <  <  sequence  of  the  vortex  stretching 
hypothesis  would  be  the  existence  of  larger  6  values  as  well 
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as  larger  mass,  momentum,  and  energy  flux  thicknesses  at 
z/a  B  +2  as  compared  with  zl  a  =  0  for  a  given  x.  Also,  a 
ratio  formed  in  this  comparison  should  have  a  larger  value 
for  larger  x/  a  stations,  e.  g. 


I 


z/  a  =  +2 


e/  a  =  0 


30  ^  *b 


b/  a  =  +2 


e/  a  =  0 


) 


c/  a  =  15 


Figures  15,  17,  and  18  show  the  0(y/ (x-x  ))  values  for 

o 

z/ a  c  0  and  +2,  for  the  x/  a  values  of  5,  15,  and  30.  Figure 

16  shows  0(y/  (x-x  ) )  for  z/a  =  0  and  x/  a  =  10.  The 
o 

measured  0  values  for  a  two-dimensional  jet  (14)  are  also  shown 
on  the  figures.  *  The  predicted  differences  in  the  distribution 


*  The  similarity  variable,  y/ (x-Xq),  was  used  because 
the  data  of  ( 14)  are  presented  with  0  as  a  function  of  y/  (x-Xq)  . 
The  intent  of  using  the  non-dimenionalizing  factor  x-Xq  makes 
a  jet  of  finite  width  appear  to  emanate  from  a  nozzle  of  infini¬ 
tesimal  width.  However,  there  are  two  commonly  used  ways  to 
define  XqJ  they  are:  (i)  the  location  where  the  straight  line 
drawn  through  the  u^./ Uq  points  (for  large  x/a)  crosses  the 
value  1.0  on  Figure  8,  or  (ii)  the  location  where  the  measure-  ' 
of-the-jet-width  is  equal  to  zero  when  extrapolated  upstream 
of  the  core' s  closing  point  (e.  g.  Xq  =  -  1.  572a)  for  (9)  as 
indicated  on  Figure  6.  From  Heskestad' s  data  the  two  defini¬ 
tions  give  the  same  x^  of  6a;  however,  it  should  be  noted  that 
the  velocity  distribution  was  not  measured  closer  to  the  nozzle 
than  x/  a  =  47,  so  that  the  value  of  x^  was  obtained  by  a  very 
long  estrapolation  of  the  linear  curve  passed  through  the  data 
points.  Neither  of  the  two  distributions  follows  a  linear 
relationship  for  the  bounded  jet  nor  do  they  lead  to  the  same 
Xq  (the  latter  item  is  in  agreement  with  the  results  of  (9)). 

For  the  present  study  Xq  was  taken  as  -3a;  it  was  arrived  at 
by  extrapolating  the  line,  which  was  faired  through  the  four 
points  of  the  present  study,  to  the  x-axis.  It  is  felt  that  the 
intent  of  using  the  variable  x-x^  is  better  served  if  the  infini¬ 
tesimal  jet  width  criterion  is  satisfied,  as  opposed  to  the 
maximum  x  location  where  u^  =  u^.  Finally,  the  probleni 
of  choosing  the  variable  x^  offers  a  good  reason  for  using 
a  more  well  defined  non-dimensionalizing  variable  such  as  b. 
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are  not  present  for  x/  a  =  5  and  15;  however,  the  6  distribution 
for  x/ a  s  30  clearly  shows  the  effect  of  the  vortex  stretching. 
The  u(z/a)/u^  profiles  at  x/ a  s  15,  Figure  9,  indicate  that 
the  velocity  bumps  appear  closer  to  the  plate  than  z/a  =  2. 
Therefore,  the  similar  6  distributions  at  x/ a  ^  5  and  15  for 
z/a  =  0  and  +2  do  not  provide  a  rejection  of  the  hypothesis: 
they  merely  indicate  that  the  major  part  of  the  outflow  from 
the  y  r  0  plane  occurs  for  z/  a  greater  than  +2  at  these 
xl  a  =  locations. 

The  values  of  the  three  flux  thicknesses  shown  in 
Table  3  also  support  the  hypothesis. 

Table  3 

Values  of  m,  b,  and  e  for  the  Single  Bounded  Jet 


x/  a 


z/  a 

5 

10 

15 

30 

m 

0 

1.48 

2.  20 

2.92 

5.39 

m. 

+  2 

1.50 

3.05 

6.  76 

b 

0 

0.  534 

0.  786 

1.06 

2.02 

b 

+  2 

0.  564 

-- 

1.  10 

2.  51 

e 

0 

0.  869 

1.  26 

L 

1.  73 

3.  33 

e 

+  2 

0.951 

1.  80 

4.  12 

The  comparison  of  the  ratio' s  defined  above,  indicate  that  the 
predictions  based  on  a  vortex  stretching  model  are  correct 
with  the  exception  of  the  b  and  e  ratio  comparisons  between 
x/  a  =  5  and  15.  As  noted  above,  the  u(z)  profiles  for  x/ a  = 

15  indicate  velocity  bumps  for  |  z/ a  |  greater  than  2  indi¬ 
cating  that  the  convection  effects  of  the  vortex  stretching  are  not 
yet  felt  at  z/  a  r  +2.  Consequently,  the  lack  of  agreement 
between  prediction  and  observation  is  not  taken  as  a  rejection 
of  the  hypothesis. 
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Another  comparison  which  also  supports  the  hypothesis 
is  shown  in  Figures  12  to  14,  viz.  the  u(y/ a)  /  u^  velocity 
profiles  for  z/ a  s  0  and  +2.  As  in  the  two  cases  above,  only 
slight  differences  appear  for  x/a  =  5  and  15;  however,  the 
profiles  at  x/ a  =  30  clearly  indicate  the  effect  of  the  outflow 
from  the  y  =  0  plane  at  z  s  1-2  as  well  as  the  convection  of 
momentum  from  z/a  =  0  to  z/a  s  2. 

The  static  pressure  contours,  shown  in  Figure  10 
were  not  expected  a  priorL  From  an  examination  of  the 
figure  it  appears  that  the  cross  section  may  be  broken  into  three 
regions;  they  are: 

i)  the  central  portion,  approximately  bounded  by 
8P 

the  planes  z  =  +2,  where  the  component  is  the  dominant 
part  of  the  pressure  gradient, 

ii)  a  region  where  the  isobars  are  parallel  to  the 

0jp 

bounding  plates  and  ~  is  not  only  the  dominant  term  of 
o  z 

the  pressure  gradient  but  has  a  rather  large  value, 

iii)  a  transition  region  between  i)  and  ii) . 

Athough  the  following  argument  is  not  amenable  to 
quantitative  proof  at  present,  it  is  advanced  as  an  explanation 
of  the  observed  phenomenon  and  is  based  on  the  vortex 
stretching  hypothesis. 

The  static  pressure  distribution  in  region  (i)  is  of 
the  same  general  shape  as  that  for  a  two-dimensional  jet  (for 
large  x/  a) .  (This  comparison  is  discussed  in  a  later  section.  ) 
The  distribution  in  the  second  region  (ii)  cannot  be  attributed  to 
either  a  jet  or  a  boundary  layer  type  of  flow.  However,  if 
the  convected  flow  due  to  the  vortex  stretching  is  considered 
as  being  directed  toward  the  plate  with  sufficient  velocity,  the 
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large  (and  positive)  8P/8z  pressure  gradient  is  simply  the 
decelerating  agent  for  the  convected  flow.  If  more  complete 
aata  were  available,  the  aP/8z  term  could  be  employed  as 
a  part  of  an  analysis  to  evaluate  the  strength  of  the  convective 
flow. 

Static  pressure  distributions,  presented  as  C  (y) , 

P 

for  z/a  =  0  and  2  at  x/a  =  5,  15,  and  30  are  shown  in 
Figures  23  to  25.  The  important  feature  of  Figures  24  and 
25  is  the  high  velocity  portion  at  z/a  =  +2  wherein  the 
pressure  is  constant  for  a  width  of  (approximately)  +  1  inch 
and  +  2  inches.  As  noted  in  the  discussion  of  the  static 
pressure  contours  for  the  x/a  =  15  location,  the  higher 
pressure  near  the  wall  is  probably  the  decelerating  agent 
for  the  convected  flow  from  the  z  =  0  plane. 

Turbulence  Measurements 
Measurements  of  S*  were  made  in  the  z/a  *  0 
and  +2  planes  with  the  intent  of  providing  data  for  comparison 
with  the  two-dimensional  jet  and  possible  support  data  for  the 
proposed  vortex  stretching  hypothesis.  Figure  26  presents  the 
u  data  which  is  available  for  the  two-dimensional  jet;  data 
from  the  present  study  for  z/a  =  0  is  presented  on  the  same 
figure.  It  is  obvious  from  the  figure  that  there  is  no  general 
agreement  for  the  two-dimensional  jet  case.  There  are 
several  possible  explanations  for  the  differences  among  the 
results  of  (9),  (lo),  and(14).  It  is  unlikely  that  the  Reynolds 
number  differences  will  account  for  the  different  magnitudes  and 
curve  shapes  since  all  tests  were  conducted  in  the  subsonic 
range  and  the  Reynolds  numbers  are  within  the  same  order  of 
magnitude.  It  would  not  seem  likely  that  the  aspect  ratio 
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differences  would  be  significant:  the  only  discrepancy  may  lie 
in  the  work  of  (13)  wherein  the  ratio  of  25  may  be  too  low  to 
obtain  a  two-dimensional  flow.  The  instruments  used  are 
noted  for  reference  puposes.  The  symbols  I  and  T  on  the 
figure  denote  the  type  of  instrument  used  (I  =  consta  denotes 
a  constant  current  device,  T  =  constant  denotes  a  cons‘'ant 
temperature  device) .  The  linearity  of  the  constant  current 
instrument  may  be  in  question  for  high  relative  intensities: 
near  the  outer  portion  of  the  jet  the  local  intensity  approaches 
values  of  50  to  60  per  cent  /  u*  ) .  The  two  most  likely 
factors  that  might  cause  the  observed  differences  are  the 
nozzle  geometry  used  and  the  turbulence  intensity  at  x/ a  =  0, 
The  Sq/uq  factor  was  not  reported  for  (13)  or  (14);  however 
the  Sq/  u^  for  ( 14)  will  probably  be  of  the  same  order  as 
the  value  for  fully  developed  flow  between  parallel  plates.  Of 
the  three  nozzles  used,  the  one  from  (9)  would  appear  to  be 
the  form  most  commonly  used  for  applications. 

With  the  lack  of  a  known  reference  curve  for  a  two- 

dimensional  jet,  it  is  difficult  to  make  a  meaningful  comparison 

with  the  two-dimensional  case;  however,  certain  differences 

are  pertinent.  With  reference  to  the  curves  from  (9)  and  (14), 

it  appears  that  there  is  less  ti.rbulence  energy  in  the  central 

portion  of  the  bounded  jet  than  in  the  two-dimensional  jet. 

Since  the  velocity  gradient,  u(y),  is  similar  to  that  for  the  two- 

dimensional  jet  (in  lact,  it  is  slightly  steeper),  one  would  expect 

at  least  as  much  turbulence  energy  for  the  central  portion  of  the 

bounded  jet.  If  this  is  indeed  the  case,  three  other  effects 

could  cause  the  lower  5^  /  u^  curves  for  0  <  ly/  x|  <  0.  06. 

If  a  convective  action  were  present,  the  lower  in 

m 
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this  region  could  be  the  result  of  a  convection  from  the  z  =  0 

plane  to  the  +  z  planes.  Alternatively,  if  a  non-random  flow 

pattern  existed  with  its  major  flow  components  in  the  y  and/  or 

z  directions,  the  energy  apparently  missing  from  the  A*  term 

could  be  present  in  the  P  and/ or  w^  terms.  Therefore,  if 

the  total  turbulence  energies  for  the  two  cases  were  equivalent, 

this  would  explain  the  lower  values  of  /  u^  These  two 

c 

effects  are  natural  results  of  the  vortex  stretching  hypothesis. 
The  third  possibility  is  an  increase  in  the  dissipation  of  the 
turbulence  energy.  An  energy  balance  in  (14)  indicates  that 
the  dissipation  term  is  fairly  large  in  the  neighborhood  of  y  = 

0.  Therefore,  if  some  action  increased  the  dissipation,  this 
would  also  explain  the  lower  /  u*  .  A  higher  dissipation 
would  result  from  a  larger  fraction  of  the  energy  in  the  small 
eddies;  however,  there  is  no  reason  to  attribute  this  effect  to  the 
vortex  stretching  action.  The  generation  of  small  eddies  by 
the  plates’^  is  not  a  reasonable  explanation  unless  there  is 
a  strong  transport  effect  that  would  carry  them  from  the 
relatively  thin  boundary  layer  up  to  the  z/ a  =  0  plane.  The 
convection  due  to  the  vortex  stretching  will  not  do  this;  this 
convective  action  will  primarily  transport  large  eddies  from 
the  outer  portion  of  the  jet. 

With  regard  to  Figure  27  which  compares  the 


^Although  there  is  no  way  to  compare  diref-tly  free 
shear  and  boundary  layer  type  flows,  the  boundary  layer  flows 
will,  in  general,  have  more  energy  in  their  small  eddies; 
this  is  due  to  the  existence  of  the  very  high  velocity  gradients 
near  the  wall  and  the  small  length  scale  of  the  boundary  layer 
thickness . 


58 


J*  /  values  of  z/  a  =  0  and  +2,  three  regions  are  apparent: 

i)  region  (i) ,  for  0  jc  y/  x  <  0.  07 

ii)  region  (ii) ,  for  0.07  Ky/  x  <  0.  14 

iii)  region  (iii),  for  y/ x  >_  0.  14. 

The  significance  of  (i)  is  that  the  S*  for  z/a  =  +2  is  higher 
than  for  z/ a  s  0.  Since  the  u{y)  gradient  for  z/a  =  +2  is 
less  steep  than  that  of  z/a  =  0.  the  (apparently)  greater 
turbulence  energy  for  z/  a  =  +2  requires  explanation.  A 
convection  effect  from  z/a  *  0  to  1  z/  a  values  for  the 
neighborhood  of  y  =  0  is  predicted  by  the  hypothesis;  there¬ 
fore,  the  turbulence  energy  may  be  convected  toward  the 
boundary  plates  and  this  would  explain  the  greater  at 
z/  a  =  +2.  However,  the  apparent  energy  is  still  much  less 
than  that  for  (9)  and  (14) .  The  energy  generated  in  the  wall 
shear  layer  and  diffused  outward  may  account  for  all  or  part 
of  the  difference  and  this  explanation  would  be  independent 
of  the  hypothesis.  Without  further  measurements  these  two 
effects  cannot  be  accurately  separated.  There  is  however 
no  production  of  turbulence  from  the  -pTRT  stresses,  since 
8u/  dz  s  0  as  seen  in  Figure  9. 

Region  (ii)  is  undoubtedly  a  direct  result  of  the 
steeper  velocity  gradient  at  z/a  s  0  leading  to  a  higher  pro¬ 
duction  of  turbulence  energy. 

Region  (iii)  is  interesting  but  nebulous.  The  jet  is 
physically  wider,  and  this  alone  may  account  for  the  difference 
in  /  u^.  Also,  production  of  turbulence  energy  will  result 
from  the  term  aw  8u/  8z.  The  much  higher  intermittency 
factor  for  z/ a  =  t-2  means  that  the  relative  intensity  of  tur¬ 
bulence  in  the  turbulent  fluid  at  z/  a  s  +2  is  much  lower  than 
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for  z/  a  r  +2.  This  result  is  incompatible  with  the  vortex 
stretching  hypothesis  as  stated  because  the  convective  action 
of  the  secondary  flow  should  introduce  high  turbulence  in¬ 
tensity  fluid  for  large  |z|  (from  the  central  region)  and 
low  turbulence  intensity  fluid  for  an  all  |z|  (which  has 
been  entrained  from  the  atmosphere).  However,  there  are 
three  possible  explanations  for  this. 

i)  If  the  secondary  flow  occurs  entirely  within 
the  jet  as  bounded  by  the  viscous  super  layer  *  proposed 
by  Corrsin  and  Kistler,  then  the  secondary  flow  will  not 
bring  in  ambient  fluid,  but  will  convect  fluid  from  the 
region  near  the  plates  to  the  region  of  z/a  e  0.  Since  some 
of  this  fluid  will  have  come  from  the  shear  layer  near  the 
plate  it  will  have  a  relatively  high  turbulence  intensity. 

If  this  explanation  describes  the  actual  flow 
conditions,  then  the  explanation  of  a  higher  dissipation 
for  small  |z|  and  small  |y|  (as  compared  with  large  |z| 
and  small  |y()  would  be  supported.  The  original  discussion 
of  this  point  is  on  page  57. 

ii)  Diffusion  in  the  +  y  directions  of  the  relatively 
high  intensity  turbulence  from  the  middle  section  of  the  jet 
at  z/ a  s  0  may  be  strong  enough  to  overcome  the  convective 
action  of  the  vortex  stretching. 


The  viscous  super  layer,  in  the  model  of  {  26  ) , 
is  a  ve  .-y  thin  region  (wherein  viscous  forces  predominate) 
which  is  the  actual  and  instananeous  boundary  of  the  jet.  It 
is  the  contention  of  (26  )  that  the  entrainment  of  ambient 
fluid  occurs  in  the  interaction  between  this  layer  and  the 
ambient  fluid,  that  is,  vorticity  is  transmitted  only  through 
this  interaction. 
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iii)  The  constant  current  device  used  is  also  unreli¬ 
able  for  large  turbulence  intensities  ((6^  /  u  *  30  per  cent) 

and  the  internnittency  factor  was  determined  simply  from  a 
visual  inspection  of  the  oscilloscope.  Both  factors  make  the 
accuracy  of  the  data  somewhat  questionable.  It  should  be 

noted  that  whereas  the  actual  di^'ference  in  /  u^  values  for 

c 

the  two  z  values  is  no  greater  than  for  region  (i) ,  the  value 
at  z/ a  s  +2  is  50  per  cent  higher  than  that  for  z/a  s  0  in  region 
(iii)  .  Although  it  is  also  true  of  the  turbulence  data  in  general, 
no  concrete  statements  may  be  made  about  (iii)  without  more 
information;  however,  for  a  reasonable  understanding  of  the 
jet  flow  between  bounding  plates,  these  and  other  phenomena 
should  be  investigated  further. 
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Comparison  with  the  Two-Dimensional  Jet 

It  should  be  borne  in  mind  throughout  this  discussion 
that  the  existence  of  the  vortex  stretching  action  prohibits  the 
flow  from  assuming  a  truly  two-dimensional  character.  However, 
certain  properties  may  exhibit  behavior  which  are  sufficiently 
two-dimensional  in  nature  to  assist  in  the  understanding  of 
jet  flow  between  parallel  plates  as  well  as  to  provide  certain 
bases  for  prediction  and  design  efforts. 

From  the  literature  survey,  the  Gaussian  error  curve 
provides  an  accurate  representation  of  the  non-dimensional  ve¬ 
locity  profile.  From  Figure  11  it  can  be  seen  that  the  non- 
dimensional  profiles  of  the  present  study  agree  quite  well  with 
those  of  the  two-dimensional  jet.  However,  this  is  only  part  of  a 
meaningful  comparison.  Two  parameters  are  employed  to 
non-dimensionalize  the  profiles,  u^  and  b.  Consequently, 
if  a  comparison  of  the  actual  velocities  were  of  interest,  it 
would  also  be  necessary  to  compare  the  u  /  u  and  b(x) 
distributions.  From  Figures  6  and  8  it  may  be  seen  that  b  is 
less  and  u^  is  greater  for  the  bounded  jet  at  a  given  x  in  the 
region  downstream  of  the  core's  closing.  Considering  these 
two  effects,  it  is  only  fortuitous  that  they  combine  to  yield  a 
non-dimensional  distribution  which  agrees  so  well  with  the 

Gaussian  curve  of  (9) .  These  two  effects,  higher  u  and 

m 

lower  b,  are  compatible  with  the  vortex  stretching  hypothesis. 

Because  of  the  inflow  at  z  *  0  the  jet  should  be  narrower 

(lower  b) ,  and  from  continuity  considerations  the  central  region 

velocity  should  be  higher,  (higher  u  )  .  Another  observation 

m 

with  regard  to  Figure  8  is  that  u^  approaches  the  two- 
dimensional  value  for  large  x/a.  It  is  assumed  that  this  is  . 
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natural  result  of  the  vortex  stretching  convective  effects. 

The  vortex  stretching  convective  action  is  discussed  in 

detail  in  the  previous  section  on  the  vortex  stretching  hypothesis. 

Another  example  of  the  importarce  of  comparing 

more  than  the  non-dimensional  profiles  is  available  from  the 

literature.  Olson  and  Miller  ( 16)  state  that  their  velocity  profiles 

follow  a  Gaussian  distribution.  This  might  lead  one  to  presume 

that  the  jet  flow  between  bounding  plates  (aspect  ratio  12) ,  and 

at  a  Mach  number  at  thep<Mer  jet  exit  of  0.66,  was  equivaient  to 

the  incompressible  two-dimensional  jet  flow.  However,  Figure  8 

shows  the  values  of  u  /  u.  for  two  two-dimensional  iet  flows 
m  0  ■' 

and  the  values  from  Olson  and  Miller  (16);  the  much  higher 
values  from  ( 16)  indicate  that  the  jet  is  spreading  less  rapidly 
(since  the  non-  dimensional  profiles  are  both  Gaussian  distri¬ 
butions)  .  A  measure  of  the  rate  of  spread  may  be  taken  from 
the  angle  that  a  line  connecting  the  u^/  2  velocity  points  makes 
with  the  X  axis.  This  angle  for  Olson' s  data  was  4.  3  degrees; 
it  was  5.  5  degrees  for  the  data  of  Miller  and  Comings. 

The  mean  velocity  profiles  provide  an  example  of  a 
comparison  based  on  the  distribution  of  a  given  quantity.  Other 

comparisons  of  this  type  are  presented  later.  The  u  and  b 

m 

comparisons  are  examples  of  two  comparisons  based  on  the 
characteristics  of  the  jet. 

A  third  comparison  between  tlie  characteristics  of  the 
two  flow  cases  is  offered  by  Figure  7  which  shows  the  mass, 
momentum,  and  energy  flux  ratios.  (Actually,  physically  meaning¬ 
ful  values  of  the  mass,  momentum  and  energy  fluxes  for  the  bounded 
jet  at  a  given  x/  a  location  should  be  determined  as  an  average  of 
measurements  several  z/a  locations.  However,  the  intent 


here  is  to  compare  the  center  plane  flow  with  the  two- 
dimensional  case;  consequently,  the  single  traverse 
evaluation  for  the  ratios  is  a  suitable  one  to  use. )  The 
higher  ratios  for  the  bounded  jet  up  to  x/ a  =  15  indicate  that 
a  greater  entrainment  is  occurring  up  to  this  x  location.  The 
ratio  values  at  x/a  =  30  are  less  than  would  be  expected 
on  the  basis  of  the  values  up  to  x/ a  s  15.  The  lower  values 
at  x/ a  s  30  may  be  explained  by  the  convective  effects  of 
the  vortex  stretching  hypothesis.  Lower  ratio  values  for 
z/  a  =  0  would  occur  if  a  convective  action  swept  the  high 
velocity  fluid  from  the  central  region  of  the  jet.  As  observed 
in  the  earlier  section,  the  convective  action  appears  to  be  much 
stronger  for  x/ a  s  15  to  30  than  for  x/a  =  0  to  15.  This  is 
also  compatible  with  the  hypothesis  in  that  it  appears  to  be 
a  developing  effect.  Also,  in  view  of  the  definitions  of  the 
three  ratios,  the  energy  flux  ratio  should  be  the  most  sensitive 
to  this  loss  in  high  velocity  fluid  and  the  mass  flux  ratio 
should  be  the  least  sensitive.  From  the  figure,  it  is  seen  that 
this  is  indeed  the  case. 

Two  other  items  are  noteworthy.  The  first  of  these 

involves  a  comparison  between  the  three  flux  ratios  for  the 

bounded  jet  as  compared  with  the  two-dimensional  jet. 

Specifically,  the  difference  between  the  E/ values  is 

greater  than  the  difference  between  the  2bq  /  aq  values 

m  ^0 

or  the  M/  values.  The  reason  for  this  is  two-fold; 

(i)  E  was  calculated  in  ( 12)  by  integrating  the  term  u’  ; 

in  the  present  study,  it  was  evaluated  by  integrating  (q)  (u)  (Note 

q  =  y  p  ^  -T  p  u^  )  and  (ii)  the  higher  u  /  u.  as 

shown  in  Figure  8  leads  to  greater  E  values  since  the  integral 

is  dependent  on  the  third  power  of  the  velocity. 
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The  second  itenn  is  the  existence  of  momentum  flux 
ratio  values  which  are  greater  than  1.0.  This  can  be  partly 
explained  by  the  static  pressure  distribution  across  the  jet. 
Consider  a  control  volume  with  surfaces  that  ere  the  nozzle 
exit,  a  downstream  plane  of  constant  x,  two  side  planes  which 
are  placed  such  that  the  entrainnoent  flow  enters  the  control 
volume  at  a  6  of  +  90  degrees,  and  two  planes  of  constant 
z  placed  at  z  s  +Az  and  z  =  -Az  (where  Az  is  small  with 
respect  to  the  |z|  of  the  bounding  plates).  The  x-direction 
momentum  equation  can  then  be  written  as: 
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where  W  is  the  width  of  the  control  volume  and  the  viscous 
shear  is  assumed  to  be  negligibly  small.  It  is  reasonable  to 
assume  that  is  approximately  equal  to  atmospheric 

pressure;  consequently,  if  P  is  considered  to  be  a  sase 

X  =  X  ®  ® 

pressure  the  first  left  hand  side  term  is  zero. 

The  effect  due  to  the  x-momentum  which  is  carried 
in  by  the  (three  dimensional)  velocity  component  w  is  probably 
a  significant  factor;  however,  there  is  no  data  from  which 
this  term  may  be  evaluated. 

The  shear  stress  terms  (on  the  right  side)  are  of 
unknown  importance.  It  is  reasonable  to  assume  that  the 
correlation  term  (KH)  is  small  because  of  the  flat  u(z) 
velocity  profile,  see  Figure  9;  however,  the  area  over  which 
this  shear  stress  acts  is  quite  large,  viz.,  xW. 

Although  four  of  the  eight  terms  of  the  equation  cannot 
be  numerically  specified,  the  equation  can  be  manipulated  to 
yield  the  effect  of  the  static  pressure. 

Dividing  by  4aqQAz  and  rewriting  the  equation  yields 
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for  the  first  two  terms  on  the  left  side  of  the  equation.  (Note  that 

the  first  integral  on  the  left  is  equivalent  to  8bq  Az.  )  This 

m 

quantity  can  be  formed  from  the  sum  of  three  similar  terms  that 
result  from  the  two  independent  sets  of  terms  on  the  right  side 
and  the  last  .wo  terms  on  the  left  side,  viz.,  the  pressure,  the 
shear  stress  and  the  negative  of  the  other  momentum  flux  terms. 
For  example  the  pressure  effect  may  be  solved  for  as 
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In  order  to  obtain  a  quantitative  evaluation  for  the 

magnitude  of  the  pressure  effect,  the  above  integral  was 

evaluated  at  x/ a  =  15  using  the  daU  of  Figure  24  and 

the  value  was  found  to  be  -O.OlSaqQ. 

The  value  of  2bq  / aq  that  was  computed  from  the 
mu 

experimental  data  was  1.07;  consequently,  the  static  pressure 

effect  accounts  for  approximately  one-fourth  of  the 

difference  between  2bq  /aq.  and  1.0.  This  implies  that 

either  one  or  both  of  the  terms  resulting  from  the  other  two 

effects  are  quite  important.  It  is  also  possible  that  the  sUtic 

pressure  assumption  is  in  error,  viz.,  P  =0,  however 

X  =  0  ’ 

it  is  not  expected  that  this  is  the  case. 

In  addition  to  the  velocity  profiles,  three  other 

distribution  comparisons  may  be  employed;  they  are 

9  >  and  C  .  The  distribution  of  /  u*  from 
tn  p  m 

the  present  study  has  been  compared  with  the  distributions 
of  three  two-dimensional  jets  in  the  preceding  section  (VI  -  B  -  1) . 
Figures  15  and  18  present  the  e(y/[  x-x^  ]  )  distribution  from  the 
present  study  as  compared  with  the  data  of  Heskestad.  The 
choice  of  x^  for  the  present  study  is  discussed  in  the  previous 
section.  Although  the  data  from  { 14)  is  Uken  at  an  x/ a  =  100, 
the  comparison  with  the  present  study  is  still  valid  since  the  0 
solution  is  similar  (for  the  two-dimensional  jet)  .  The  striking 
difference  between  the  distributions  for  the  two-dimensional 
and  the  bounded  jet  is  the  magnitude  of  0  for  small  y  values. 

This  may  indicate  that  the  vortex  stretching  action  is  not  the 
sole  difference  between  the  two  flows  since  the  0  differences 
cannot  be  explained  by  the  hypothesis. 
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The  Cp  distributions  as  a  function  of  IJ  are  shown 
in  Figures  19  to  2Z  along  with  the  data  from  (9).  The  Cp 
profiles  for  x/a  =  5  and  10  are  quite  different  from  those 
of  the  two-dimensional  jet.  The  distribution  at  x/a  =  5  resem¬ 
bles  the  shape  of  the  distributions  for  the  developing  flow 
region  of  (9).  This  may  indicate  that  the  bounded  jet  develops 
much  more  slowly  even  though  the  non-dimensional  velocity 
profile  at  x/a  =  5  appears  to  possess  a  similarity  solution, 

(see  Figure  11).  The  Cp  distribution  at  x/a  *  10  for  the 
bounded  jet  appears  to  follow  a  fairly  smooth  transition  in 
shape  from  x/a  =  5  to  x/a  =  15  but  it  does  not  correspond 
to  the  two-dimensional  distribution.  The  shape  of  the  curve 
and  the  magnitudes  of  the  values  for  the  bounded  jet  are  again 
quite  different  from  the  two-dimensional  jet  at  x/a  =  15, 

The  magnitude  of  the  difference  is  increased  but  the  curve 
shape  is  becoming  more  like  the  two-dimensional  case  for 
x/a  =  30.  If  the  absolute  values  of  the  static  pressures  were 
compared,  the  magnitude  differences  would  be  greater 
for  x/a  =  15  and  30  because,  with  reference  to  Figure  8, 
Um/uQis  greater  for  the  bounded  jet  at  the  same  x/a 
location.  Two  factors  are  of  importance  in  the  over-all 
comparison  between  the  bounded  and  the  two-dimensional  jet: 

(i)  the  lower  pressures  and  (ii)  the  wider  regions  about 
y  =  0  of  these  lower  pressures.  In  a  previous  section,  it 
was  stated  that  the  lower  vT^  /u*  values  of  the  bounded  jet 
for  the  neighborhood  of  y  =  0  might  be  a  result  of  more 
turbulence  energy  being  present  in  the  and  components. 

If  it  is  assumed  that  the  relationship  for  the  two-dimensional 

jet  (proposed  by  Townsend  (10)  and  verified  by  Miller  and 
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Comings  (9) )  viz. , 

P  +  p  =  P 

atm, 

provides  at  least  an  approximate  description  for  the  flow  at 
z/  a  s  0  for  the  bounded  jet,  then  the  large  differences  in  the 
static  pressure  magnitudes  would  result  from  a  relatively 
higher  .  This  would  support  the  earlier  contention  that 
more  turbulence  energy  is  contained  in  the  component. 

Conclusions  From  the  Study  of  the  Single  Bounded  Jet 
An  incompressible  turbulent  jet  issuing  from  a  nozzle 
of  rectangular  cross-section  with  an  aspect  ratio  of  6  and  bounded 
by  parallel  plates  normal  to  the  longer  sides  of  the  nozzle  is  highly 
three-dimensional  and  has  the  following  characteristics  (where  the 
jet  axis  is  in  the  x-direction  and  the  bounding  plates  are  planes  of 
constant  z) : 

i)  Complete  velocity  profiles  15  nozzle  widths  downstream 
from  the  nozzle  exit  show  that  in  planes  of  constant  y  removed  from 
y  =  0  the  maximum  velocity  occurs  near  the  bounding  plates  and  not 
midway  between  them.  For  example,  the  ratio  of  the  maximum 
velocity  to  the  centerplane  velocity  was  approximately  2.  5  at  a  y 
distance  of  2.5  nozzle  widths. 

ii)  The  non-dimensional  velocity  profiles  in  xy  planes  of 
the  bounded  jet  agree  with  those  of  the  two-dimensional  jet  of 
Reference  9. 

iii)  The  centerline  velocity  is  higher  and  the  momentum 
flux  width  is  lower  (downstream  of  the  potential  flow  core)  than 
for  a  two-dimensional  jet. 

iv)  In  the  high  velocity  region  of  the  jet  in  an  xy  plane, 
the  flow  angle  of  the  mean  velocity  (directed  away  from  the  y  =  0 
plane)  is  larger  for  the  bounded  jet  than  for  the  two-dimensional  jet. 
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v)  For  sufficiently  large  downstream  distances,  the 
flux  ratios  (mass,  momentum  and  energy  with  respect  to  the 
nozzle  exit)  in  xy  planes  near  the  bounding  plates  are  greater 
than  those  in  the  midplane  between  the  bounding  plates. 

vi)  For  all  x/  a  locations,  the  static  pressure  dis¬ 
tribution  (as  a  function  of  y)  shows  a  marked  change  from  that 
of  the  two-dimensional  jet.  Near  the  bounding  plates,  9P/8  y  =  0 
in  the  high  velocity  region  of  the  jet. 

vii)  The  distribution  of  across  i  bounded  jet  is 
unlike  that  across  a  two-dimensional  jet  (even  though  there  is 
lack  of  agreement  in  the  literature  for  the  two-dimensional 
distribution  of  d^  . 

The  vortex  stretching  hypothesis  which  was  formulated 
to  explain  the  above  item  i)  has  been  verified  by  or  found  to  be 
compatible  with  items  ii)  to  vii). 
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FIG.  6 


b/o  VERSUS  x/o 


FIG.  7  MASS,  MOMENTUM,  AND  ENERGY  FLUX 
RATES  FOR  THE  BOUNDED  AND  THE  TWO 
DIMENSIONAL  JET 


FIG.  8  uw/uo  FROM  (9),  (12)  , ( 16 ) ,  AND  THE  PRESENT  STUDY 
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NOTE:  I)  NUMERICAL  VALUES  ARE  NEGATIVE  INCHES  OF  WATER 
GAGE  PRESSURE  «  10*  ,  rg.  80  »  - (  P-  J  »I0*  in.  wot«r 

ID  CONTOURS  ARE  BASED  ON  420  DATA  POINTS 
ill)  1/2  1 1  04  I  10'*  in.  wattr 


*S  -2  -I  0  ►  2  3 


y /a 


FIG.  10  STATIC  PRESSURE  CONTOURS,  x/o  » 15 


u{TJ)/\i^  VERSUS  7J  WITH  THE  REICHARDT  SOLUTION 
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FIG.  13  u(y/a)/ue  VERSUS  y  FOR  x/o  =  15,  z/a  =  0,*2 
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FIG.  14  u(y/a)/Uc  VERSUS  y  FOR  x/o=30,  z/a=0,^2 
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FIG.  19  Cp  VERSUS  77  FOR  x/o»5 


FIG.  20  Cp  VERSUS  7]  FOR  x/a=IO,zA)=0 
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FIG.  22  Cp  VERSUS  rj  FOR 


X  /o  =  30 ,  z/a-0 


FIG.  24  Cp(y)  VERSUS  y  FOR  x/a  =  15 ,  z/o  =0,*2 
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FIG.  25  Cp(y)  VERSUS  y  FOR  x/a  =30 ,  z/a  =0,»2 
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FIG.  27  6*/u*  VERSUS  y/x  FOR  x/of  15  ,  zAj  =0,*2 


90 


JET  DEFLECTION  STUDY 


Analytical  Prediction  of  a 


and  P^. 


As  noted  in  the  literature  survey,  references  17  and 
18  present  analyses  which  predict  tana  from  known  input  and 
geometric  conditions  for  the  "defined  region  jet  interaction 
problem.  "  A  shortcoming  of  the  more  involved  analysis  in 
17  is  the  empirical  value  of  the  far  side  pocket  pressure  since 
no  reliable  relationship  (empirical  or  analytical)  between  this 
pressure  and  the  other  geometric  and  flow  variables  is  known. 

It  is  also  felt  that  the  computation  technique  for  the  near  side 
pocket  pressure  involves  rather  tenuous  assumptions  which  were 
unavoidable  because  of  the  method  of  attack  on  the  problem. 

In  the  analysis  presented  below,  the  problem  of  the 
apparent  piVot  point  is  skirted  by  using  an  empirical  formula 
(deduced  from  the  present  study).  By  so  doing,  the  formu¬ 
lation  for  the  near  side  pocket  pressure  of  (17)  was  avoided 
and  the  empirical  far  side  pocket  pressure  of  (17)  was  replaced 
by  the  empirical  apparent  pivot  point  formulation. 

Figure  29  indicates  the  geometry  considered  and 
defines  the  terms  to  be  used. 
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FIGURE  29  Definition  of  terms  used  in  analysis 

For  a  curved  flow  the  equation  of  motion  normal  to  a 
streamline  may  be  written  as 

dn  "  R  [  1  ] 

where:  R  is  the  radius  of  curvature  at  the  location  under 

investigation:  n  is  the  co-ordinate  normal  to  the  streamline. 

The  following  assumptions  are  made: 

i)  the  flow  is  steady  and  two-dimensional; 

ii)  there  is  no  jet  spreading  before  point  P; 

iii)  atmospheric  pressure  acts  on  the  open  surfaces 

of  the  control  volume; 

Pi  ■  P  * 

L _ atm 

a2 


iv)  ^ 
dn 


v| 

R* 
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dP, 


=  P. 


li 

Rl 


dn  a  1 

v)  some  circulatory  flow  exists  in  the  far  side  pocket 
such  that  the  pressure  on  the  walls  of  the  pocket  may  be  different 
from  atmospheric  pressure. 


Then  for  P  and  P  as  gage  pressures, 

1j  r 


p  Vi 


aa 


R  2 

=  p  V/  aj 


Ra 

Ri 


[2] 


Solving  for  P  gives 
R 

Pr  =  ^-31 

Ra  Ri  [3] 

where  q  is  the  momentum  flux  per  unit  depth. 

In  the  analysis  of  control  volume  1,  the  following 

assumptions  are  made; 

i)  constant  pressures  P  and  P  act  on  the  control 

Lj  r 

volume  surfaces  defined  by  the  setback  and  standoff  walls 
plus  a/  2. 

ii)  atmospheric  pressure  acts  on  the  other  half  of  the 
control  jet  nozzle  and  on  the  open  surface  of  the  control  volume. 
F rom  the  general  equation,  * 

>  >  rr>>> 

^  ^s  ®c.v.  J  )  V(pV.dA)  +  ^  j  I  j  p  V  dvol. 

the  component  equations  may  be  written  as; 
y  direction 

[Pj^  -  Pj^l  X  =  qj  sin  a  -  qa 


5  > 

*  F  represents  the  surface  forces  and  B  represents 
the  body  forces  acting  on  the  control  volume. 

c.  s.  stands  for  "control  surface" 
c.  V.  stands  for  "control  volume" 
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X  direction 

[  ^r]  «  -  qi 

and  solving  these  one  obtains 

t.n  „  .  |Pl  -  Pgl  X-  M.  .  [<] 

I  ’“l  *  '’rI  *"■ 

Static  pressure  measurements  made  in  the  present 
study  show  that  the  assumption  of  constant  pressure  on  the 
standoff  and  setback  walls  is  quite  reasonable  and  that  the 
assumption  of  atmospheric  pressure  on  the  open  surface  is 
satisfactory  (see  Figure  40) .  However,  the  assumption  of  con¬ 
stant  values  of  and  to  develop  equations  2  and  3  is 
probably  in  serious  error  near  the  point  P.  To  compensate 
for  this,  control  volume  2  can  be  analyzed  as  if  the  net 
pressure  forces  are  zero.  This  assumption  then  allows 
for  the  momentum  flux  components  to  be  added  (vectorially) 
yielding  a  second  equation  for  tan  a. 

The  y  component  of  the  momentum  equation  yields 
q3  sin  a  -  qi  sin  y  +  qz  cos 
For  the  x  component  we  obtain 

qs  cos  a  =  qi  cos  y+  qz  sin  . 

Combining  these  gives 

tan  a  =  qi  sin  y  i-  qz  cos  p 

qi  cos  y  +  qz  sin  ^  [5] 

Two  more  equations  may  be  obtained  from  purely 
geometric  conditions 

Rl  +  Y'  =  Rj  cos  y+  Rz  sin  ^ 
or 


1  +  V/R  1  r  COS')  +  Rz/R  1  sin  p 


[6] 
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and 

R2  +  X*  r  Rising  +  R2COS 
or 

1  +  X'/R  =  (R1/R2)  sin-yf  cos  [7] 

At  this  point,  6  independent  equations  for  7  variables 
are  available.  A  relationship  for  P  could  be  introduced 
which  would  be  patterned  after  the  analysis  of  ( 17)  .  However, 
since  this  is  felt  to  be  a  weakness  in  the  analysis  of  ( 17) ,  this 
will  not  be  done.  One  new  variable,  the  apparent  pivot  point, 
will  be  introduced  which  will  provide  an  additional  geometric 
equation  and  which  will  necessitate  the  introduction  of  an 
empirical  relationship.  This,  of  course,  sacrifices  the  purely 
analytical  character  of  the  analysis.  The  sketch  below  illustrates 
the  geometry  for  the  equation  involving  the  apparent  pivot  point. 


FIGURE  30  Definition  of  the  apparent  pivot  point 

'  sin  a  =  R  1  ( 1  -  cos  y  ) 

cos  a  +  X  =  R 1  sinv 
app 

tan  a  =  R  1  ( 1  -  cos y  ) 

R  I  sin  Y  -  X 

app 


[8] 
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The  empirical  relation 


X 

_i£E 

X' 


=  0.25  k 


Y 

Y' 


+  0.875 


X 

X' 


[9] 


where:  k  e  ( Y/  X) 

has  been  deduced  from  experimental  data  of  the  present  study. 

Thus,  there  are  8  equations  involving  the  8  variables 
<vU..  R.  R.  Pj_. 

An  outline  of  the  reduction  of  these  equations  is 
presented  in  Appendix  C;  the  8  equations  are  reduced  to  2 
equations  as  R/  =  Rl(y)  where  R'l  =  Rj/X*.  A  computer 
solution  was  employed  to  determine  the  correct  Ri  as  a 
function  ofy  for  a  given  k  and  mfr,  and  the  procedure  was 
then  repeated  for  mfr  values  of  0.05,  0.  1,  0.  15,  0.2,  0.25 
as  well  as  k  values  of  0.5,  1.0  and  1.5.  The  results  for 
tan  a,  C  ,  C  ,  R'l,  and  R  2  from  these  calculations  are 

Pr 

shown  in  Table  6.  C  and  C  are  defined  in  equations  2a 

Pr 

and  3a  below;  they  provide  a  dimensionless  representation  of 

P  and  P_ . 

L  R 


C 


X'  P. 


q  I 


mfr 

Ri 


[  2a  ] 


C  =  X'P_  =  mfr  -  1 

Pr  ^  -rT 


[  3a  ] 
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Method  of  Testing  and  Data  Evaluation 
For  each  geometry,  the  following  experimental  method 
was  employed  to  determine  the  value  of  a  as  a  function  of 
the  momentum  flux  ratio: 

i)  The  mfr  value  was  fixed  by  setting  the  dynamic 
head  at  the  control  jet  nozzle  exit  at  a  given  fraction  of  the 
power  jet  dynamic  head.  Since  both  nozzle  exits  were  of 
the  same  exit  area,  the  dynamic  head  ratio  was  equivalent 
to  the  mfr.  The  control  jet  setting  was  made  with  no  flow 
through  the  power  jet  nozzle. 

ii)  The  x  -  y  -  0  device  was  then  used  to  determine 
the  position  and  flow  direction  of  the  maximum  velocity  of 
the  resultant  jet  (in  the  z  =  0  plane) .  The  m.'  ximum  velocity 
location  was  determined  from  a  null  reading  of  a  dual-impact 
probe  (  spacing).  These  data  were  taken  at  x/a  values 
of  5,  10,  15,  20,  30. 

iii)  A  new  mfr  value  was  then  set  and  the  steps  (i) 
to  (iii)  were  repeated  until  all  the  mfr  values,  viz.  0.05, 

0.  10,  0.  15,  and  0.  20  were  investigated. 

iv)  The  control  jet  flow  was  then  turned  off  and  the  un- 
deflectsd  values  of  y  and  0  were  recorded. 

Plots  of  (  a)  vs  (x/  a)  were  then  prepared  for 

each  mfr  setting.  From  these  the  value  of  tan  a  and  the 
apparent  pivot  point  could  be  determined;  a  sample  plot  is 
shown  in  Figure  3  1. 

The  above  technique  was  employed  in  the  z/  a  s  0 
plane.  However,  this  does  not  limit  the  generality  of  the  results 
since  it  was  determined  that  the  resultant  jet  was  deflected 
evenly  as  a  function  of  z.  This  observation  is  in  agreement 
with  the  results  of  ( 17)  . 
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Results 

Values  of  Cp.  *j^pp  from  the  present  study 

are  given  in  Table  4.  *  Four  definitions  of  Cp  are  presented  in 
order  to  ascertain  whether  any  of  these  would  exhibit  the  sanne 
functional  character  for  the  different  geometries.  Experimental 
values  of  the  tan  have  been  determined  by  the  methods  described 
in  the  previous  section.  The  agreement  between  these  two  inde> 
pendent  methods  of  measuring  tan  is  excellent.  Table  5  pre¬ 
sents  a  comparison  of  the  experimental  values  of  Cp  ,  Cp  and 

Lj  r 

tan  with  those  predicted  by  the  analysis  of  reference  ( 17)  and 
the  analysis  of  section  VII-A.  A  comparison  is  also  presented 
wherein  the  value  of  tan  is  computed  by  the  control  volume 

equation  of  (18)  using  the  measured  values  of  Cp  and  Cp  . 

Li  R 

Figures  32  to  37  offer  a  graphical  comparison  of  the  tan 
and  the  Cpj^  values. 

Table  6  is  included  to  present  the  complete  results 
from  the  computer  solution  of  the  analytical  prediction  equations. 
As  shown  in  this  table,  there  Mere  two  possible  solutions  for 
most  values  of  k  and  the  mfr.  The  "correct"  one  was  chosen 
as  that  which  gave  the  lower  value  of  tan  The  other  solution 
may  represent  a  physically  unstable  case. 

Manion  and  Goto  determined  tan  as  a  function  of  the 
mfr  for  three  geometries;  X  e  Y  =  la,  X  =  Y  =  3a  and  X  =  la, 

Y  =  3a:  three  power  jet  plenum  pressures  of  5,  10  and  20 

psig  ..ere  considered  fo.-  each  geometry.  Since  the 
geometries  differ  from  those  of  the  present  study,  a  direct 
comparison  is  not  possible;  however,  a  valid  comparison  may 

*  The  tables  that  present  the  results  of  the  jet 
deflection  study  can  be  found  in  the  "Discussion  of  Results" 
section. 

•  *  The  figures  that  present  the  results  of  the 
jet  deflection  study  can  be  found  at  the  end  of  the  section, 
see  page  112  ff. 
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be  to  compare  the  results  for  X  =  Y  =  2a  (present  study)  with 
those  for  X  =  Y  =  la  and  3a  of  ( 17)  .  Figures  and  39  show 
this  comparison. 

A  static  pressure  survey  was  made  with  the  plate 
static  probe  (from  Reference  9)  in  the  interaction  region  for 
the  condition  X  =  Y  =  2a.  and  a  mfr  =  0.  1.  The  results  of 
this  survey,  which  are  shown  as  h^  (x,  y)  for  z/asO  in  Figure 
40,  should  be  considered  as  a  presentation  of  approximately 
true  values  of  the  static  pressure  with  the  error  being  depen¬ 
dent  on  the  proximity  of  the  probe  to  the  jet.  This  probe  is 
highly  sensitive  to  pitch  and  there  was  no  reliable  way  to 
align  the  probe;  however,  the  probe  and  the  wall  taps  did 
give  the  same  reading  for  a  probe  position  ■-  -  inch  from 
the  wall.  This  indicates  that  in  the  regions  which  are  not  in 
the  jet  flows  the  probe  readings  should  be  rather  accurate. 

The  probe  was  also  used  to  traverse  the  z  direction;  no 
significant  changes  in  the  h^  reading  were  noted  with 
respect  to  z. 


Discussion  of  Results 
Effect  of  Probe  Spacing 

With  regard  to  the  experimental  results  of  this  section 
it  is  important  to  note  an  arbitrary  definition  which  necessarily 
was  made  at  the  beginning  of  the  tests.  Since  a  dual  pitot 
probe  with  a  fixed  spacing  of  0.  25  inches  was  used  for  all  x 
staf  ons,  it  must  be  acknowledged  that  the  quantity  measured 
i  simply  the  y  location  where  3u<  h  a  probe  gives  a  null 
reading.  The  intent,  of  course,  is  to  locate  the  y  position 
of  maximum  velocity.  If  the  total  pressure  profile  and  the 
static  pressure  profile  were  symmetric  about  the  maximum 
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velocity,  the  spacing  of  the  dual  probe  would  be  immaterial. 
However,  neither  of  these  two  conditions  is  realized;  the 
asymmetries  are  most  pronounced  for  high  mfr  and  low 
x/  a  values.  It  can  be  reasonably  argued  that  no  single  quan¬ 
tity  provides  a  complete  description  of  the  jet  deflection. 

For  example,  the  maximum  velocity  location  and  the  loca¬ 
tions  which  divide  each  of  M,  2bq^,  and  E  into  two  equal 
parts  would  yield  four  unequal  values  of  y  for  the  center 
line  of  the  resultant  jet.  Consequently,  the  use  of  a  dual 
probe  is  an  arbitrary  selection  of  the  location  of  the  resultant 
jet  and  this  particular  selection  is  a  matter  of  convenience. 

Because  of  the  small  size  of  the  experimental  flow 
system  and  physical  limitations  on  the  smallness  of  their  probe, 
Manion  and  Goto  employed  a  dual  probe  with  a  spacing  of  2.  88a. 

It  would  seem  that  the  probe  response,  for  the  2.  88a  spacing, 
would  be  insensitive  to  positions  at  low  x/ a  values  since  the 
majority  of  the  jet  is  contained  within  the  probes.  The  dual 
probe  of  the  present  study  (spacing  of  0.  25a)  was  rather  in¬ 
sensitive  to  location  changes  at  x/ a  s  30  because  of  the  flat 
velocity  profile.  It  would  also  appear  that  a  probe  with  a 
spacing  of  2,  88a  would  be  more  susceptible  to  a  skewed  velo¬ 
city  profile  if  such  a  profile  resulted  from  the  interaction  process. 

Interpretation  of  Results  from  the  Present  Study 

The  most  significant  result  to  be  gained  from  the 
present  study  (see  Table  4)  is  the  pronounced  effect  of  the 
geometry.  This  can  be  seen  by  noting  the  mfr  necessary  to 
obtain  the  same  value  of  tan  a  for  the  three  geometries.  A 
mfr  of  0.  10,  0.  15,  and  0.  20  will  yield  (approximately)  the  same 
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TABLE  4 

Results  of  the  Jet  Deflection  Study 


mfr 

X 

Cp 

s 

Cp 

Tan  a 

app 

i^-cj  Lrpj 

R-cj 

R-pj 

F  rom 

F  rom 

A  e 

X/a  = 

2 

Y/a 

=  1 

0.  05 

1.9 

+  0.  54 

+  0.  024 

-0.  55 

-0.  028 

0.  120 

0.  121 

0.  10 

2.  0 

+  0.  54 

+  0.  065 

-0.  33 

-0.  034 

0.  235 

0.  236 

0.  15 

1.9 

+  0.  60 

+  0.  09 

-0.26 

-0. 038 

0.295 

0.  302 

0.  20 

2.  0 

+  0.  59 

+  0.  120 

-0.  24 

-0.  047 

0.  362 

0.  375 

X/a  = 

2 

Y/a 

=  2 

0.  05 

2.4 

0.  24 

0.  013 

-0.29 

-0.  014 

0.  081 

0.  082 

0.  10 

2,2 

0,  37 

0.  037 

-0.  14 

-0.  014 

0,  149 

0,  15 

2.  4 

0.41 

0.  061 

-0. 095 

-0.  014 

0.23 

0.  233 

0.  20 

2.  2 

0,  41 

0.  082 

-0. 072 

-0. 014 

0.29 

0.  294 

X/a  = 

2 

Y/a 

=  3 

0.  05 

2.  2 

0.  170 

0.  0085 

-0. 142 

-0.  007 

0. 0625 

0.  061 

0,  10 

2.4 

0.  254 

0.  0254 

-0. 066 

-0. 0066 

0.  119 

0.  118 

0.  15 

2.6 

0.  287 

0.  0430 

-0. 053 

-0.  0080 

0.  192 

0.  186 

0.  20 

2.8 

0.  294 

0.  0586 

-0. 047 

-0.  0093 

0.247 

0.  233 

Cp  -  P.  .  1/2  pvj  where  Pj  =  Average  Wall  Pressure 

(i*k) 
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tan  a  *  for  the  geometries  X/  a  =  2,  Y/ a  =  1,  2,  and  3 

respectively.  As  previously  noted,  the  actual  jet  deflection 

depends  not  only  on  tano  but  also  on  the  apparent  pivot 

point.  Comparing  the  approximate  x  values  for  the 

app 

above  three  cases,  one  sees  that  for  smaller  Y  values, 
not  only  is  tana  larger,  but  the  apparent  pivot  point  is 
closer  to  the  power  jet  nozzle. 

In  the  analysis  of  the  data  in  Table  4  it  should  be 
borne  in  mind  that  if  the  momentum  of  the  control  jet  were 
the  only  deflecting  agent,  then  the  value  of  the  mfr  and  the 
tan  a  would  be  equal.  The  difference  between  these  two  can 
then  be  used  as  a  measure  of  the  static  pressure  effect 
induced  by  the  defined  region  geometry.  For  example, 
comparing  tan  a  at  mfr  values  of  0.  05  a  d  0.  2  for  the  three 
geometries,  one  obUins  for  the  ratio  (tan  a  -  mfr)  /  mfr 
the  following  values; 


Table  7 

Values  to  Indicate  the  Effect  of  the  Static  Pressure 

tan  a  -  mfr 
mfr 


mfr 


X/a _ Y/a _ 0^.  05  0.  20 


1 - 

2 

1 

1.  40 

>  o 

1  00 

>  o 

2 

2 

0.60 

0.  45 

2 

3 

0.  25 

0.  24 

+  For  this  and  subsequent  comparisons  the  value  of 
tan  a  is  taken  from  the  A  y  data  except  as  noted. 
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Another  method  by  which  to  gain  an  appreciation  of 

the  magnitude  of  the  pressure  effects  is  to  compare 

^  “  C  values  for  the  different  geometries. 

-  pj  Pr  -  pj 

The  difference  is  important  because  it  is  the  net  pressure  force 
which  contributes  to  the  actual  deflection.  From  the  deflection 
angle  results  it  is  apparent  that  the  geometry  X/a  =  2  and 
Y/a  =  1  is  the  most  effective  in  obtaining  a  maximum  de¬ 
flection  for  a  given  mfr.  These  results  indicate  that  the  gain 
of  the  proportional  amplifier  increases  with  decreasing  Y/ X 
values;  however,  the  range  over  which  this  generalization  is 
valid  has  only  been  investigated  for  Y/X  values  greater  than 
or  equal  to  —  and  was  made  for  data  which  employed  a  con¬ 
stant  X  value.  Further  investigations  are  necessary  to 
reliably  extend  the  range  of  the  generalization. 

With  regard  to  the  values  of  C  *,  it  appears 

that  no  simple  method  of  non-dimensionilizytion  offers  a 
reliable  prediction  of  the  pocket  pressures.  If  a  non- 
dimensionalization  technique  is  to  yield  a  predictable  C 

P 

value,  it  is  obvious  that  the  Y  length  as  well  as  some  com¬ 
bination  of  the  power  and  control  jet  dynamic  head  or  mo¬ 
mentum  fluxes  must  be  incorporated  in  the  definition.  This 
was  not  pursued  because  of  the  rather  limited  experimental 
data  with  regard  to  X,  X',  Y,  and  Y' . 

From  other  literature  sources  (e.  g.  (1)  to  (3)) 
the  mfr  range  covered  in  the  present  study  appears  to  encom¬ 
pass  the  range  of  current  interest  for  amplifier  applications. 

Comparison  with  Published  Data 

Manion  and  Goto  ( 17)  present  the  only  data  which 
can  be  utilized  for  a  direct  comparison  with  the  results  of  the 

*  C 

P(i-k) 


is  def  ined  in  Table  4. 
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present  study.  As  noted  in  the  literature  survey,  their  work 
was  for  compressible  flow  and  an  aspect  ratio  of  8.  With 
reference  to  Figures  38  and  39,  it  can  be  seen  that  the  de¬ 
flection  angles  of  the  present  study  are  of  the  same  order 
as  the  maximum  angles  found  in  (17).  Furthermore,  the 
data  of  (17)  are  for  three  power  jet  plenum  pressures  and 
the  pattern  "a  greater  plenum  pressure  yields  a  greater  angle 
a  for  the  same  mfr"  is  clearly  established.  Three  state¬ 
ments  can  be  made  with  regard  to  this  observation;  it  cannot 
be  ascertained  which  statement  or  which  combination  of  state¬ 
ments  provides  the  most  accurate  representation  of  the 
physical  case  without  further  experimental,  or  more  refined 
analytical,  investigations.  However,  the  statements  do 
have  significance  for  the  designer.  They  are: 

i)  Assuming  that  the  trend,  in  (17),  of  increasing 
a  for  an  increasing  at  a  given  mfr  may  be  extrapo¬ 

lated  to  the  incompressible  regime,  then  the  large  angles 
(  Of  )  of  the  present  study  are  attributable  to  the  strong 
effect  of  the  aspect  ratio.  That  is,  the  gain  of  the  propor¬ 
tional  amplifier  is  higher  for  smaller  aspect  ratios. 

ii)  Assuming  that  the  aspect  ratio  effect  between 
6  and  8  is  of  negligible  importance,  the  incompressible 
regime  is  seen  to  yield  better  amplifier  performance  than 
the  compressible  regime. 

iii)  The  comparison  with  (17)  is  predicated  on  the 
assumption  that  the  performance  of  X  =  Y  =  2a  should  lie 
between  X  r  Y  =  la  and  X  =  Y  =  3a,  all  the  other  factors 
being  equal.  The  physical  phenomena  involved  may  not  be  amen¬ 
able  to  such  a  simple  analysis.  This  would  mean  that  the 


designer  must  be  rather  careful  in  extrapolating  and  interpo¬ 
lating  his  experience  with  previous  configurations. 

An  insight  into  the  jet  deflection  phenomenon  can  be 
gained  if  the  results  of  the  present  study  are  compared  with 
those  of  Olson  (22).  The  apparent  pivot  point  for  the 
majority  of  the  tests  in  the  present  study  was  upstream  of  the 
geometric  intersection  whereas  in  (22)  it  lay  downstream  of 
this  intersection.  Olson  attributed  the  observed  displace¬ 
ment  to  a  time  delay  in  the  momentum  exchange  process. 

This  explanation  is  not  dependent  on  the  geometry  considered 
and  therefore  the  same  time  delay  should  occur  in  the  inter¬ 
action  process  of  the  present  study.  The  effect  is  not 
observed  in  the  majority  of  the  runs;  however,  for  Y  s  3a 
the  apparent  pivot  point  does  move  downstream  with  in¬ 
creasing  mfr  values.  These  two  observations  suggest 
that  the  apparent  pivot  point  is  strongly  related  to  the 
relative  magnitude  of  the  pocket  static  pressure.  As  in¬ 
dicated  in  Table  7  the  static  pressure  induced  deflection 
(as  a  percentage  of  the  deflection  due  to  the  momentum 
flux)  is  25  per  cent  for  Y  =  3a  and  80  to  140  per  cent  for 

Y  =  la.  Consequently  the  larger  x  values  for  Y  =  3a 

app 

are  in  agreement  with  Olson's  results  since  the  momentum 
interchange  is  more  important  for  this  geometry.  However, 
Olson  also  notes  that  the  displacement  is  reduced  for  in¬ 
creasing  mfr  values;  a  trend  which  is  the  inverse  of  the 
trend  observed  in  the  present  study.  A  possible  conclusion 
from  this  is  that  the  presence  of  the  walls  and  the  greater- 
than-atmospheric  pressure  in  the  near  side  pocket  of  the 
present  study  induces  a  change  in  the  momentum  exchange 


process  of  the  intersecting  jets.  Another  important  factor 

in  the  comparison  is  that  a  power  jet  exit  Mach  number  of 

0.  66  was  used  in  (22) .  However,  there  is  no  apparent 

reason  why  this  would  induce  the  observed  trend  of  the 

X  with  the  mfr. 
app 


Comparison  with  the  Analytical  Prediction 

A  principal  requirement  of  a  successful  analytical 
model  is  that  it  provide  an  accurate  representation  of  the 
physical  phenomenon  described.  Figures  32  to  37  and  Table 
5  offer  this  comparison  for  the  model  proposed  here  and  also 
present  the  results  of  the  Manion  and  Goto  analysis  (  17) .  The 
latter  were  obtained  from  a  computer  solution  of  the  equations 
in  ( 17) ;  the  nozzle  geometry,  the  mfr,  and  the  right  side 
pocket  pressure  were  the  required  data  inputs  for  the  program* 
The  present  analysis  for  tan  a  gives  a  rather  good 
approximation  for  k  =  0.5  ;  the  agreement  is  within  10  per  cent 
up  to  a  mfr  of  0,  2.  However,  for  a  k  of  1.0  the  agreement  is 
only  within  25  per  cent  up  to  a  mfr  of  0.  20.  It  is  fortuitous 
that  the  agreement  is  best  for  a  k  of  0.  3  since  this  geometry 
offers  the  best  performance  of  the  three  tested.  The  agree¬ 
ments  shown  for  k  =  1.0  and  1.  5  are,  of  course,  better  than 
no  prediction  whatever,  but  would  probably  not  be  sufficient 
for  the  design  of  an  operational  amplifier. 

*  The  writer  is  indebted  to  che  Harry  Diamond 
Laboratories  for  the  computer  solutions  of  the  Manion  and 
Goto  analysis. 


were  computed  by  the  control  volume  equation  of  (18)  using  the 
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Tan  a  and  the  location  of  the  apparent  pivot  point 
will  Define  the  location  of  the  resultant  jet;  consequently, 
agreement  of  the  tan  a  values  is  the  most  important  cri¬ 
terion  for  the  analytic  lodel.  However,  an  appreciation 
for  the  accuracy  of  th  lysis  may  also  be  gained  from 
the  predicted  pocket  pr  ,«res.  This  comparison  is  offered 
in  Figures  35  to  37  and  Table  5,  The  Cpj^  values  were 
not  plotted  because  of  the  wide  discrepancy  between  the  • 
observed  and  predicted  values,  a  difference  that  is  obvious 
from  the  table.  Jt  is  interesting  to  note  that  the  Cj, 
value  that  shows  the  largest  discrepancy  is  for  k  =  0,  5, 
the  geometry  for  which  tan  a  showed  the  best  agreement. 

The  agreement  of  the  experimental  values  for  k  b  0.5  and 
1.0  with  the  analysis  for  k  *  1.0  and  1.5  (respectively)  is 
interesting;  but,  it  is  not  believed  to  be  indicative  of  any 
significant  phenomenon. 

A  striking  feature  that  is  apparent  from  the  figures 
is  the  good  agreement  between  tan  o  values  predicted  by  the 
Manion  and  Goto  analysis  and  measured  values  for  the  geo¬ 
metry  of  k  r  1.  5  and  the  poorer  agreement  for  lov/er  k  values. 
The  analysis  of  ( 17)  gave  the  best  P  prediction  for  k  s  0.  5 

Xx 

and  poorer  agreement  for  increasing  k  values.  This  pattern 
is  the  inverse  of  that  between  the  present  analysis  and  obser¬ 
vations. 

The  analysis,  as  presented  in  the  earlier  section, 
has  several  arbitrary  steps  and  assumptions  which  were 
necessary  to  obtain  a  complete  formulation  of  the  problem. 
The  test  of  the  **correctnes8 ■'  of  these  is,  of  course,  the 
agreement  w  th  experimental  results.  After  one  finds  the 
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d'  iree  of  agreement,  he  can  then  either  chana 
easily  done.  Reco»«^  j 

ar.  .,v.„  ■»«  a..™p. 

8  irst,  the  problem  of  modifyino  k  • 
•-alyi.  i.  U,e„  co„.ida„^, 

-u.d  b.  „p..c.d  W  a  cap.  •"' 

at-  ^  ‘^onatant  pressure  over  X  ai 

the  near  side  nock-* 

*n....d  r  »'  “—'O  >■>  ao,  „.P  V  . 

aP  .0  pain.  p',.ad, 
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valaa.  « 

•implinad;  hawav.a 

•■ant  P  and  P  “•"“Pt--  »■  ■<>.- 

„a. 

a.aociatadwithplin.p”'^,^''''  °'®'’  ^ 
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(2c) 

(30 

To  modify  the  basic  analysis  ‘i, 

coiiiH  K  .  y®ia>  -he  method  of  1171 

could  be  employed  to  predict  th-  ' 

and  this  would  eliminate  the  needTr^thT' 

----^Vaes  have  inverse  ruccr::^ 

combination  of  the  two  analytical  model  •  ^ 

more  satisfactory  set  of  h-  •  *  Provide  a 

if  more  were  known  aboufthe  in^” 

-  questionable  assumptions^^ ^  ~ 

(17).  of  no  entrainment  and 
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no  penetration  of  the  control  jet  could  be  eliminated,  leaving 
the  basic  part  of  the  analysis  which  is  quite  plausible. 

The  analysis  produced  two  possible  solutions  for 
most  combinations  of  k  and  the  mfr.  If  there  is,  in  the  physical 
sense,  another  solution  that  is  possible,  it  could  lead  to 
important  implications  for  the  operational  amplifier.  The 
tan  a  values  for  the  second  case  was  much  higher  than  the 
solution  chosen  as  the  correct  one.  It  is  assumed  that  the 
second  solution  represents  an  unstable  case.  The  signifi¬ 
cance  to  the  operational  amplifier  would  be  that  it  may  be 
possible  to  bias  the  amplifier  with  a  low  pressure  in  the 
far  side  pocket  and  have  it  operate  cn  the  second  or  higher 
gain  solution  curve. 
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TABLE  6 


Complete  Solution  for  the  Jet  Deflection  Analysis 


Possible  Solutions  for  Equations  2a  -  9  for  6  =  0.  2k  +  0.  7 


k 

mfr 

y 

^2* 

tan  a 

^Pl 

Cp 

Pr 

0.  5 

0.  05 

6.8 

12.5 

0.  602 

0.  130 

0,  0829 

0. 0026 

0.  5 

0.  05 

34.4 

1.97 

3.22 

1.  10 

0.  0156 

-0. 492 

0.  5 

0.  10 

11.  7 

7.  35 

0.682 

0.  233 

0.  146 

0.  0102 

0.  5 

0.  15 

15.  6 

5.55 

0.  752 

0.  324 

0.  200 

0.  0193 

0.  5 

0.  20 

18.  9 

4.61 

0.815 

0.405 

0.  246 

0.  028 

0.  5 

0.  25 

21.8 

4.01 

0.875 

0.  481 

0.  286 

0.  037 

1.0 

0.05 

4.  68 

22.  3 

1.  13 

0.  096 

0.  0442 

-0. 00138 

1.0 

0.  05 

15.9 

4.64 

2.72 

0.  441 

0.  0184 

-0. 198 

1.0 

0.  10 

8.  36 

12. 1 

1.22 

0.  183 

0.  0819 

-0.  0097 

1.  0 

0.  10 

18.8 

4.  11 

2.47 

0.  579 

0.  0409 

-0.203 

1.0 

0.  15 

11.6 

8.68 

1.31 

0.266 

0.  115 

-0.0088 

1.0 

0.  15 

21.2 

3.  76 

2.  37 

0.  700 

0.  0635 

-0.202 

1.0 

0.  20 

14.  6 

6.81 

1.40 

0.  35 

0.  143 

-0. 0036 

1.0 

0.  20 

23.  3 

3.55 

2.  32 

0.  807 

0. 0844 

-0.  1987 

1.0 

0.25 

16.  9 

5.86 

1.42 

0.  421 

0.  171 

-0.  0003 

1.0 

0.  25 

25.  3 

3.  31 

2.31 

0.  874 

0.  108 

-0. 194 

1. 50 

0.  05 

3.  79 

29.8 

1,67 

0.  0818 

0. 0303 

-0. 0013 

1. 50 

0.  05 

12.  6 

6.23 

4.04 

0.  382 

0. 0124 

-0. 150 

1. 50 

0.  10 

6.95 

17.  1 

1.77 

0.  159 

0. 0565 

-0. 020 

1.  50 

0.  10 

15.  0 

5.61 

3.60 

0.  535 

0. 0287 

-0. 155 

1. 50 

0.  15 

9.82 

11.9 

1.88 

0.  238 

0.  0/97 

-0. 0056 

1. 50 

0.  15 

17.  1 

5.20 

3.  10 

0.  603 

0.  045 

-0.  149 

1. 50 

0.20 

12.  4 

9.  23 

1.99 

0.  317 

0.  10 

-0. 0083 

1. 50 

0.  20 

18.8 

4.94 

2.96 

0.  676 

0. 0682 

-0. 135 

1. 50 

0.  25 

14.8 

7.66 

2.  04 

0.  375 

0.  124 

-0. 0125 

1. 50 

0.25 

20.  2 

4.75 

2,83 

0.  709 

0. 0836 

-0. 122 
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FIG.  32  tan  a  VERSUS  mfr  ,  EXPERIMENTAL 
AND  ANALYTICAL  VALUES  ,  k  «0.5 
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FIG.  33  ton  Q  VERSUS  mfr  ,  EXPERIMENTAL 
AND  ANALYTICAL  VALUES  ,  K  =  1.0 


114 


FIG.  34  tan  a  VERSUS  mfr  ,  EXPERIMENTAL 
AND  ANALYTICAL  VALUES  ,  k  =  1.5 
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FIG.  35  Cp^  VERSUS  mfr  ,  EXPERIMENTAL 
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FIG.  39  COMPARISON  OF  tan  Q  FROM 
AND  THE  PRESENT  STUDY 


POWER  JET  FLOW 

FIG.  40  STATIC  PRESSURE  DISTRIBUTION  FOR  mf  r  =0.1  ,  X=Y=  2a 
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DENOTES  MEASURING 


DETAILED  MEAN  FLOW  SURVEY  -  RESULTANT  JET 


Methods  of  Measurement  and  Data  Evaluation 


Flow  Angle  Measurennents 

Two  methods  of  measurement  were  necessary  to 
>(■ 

obtain  complete  u(y)  data.  The  x  -  y  -  G  traverse  could 
only  be  used  for  direct  measurements  at  z/  a  =•  0  and  +2. 
Consequently  another  method  had  to  be  devised  to  obtain 
the  flow  direction  for  z/a  =  +  1,  -2,  and  +  2j.  The  tri¬ 
dent  probe,  mounted  in  the  6  =  0  position,  was  used  for 
these  angular  measurements. 

After  mounting  the  trident  probe  on  the  x  -  y  -  z 
traverse  device,  a  traverse  was  made  to  obtain  the  dif¬ 
ferential  reading  from  the  two  46  degree  impact  tubes.  The 
total  pressure  was  also  recorded  for  this  traverse.  This 
differential  reading  was  denoted  as  Ah  .  In  a  separate 

”45 

series  of  calibration  runs  (see  Appendix  B)  the  relation¬ 
ship 

e  =  ©{Ah  ) 

45/  q 

was  determined.  Consequently,  if  the  values  of  q  and 

A  h  are  known  as  a  function  of  y,  the  6(y/a) 

45 

distribution  can  be  considered  to  be  known.  A  discussion 

of  the  calibration  technique  and  the  presumed  accuracy  of  this 

method  as  well  as  the  comparison  of  the  u(y  /  a)  /  u  and 
*  •- 
0  (y)  distributions  from  the  two  methods  are  presented  in 

Appendix  B. 


122 


Data  Processing 

Since  two  measurement  methods  were  used,  the  data 
processing  was  similarly  of  two  distinct  types.  These  will  be 
discussed  separately:  they  are  referred  to  as  method  1  for 
the  data  collected  by  the  x  -  y  -0  device  and  method  2  for 
the  data  collected  by  the  x  -  y  -  z  device. 

Method  1  was  the  more  direct  in  that  the  flow  angle 
and  the  total  pressure  could  be  measured  directly.  The  static 
pressure  was  measured  in  a  separate  traverse.  In  order  to 
ensure  a  matching  of  these  separate  traverses,  h^  was 
measured  in  the  same  traverse  as  h  and  the  h_(y*)  from 

d  X 

the  X  -  y  -  0  device  was  then  compared  with  the  h,j.(y*)  from 
the  X  -  y  -  z  device.  By  this  procedure,  the  matching  of 
the  corresponding  y  values  depended  upon  matching  15  to 
20  measurements  from  each  traverse  and  is  considered  to 
be  quite  reliable  in  a  statistical  sense.  The  liata  were  then 
considered  to  be  ready  for  the  computer.  A  discussion  of  the 
computations  performed  will  be  deferred  until  the  description  of  method 
2  is  given. 

Two  traverses  were  also  required  for  method  2.  The 
first  traverse  was  made  with  the  tapered  total  and  the  forward 
facing  static  probes  mounted  on  the  x  -  y  -  z  traverse  device; 
h^  and  were  measured  on  the  second  using  the  tri¬ 

dent  probe  m  ^he  0=0  position.  A  matching  of  the  h^(y*) 
profiles  similar  to  ihat  of  method  1  ensured  a  matching  of  the 
y*  values.  The  data  for  method  2  was  then  compiled  as  h^, 

h  ,  and  Ah  as  a  function  of  y*. 

^  ^  45 

The  quantities  computed  from  all  traverses  were; 


(a)  as  a  tunction  of  y* 

i)  Q  V,  u,  V,  P/-i  pu  *  ,  q 

4  C 

ii)  u/  u  ,  (u/  u  )  (v/  u  ) 

c  c  c 

iii)  u  /u  ,  (u  /u 
^ 

(b)  for  each  traverse 


max 

i)  ^  (u/u^)^  dy 


and 


max 


max 


"y  , 


max 
^x 


im 


ii)  \  ^4 


and 


max 


max 


-L 


for  N  =  1  and  2 


is  the  mean  velocity  component  parallel  to  the 
axis  of  the  resultant  jet;  u  is  the  maximum  value  of  u  . 

The  integrations  in  part  (b)  above  yield  mass,  mo¬ 
mentum  and  energy  flux  thickness  with  respect  to  the  fixed 
geometric  co-ordinates  x,  y,  z  ai.d  the  streamline  co-ordin¬ 
ates  [  ,  4,  z. 

The  computer  program  written  for  these  computations 
is  presented  in  Appendix  D.  Because  of  the  differences  in  the 


experimental  techniques,  the  program  has  two  input  paths; 

.however,  once  the  value  of  6  is  known,  the  computational 

procedure  for  the  values  of  V,  u,  v,  (u/  u  ) ,  (v/  u  ),  (u/  u  )  *  , 
I  c  c  c 

and  P/y  P  were  quite  straightforward. 

Since  the  measurements  were  made  with  respect  to 
the  geometric  variables  x,  y,  z,  the  velocity  components 
u  and  V  are  given  with  respect  to  these  co-ordinates. 
However,  if  one  wishes  to  compare  the  resultant  jet  profile 
with  the  undeflected  jet  or  with  a  two-dimensional  jet  this 
comparison  could  also  be  made  with  the  profile  expressed  in 
t ,  4,  z  co-ordinates.  { t  is  measured  along  the  axis  of 
the  resultant  jet  with  its  origin  at  the  power  jet  nozzle  exit, 

4  is  measured  perpendicular  to  [  and  z,  and  t,  {,  and  z 
form  a  right  hand  orthogonal  co-ordinate  system.)  The 
question  is  then,  "How  should  the  u(x,  y)  be  translated  into 
>  4)  The  following  technique  was  adopted  for  the 
computation  of  f  ,  4).  First,  u^^Jx,  y)  was  computed. 

Then  a  Taylor  aeries  was  employed,  with  only  the  first 
two  terms  being  considered,  giving 

ujf.  5)  =  uj[x,  y"^  +  -|p  ujx,  y*)  At 

Note  that  At  =  0 

The  term  u^x,  y*)  was  evaluated  by  considering 

the  ujx,  y*)  to  behave  as  a  two-dimensional  jet  (i.  e. 
the  mean  velocity  follows  a  Gaussian  distribution) .  The 
center  line  velocity  decay  was  assumed  to  be  equivalent  to 
the  decay  relation  from  the  previously  obtained  oingle 
bounded  jet  data.  The  momentum  flux  thickness  was  avail¬ 
able  from  the  calculated  value  of  b  in  an  earlier  step  of  the 
program.  The  details  of  the  computational  procedure  are 
presented  in  Appendix  D. 
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The  two  primary  objectiveu  of  the  detailed  flow  at  idy 
are  (i)  to  determine  the  character  of  the  mean  flow  in  the 
resultant  jet  and  (ii)  to  determine  the  nature  of  the  inter¬ 
action  process  by  inference  from  the  characteristics  of  the  re¬ 
sultant  jet. 

The  character  of  the  resultant  jet  can  be  inferred 
from  several  comparisons.  Three-dimensional  effects  are 
indicated  by  differences  in  the  velocity  profiles,  the  flux  ratios, 
and  the  centerline  velocity  variation  (with  x)  as  a  function  of  z. 
The  character  of  the  resultant  jet  can  also  be  assessed  by  a 
comparison  of  pertinent  itemn  with  the  single  bounded  jet. 

The  role  of  the  vortex  stretching  will  be  very  impor¬ 
tant  in  determining  the  nature  of  the  resultant  jet. 

Three  basic  tets  of  results  are  presented  for  this 
section:  they  are: 

(i)  the  complete  velocity  profiles,  see  Figures  41 

to  46, * 

(ii)  isobaric  contours  in  Figures  47  to  49, 

(iii)  characteristics  of  the  jet,  viz.,  centerline  velocity 
distribution  and  mass,  momentum,  and  energy  flux  ratios,  see 
Figures  50  to  53. 

The  velocity  profiles  are  presented  for  seven  z/  a 

and  three  x/  a  locations  and  are  shown  as  u(y  )/u.aid  u  (n  }/u 

'  c  c 

The  isobaric  contours  are  formed  from  given  values  of  the  non- 
dimensional  pressure,  (P  ^  -  P)  /  4  p  u^  .  These  are 

presented  as  a  function  of  z,  i.  e.  in  a  plane  normal  to 
the  mean  flow.  The  measurements  were  made  as  a  function 
of  y  and  z  but  the  value  of  — —  is  sufficiently  small  to  set 

♦  The  figures  that  present  the  results  of  the  detailed 
mean  flow  survey  may  be  found  at  the  end  of  this  section, 
page  13  1, 


see 


^  i|c 

P(y  )  =  P(4)  for  all  the  y  values  encountered.  The  dashed 
portions  of  the  contours  are  extrapolations  to  regions  where 
no  data  were  taken. 

In  the  presentation  of  the  results  for  the  jet  character¬ 
istics,  the  increase  in  mass  flow  rate  was  considered.  For 
example,  the  term  ( 1  +  (mfr)l)  was  used  instead  of  Mq 
to  account  for  the  increased  mass  flow;  Eq(1  +  (mfrji) 
was  used  instead  of  to  account  for  the  increased  energy 
flux.  Since  the  value  of  the  x-direction  momentum  flux  ratio 
was  used,  the  term  aq^  describes  the  initial  momentum  flux. 
The  profiles  at  x/  a  s  30  are  actually  not  sufficiently  complete 
to  yield  a  highly  accurate  measure  of  the.  flux  ratios  because 
they  do  not  extend  to  low  enough  u/ u^  values  to  consider  the 
integration  to  be  from  -  •  to  +«o.  However,  this  defect  is 
not  serious  since  the  outer  portions  of  the  profile  contribute 
only  a  small  portion  of  the  toul  integral.  This  defect  is  most 
serious  for  the  mass  flux  and  least  serious  for  the  energy 
flux  ratio. 


Discussion  of  Results 

From  an  examination  of  the  results  of  this  section, 
the  following  model  is  proposed  to  describe  and  account  for 
the  observed  phenomena. 

As  the  control  jet  flow  strikes  the  power  jet,  the 
portion  of  the  control  jet  adjacent  to  the  near  side  pocket  is 
directed  into  the  near  pocket.  As  Manion  and  Goto  (17) 
describe  it,  this  is  similar  to  a  jet  impinging  on  a  solid  wall 
where  half  the  normal  component  of  flow  goes  to  either  side; 


I 
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note  however  that  the  "solid  wall" 

“  analogy  n  only  used  for  the 

Th.  j 
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proposed  flow  pattern.  The  decelerating  agent  is  the  increased 


FIGURE  54  Definition  of  terms  used  in  jet  intersection  model 


pressure.  There  are  wo  other  direct  results  from  this 
splitting  of  the  control  jet  flow;  there  are  (i)  a  very  strong 
vortex  production  in  the  near  side  of  the  power  jet  and  (ii) 
initial  profiles  of  the  resultant  jet  which  are  skewed  to  the 
left  (+  y  )  .  The  result  of  the  vortex  production,  caused 
by  the  power  jet  -  returning  control  jet  velocity  distribution, 
will  probably  have  the  form  of  a  well  defined  vortex  sheet. 

The  initial  profiles  which  are  skewed  to  the  left  are  shown  in 
Figure  44.  The  Gaussian  error  curve  is  included  on  each  z/a 
plot  for  reference;  the  skewness  of  the  profiles  is  apparent 
by  comparison  with  the  symmetric  reference  curve. 

The  infiltration  process  by  the  control  jet  is  seen  in 
the  non-dimensional  velocity  profiles  for  the  resultant  jet  at 
x/a  =  15  and  30,  see  Figures  45  and  46.  The  original  skewed 


*  The  control  jet  flow  enters  the  defined  region 
from  the  left  control  jet.  A  distribution  is  skewed  to  the  left 
if  the  third  moment  is  positive  (with  the  co-ordinate  system 
of  the  present  study  where  +y  is  to  the  left) . 
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left  profiles  are  now  skewed  to  the  right  (-y  )  with  the  effect 
being  more  pronounced  at  x/  a  c  30  than  at  15,  The  rate  of 
infiltration  is  related  to  the  degree  of  skewness  with  respect 
to  the  downstream  distance  since  the  longitudinal  distance  is 
the  product  of  the  convection  velocity  (viz.,  some  fraction  of 
u  )  and  the  infiltration  time.^  The  presence  of  two  effects 
makes  the  net  effect  difficult  to  measure  quantitatively.  The 
+  y  side  of  the  jet  at  x/  =  30  is  closer  to  the  Gaussian  curve 
than  at  x/ a  c  15.  This  indicates  that  the  infiltration  process 
lias  passed  through  the  left  side  and  the  jet  is  returning  to  the 
form  of  a  single  bounded  jet.  Whether  the  infiltration  process 
was  completedbefore  x/  a  c  15  cannot  be  ascertained  from  the 
data  at  hand.  As  discussed  in  the  "Literature  Survey"  section, 
Olson  (22)  has  observed  that  the  resultant  jet  appears  to 
leave  the  interaction  region  from  a  point  downstream  of  the 
intersection.  He  attributes  this  to  a  time  delay  in  the  inter¬ 
action  process.  This  observation  and  the  proposed  flow  model 
are  compatible. 

That  the  convective  effects  of  the  vortex  stretching 
are  much  stronger  for  the  resultant  jet  can  be  seen  by  an  examina¬ 
tion  of  Figure  50.  Although  the  difference  in  magnitudes  of  the 

u  /  u,^  values  is  slight,  there  is  a  pattern  wherein  u  is  less 
c  0  m 

and  u^  at  z/  a  ^  +2  is  greater  when  the  resultant  jet  is  compared 
with  the  single  bounded  jet. 

This  stronger  vortex  stretching  effect  may  be  ex¬ 
plained  by  the  flow  model  and  the  observed  profiles  of  u(y*/a). 

The  part  of  the  control  jet  that  turns  upstream  produces  a  very 

*  01s;n's  flow  system  is  geometrically  similar  to 
the  present  study  except  the  standoff  and  setback  walls  are 
removed. 
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strong  vortex  sheet  on  the  left  side  (+  y  ),  see  Figure  54. 
Also,  the  right  side  of  the  jet  has  a  steeper  than  usual  pro¬ 
file  because  of  the  turning  of  the  jet.  Both  of  these  effects 
give  the  vortex  stretching  action  a  stronger-than-nornnal 
start. 

That  the  interaction  process  is  inherently  three- 
dimensional  is  inferred  from  data  at  x/ a  =  5,  a  sUtion  which 
is  assumed  to  be  close  enough  to  the  intersection  zone  that  the 
intervening  flow  will  not  induce  strong  three-dimensional 
effects.  The  u(y  /  a)  profiles  show  that  the  velocity  profiles 
differ  with  respect  to  z  (Figure  41),  the  velocity  being  higher 
for  a  given  y  as  |  +  z  |  increases.  Isobaric  contours  at 
x/ a  =  5  also  show  a  strong  three-dimensional  character.  The 
longitudinal  distributions  of  Zbq^/aq^  and  E/ [l+(mfr)^^] 
are  also  indicative  of  a  three-dimensional  phenomenon  in  the 
interaction  region.  While  these  distributions  for  the  single 
bounded  jet  stort  together  and  then  diverge,  the  distributions 
for  the  resultant  j-t  are  divergent  between  x/ a  =  5  and  30. 

This  is  interpreted  as  a  developing  three-dimensional  effect 
for  the  former  case  and  a  fully  three-dimensional  effect  for 
the  latter. 

Several  pertinent  observations  may  be  made  with 
regard  to  the  overall  data.  The  entrainment  rate  of  the 
ambient  fluid  for  the  two  flow  cases  is  equivalent  although  the 
mass  flux  for  the  resultant  jet  is  less  than  expected.  The 
isobaric  contours  for  x/ a  =  15  and  30  indicate  that  this 
aspect  of  the  resultant  jet  flow  has  not  yet  attained  the  single 
bounded  jet  distribution. 
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FIG.  51  MASS  FLUX  RATIOS  FOR  THE 
RESULTANT  AND  THE  SINGLE  BOUNDED  JETS 
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FIG.  52  MOMENTUM  FLUX  RATIOS  FOR  THE 
RESULTANT  AND  THE  SINGLE  BOUNDED  JETS 
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NOZZLE  GEOMETRY  EFFECT  ON  THE  RESULTANT  JET 

Results 

The  data  collected  to  determine  the  effect  of  nozzle 

geometry  on  the  resultant  jet  is  not  intended  to  be  exhaustive; 

its  purpose  is  to  establish  trends  from  which  the  effect  may 

be  inferred.  Table  8*  shows  the  values  of  u{y  =  0)  /  u^  for 

X/  a  =  2,  Y/  a  =  1,  2,  3  and  z  =  0,  +2.  Profiles  of  u(y/  a)  /  u 

*  c 

for  z/  a  =  +2.  x/  a  =  1  5.  X/  a  =  2,  and  Y/  a  =  1,  2,  3  at 
an  mfr  =  0  are  presented  in  Figure  55.  A  dimensional-like 
plot  for  y  was  chosen  so  as  net  to  mask  any  effects  by 
non-dimensionalization  (e.g.  u(n)).  Table  9  gives  the  three 
flux  ratios,  mass,  momentum,  and  energy,  as  determined 
from  the  above  velocity  profiles;  Table  10  presents  the  static 
pressures  that  existed  on  the  standoff  pocket  walls. 

Discussion  of  Results 

In  the  initial  evaluation  of  the  data  that  were  collected 
to  determine  the  nozzle  geometry  effect,  it  was  noted  that 
the  velocity  profiles  for  Y/a  =  2  did  not  fall  between  the 
profiles  for  Y/a  =  1  and  3.  The  pattern  for  1,  3  and  •<> 
did,  however,  follow  a  consistent  pattern.  A.i  alternative 
way  to  express  this  unexpected  phenomenon  is  by  the  mass, 
momentum  and  energy  flux  ratios  as  shown  in  Table  9. 


*  Table  8  may  be  found  on  page  151. 
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FIG.  55  u(yVa)/Ug  FDR  z/o  =  *2 ,  x/a  =  15  ,  AND 


X/a=2  ,  Y/o  =  1,2,3,  ©  TO  SHOW  THE 
EFFECT  OF  THE  NOZZLE  GEOMETRY 


Table  9 

Mass,  Momentum,  and  Energy  Flux  Ratios 
To  Show  the  Effect  of  the  Side  Nozzles, 
x/a  =  15,  z/  a  =  +2,  X/ a  s  2 


M/  M 

o 

2bq  / aq 
c  ^o 

E/  E 

o 


Y/a 


1 _  2  3 


2.  47 

2.05 

2,38 

2.  180 

1.  23 

1.  05 

1.  22 

1.  12 

0.  706 

0.  601 

0.  711 

0.656 

The  tabular  data  come  from  the  velocity  profiles  and  are 
therefore  not  independent  of  them.  Consequently,  these 
flux  ratios  also  exhibit  an  extremum. 

A  possible  flow  model  could  be  formulated  on  the 
basis  of  the  static  pressure  data;  however,  these  data 
showed  a  consistent  trend  with  Y,  viz.,  a  larger  Y 
value  leads  to  a  higher  pressure,  and  this  made  it  appear 
that  some  combination  of  effects  was  causing  the  unex¬ 
pected  extremum. 

Table  10 

Values  of  Static  Pressure  at  he  Pocket  Walls 
For  Zero  mfr  Conditions 


2^ 


Y/a 


1 

2 

3 

oo 

-2.  49 

-  1.  13 

-0.69 

0 

Therefore,  it  was  necessary  to  develop  a  flow  model 


which  would  be  compatible  with  the  monotonically  decreasing 
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pressure  and  which  would  explain  the  observations  for  Y/a  =  2. 
It  was  reasoned  that  if  an  effect  were  present  which  altered  the 
vortex  stretching  action,  this  could  explain  the  apparent  anomoly 
for  Y/  a  =  2, 

Consider  the  defined  region  and  the  power  jet  flow 
only.  The  existence  of  a  low  pressure  in  the  pocket  implies 
a  standing  vortex  with  a  rotation  such  that  the  vortex  flow 
is  counterclock-w  se  in  the  left  pocket  and  clockwise  in 
the  right  pocket,  see  Figure  5  6.  The  flow  of  the  standing 


FIGURE  56  Flow  model  for  the  nozzle  geometry  effect 
vortex  past  the  forward  wall  of  the  power  jet  nozzle  will  form 
a  boundary  layer  on  this  surface.  As  a  result  of  this,  the  en¬ 
trained  flow  will  contain  vorticity  (generated  in  the  boundary 
layer)  with  its  major  component  in  the  z  direction  (-z  for 
the  left  side  pocket  and  +z  for  the  right  pocket) .  When  this 
vorticity  is  added  to  the  power  jet,  the  net  effect  is  to  reduce 
the  vortex  stretching  actio.->  because  the  vorticity  that  is  in¬ 
troduced  is  of  the  opposite  sense  from  that  generated  in  the 
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shear  layer  of  the  jet.  Since  this  effect  tends  to  retard  the 
development  of  the  three-dimensionality,  it  is  termed  a 
negative  effect.  An  important  feature  of  the  phenomenon 
is  that  the  magnitude  of  the  effect  will  be  dependent  on  the 
length  of  the  setback  wall  (Y)  and  the  strength  of  the  stand¬ 
ing  vortex  (which  can  be  qualitatively  inferred  from  the 
pocket  pressure) .  The  length  of  the  setback  wall  will  de¬ 
termine  the  development  of  the  boundary  layer  vortex  fila¬ 
ments  and  the  strength  of  the  standing  vortex  will  determine 
the  shear  rate  near  the  wall.  Therefore,  the  strength  of 
the  entrained  vorticity  is  fixed  by  a  combination  of  two  ef¬ 
fects.  Consequently,  the  vorticity  that  is  produced  in  the 
boundary  layer  will  be 

i)  small  for  very  large  Y/X  values  because  the 
standing  vortex  will  be  relatively  weak, 

ii)  small  for  very  small  Y/X  values  because 
the  boundary  layer  will  be  quite  small. 

iii)  attain  a  maximum  for  some  Y/  X  and  approach 
the  two  extremes  asymptotically. 

A  second  feature  of  the  flow  model  results  from 
the  low  pocket  static  pressures.  The  reduced  static 
pressure  in  the  pockets  will  cause  an  early  spreading  of 
the  jet.  However,  it  will  spread  more  rapidly  near  the 
plates  because  the  boundary  layer  fluid  has  less  x  -  di¬ 
rection  momentum:  therefore,  this  fluid  is  influenced  more 
by  the  transverse  pressure  gradient  than  the  fluid  in  the 
central  portion  (i,  e.  +z  <  2.75),  This  effect  induces  an 
initial  velocity  profile  which  resembles  a  profile  resulting 
from  the  vortex  stretching  action;  consequently,  the  proposed 
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eifect  is  a  positive  effect.  Figure  57  is  an  arbitrarily  drawn 
graphical  representation  of  the  results  of  these  two  effects 
on  a  typical  jet  characteristic;  M/M^  is  used  aj  the  or¬ 
dinate.  The  algebraic  sum  of  the  two  curves  can  be  inferred 
from  the  results  of  Table  9:  the  solid  curve,  which  is  the 
algebraic  sum  of  the  two  contributions,  was  drawn  from  the 
tabular  data. 

Although  the  above  analysis  is  not  quantitative, 
it  offers  a  flow  model  and  a  plausible  starting  point  for 
further  investigations. 

Two  independent  sets  of  data  were  collected  to 
establish  the  effect  of  the  nozzle  geometry:  the  velocity 
profiles  at  x/a  =  15,  z/a  =  +2  as  discussed  above  and 
four  groups  of  u^/ug  values.  The  purpose  of  the  U|./uq 
data,  presented  in  Table  8,  was  to  establish  whether  the  effect 
of  the  nozzle  geometry  (if  any)  was  of  such  a  nature  that  it 
would  be  apparent  in  the  u^  val  ues. 

After  careful  study  of  the  results  shown  in  the  table 
(for  x/a  =  10  or  less),  it  was  concluded  that  if  a  meaningful 
pattern  was  indicated  by  the  data,  it  was  both  very  complex 
and,  as  a  result  of  too  few  data,  not  accessible  to  direct 
analysis.  However,  it  was  possible  to  determine  that  if  the 
effect  existed,  it  was  of  a  reasonably  small  magnitude  since 
the  spread  of  the  values  for  the  different  geometries 

is  small. 

For  x/a  =  25  an  interesting  pattern  is  indicated. 

For  the  three  geometries,  Y/a  =  1,  2,  3,  the  difference  in 
the  u^/uq  values  is  (approximately)  0.  04  and  0.  02  when  the 
z/a  =  +2  and  z/a  =  -2  values  are  compared  with  the  values 


1^0 


FIG.  57  M/Mo  VERSUS  k  TO  INDICATE  THE 
RELATIVE  IMPORTANCE  OF  THE  TWO  TERMS 
IN  THE  PROPOSED  FLOW  MODEL  FOR  THE 
EFFECT  OF  THE  SIDE  NOZZLES 


151 


TABLE  8 


Uc/uq  to  Show  the  Effect  of  the  Side  Nozzles 


X/a 

Y/a 

x/a 

z/a 

5 

10 

15 

25 

2 

1 

0.  993 

0.864 

0.  705 

0.  523 

+  2 

2 

2 

0.  99 

0.860 

0.  71 

0.  541 

+  2 

2 

3 

0.  988 

0.860 

0.  72 

0.  540 

+  2 

- 

00 

0.  995 

0.  870 

0.  716 

0.  552 

+  2 

2 

1 

0.  975 

0.830 

0.  6  78 

0.481 

0 

2 

2 

0.  980 

0.810 

0.  674 

0.  497 

0 

2 

3 

0.  983 

0.825 

0.689 

0.  502 

0 

- 

00 

0.  988 

0.840 

0.  705 

0.  531 

0 

2 

1 

0.  971 

0.826 

0.  668 

0.  50 

-  2 

2 

2 

0.  995 

0.835 

0.679 

0.  518 

.  2 

2 

3 

0.  987 

0.  846 

0.  686 

0,  519 

.  2 

- 

00 

0.  997 

0.843 

0.686 

0.  535 

.  2 

For  X/a=  2 

Y/a 

x/a 

z/a 

1 

2 

3 

ee 

+  2-1/2 

0.  975 

0.  993 

0.982 

0.  971 

5 

+  1 

0.  976 

1.0 

0.  987 

0.  985 

5 

-  1 

0.  972 

0.  993 

0.  986 

0.  995 

5 

-  2-1/2 

0.  975 

- 

0.  984 

0.  994 

5 

+  2-1/2 

0.  822 

0.  826 

0.  804 

0.808 

10 

+  1 

0.  843 

0.  863 

0.835 

0.850 

10 

-  1 

0.  821 

0.  842 

0.840 

0.848 

10 

-  2-1/2 

0.804 

- 

0.  796 

0.804 

10 

16*, 


at  z/a  =  0.  In  contrast,  for  Y  =  «>,  the  'values  are  0.  017 
and  0.  004  for  the  same  comparison.  Another  mode  of  com¬ 
parison  is  offered  by  the  ratios 

u  /u_  z/a  =  +2 
c  0 _ 

u^/uq  z/a  =  0 

and  similarly  for  z/a  =  -2.  The  value  of  the  ratios  for  Y  =  1, 

2,  and  3  are  virtually  equal  and  the  value  for  Y  =  *>  is  smaller. 

From  these  observations  alone,  the  conclusions  would 
bo  that  there  is  no  strong  effect  of  the  control  jet  nozzles  before 
x/a  =  25  and  for  this  x/a  location  the  nature  of  the  effect  is 
essentially  independent  of  k  for  k  less  than  or  equal  to  1.5. 

Whereas  the  velocity  profiles  showed  a  pronounced  extremum 

for  Y/a  =  2,  the  u^/uq  values  are  quite  similar.  Two  con¬ 
clusions  are  drawn  from  these  observations;  they  are,  what¬ 
ever  the  full  effect  of  the  nozzle  geometry:  i)  it  has  a  strong 
memory  characteristic  since  the  same  nozzle  geometries 
induced  two  different  types  of  effects  at  x/a  stations  of  15 
and  30;  ii)  it  is  very  complex. 
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MOMENTUM  FLUX  RATIO  EFFECT  ON  THE  RESULTANT  JET 

Results 

As  in  the  section  on  the  nozzle  geometry  effect,  IX, 

experiments  were  conducted  to  provide  trends  from  which  the 

full  effect  of  the  mfr  on  the  resultant  jet  may  be  inferred. 

One  geometry  and  three  values  of  the  mfr  were  investigated. 

Figures  58  and  59  present  the  profiles  for  u(y  /a)/u  at 

c 

z/ a  c  0  and  +2  and  for  the  mfr  values  of  0,  0.  1,  0.  2.  As 
another  mode  of  comparison,  the  profiles  of  u.  (4/a)/u 
and  u^(t}j^)/u  are  presented  in  Figures  60  to  6  3  for  z/ a  *  0 
and  +2  and  for  mfr  values  of  0,  0.  1,  0.  2,  Table  11  shows  the 
mass,  momentum,  and  energy  flux  ratios  for  the  resultant  jet, 
for  the  z/  a  and  mfr  parameters  as  above,  and  for  the  x,  y, 
z  and  t,  4i  z  co-ordinates  respectively. 

The  two  sets  of  figures  (z/  a  =  0,  +2),  each  containing 
three  plots,  were  considered  necessary  for  a  complete  pre¬ 
sentation  of  the  data  for  the  following  reasons: 

i)  The  measured  quantities,  viz.,  0,  h  ,  and  h  , 

T  s 

lead  directly  to  the  u(y*/  a)  /  u^  profiles.  Also,  since  the 

output  ports'  configuration  can  be  most  easily  described  in 

♦ 

X,  y,  z  co-ordinates,  the  u(y  )  profiles  provide  the  most 

meaningful  description  of  the  flow  into  these  ports. 

ii)  The  u^(4/a)/u  profiles  present  the  velocity 
^  * 

distribution  in  the  natural  co-ordinates  of  the  flow.  viz. ,  f  , 

4,  z.  They  also  offer  a  meaningful  comparison  with  the  u{y/ a)  /  u 
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profiles  of  a  single  bounded  jet  as  well  as  the  most  meaningful 
form  that  will  indicate  profile  variations  with  respect  to  z. 

iii)  The  profiles  of  present  the  data  in  such 

a  form  that  a  comparison  with  the  two-dimensional  jet  is 
possible.  Also,  since  the  Gaussian  curve  can  be  plotted  on  these 
co-ordinates,  the  skewness  of  the  profiles  is  easily  seen. 


Discussion  of  Results 

A  striking  feature  of  the  velocity  profiles  is  the 
regular  pattern  for  -4/ a,  viz.,  a  pattern  described  as 
"an  increasing  u_^  with  increasing  mfr  for  a  fixed  4  /  a" 
and  the  irregular  pattern  for  +4/  a.  However,  this  may  be 
explained  by  notipg  that  the  mass  flux  is  increased  for 

®  0.  2  and  the  rate  of  infiltration  is  also  increased  for 
this  higher  mfr  with  respect  to  an  mfr  of  0.  1. 

In  the  "Detailed  Mean  Flow  Study"  section 
it  was  noted  that  the  profiles  changed  from  a  skewed  left  to 
a  skewed  right  pattern  between  x/ a  =  5  and  15.  From 

Figures  62and  63  it  ;8  apparent  that  the  profiles  are  skewed 

>(< 

right  (-y  )  at  x/  a  s  15;  furthermore,  the  skewness  is  more 

pronounced  for  z/ a  =  +2  than  for  z/a  =  0.  This  difference 

is  assumed  to  result  from  the  vortex  stretching  effects. 

With  reference  to  Table  11,  the  mass  flux  ratio  drops 

as  the  mfr  increases  from  0  to  0.  1  for  z/a  =  0;  however, 

the  mass  flux  ratio  has  the  opposite  trend  for  z/a  s  +2. 

Quantitatively,  the  ratio 

mass  flux  ratio  at  a  mfr  of  0 
mass  flux  ratio  at  a  mfr  of  0.  1 
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TABLE  1 1 

Values  of  the  Flux  Ratios  to  Show  the  Effect  of  the 
mfr  (x,  y,  z  and  f  ,  4  ,  z  co-ordinates),  x/a=15,  k=1.0 


mfr 


z/a 

0.  0 

0.  1 

0.  2 

M/MQ(l  +  [mfr] 

0 

2. 23 

2.  02 

2.  24 

M#/Mo(l  +  [mfr]  1/2) 

0 

2.  23 

2.  02 

2.  24 

M/MQ(l+[mfr]  1/2) 

+  2 

2.  05 

2.  34 

2.  50 

M#/Mq(1+  [mfr]  1/2) 

+  2 

2.  05 

2.  34 

2.  51 

2  b  q^/a  qo 

0 

1.21 

1. 26 

1.  19 

2  b  q^/a  qo 

+  2 

1.  05 

1. 30 

1.  34 

E/Eo(l  +  [mfr]3/2) 

0 

0.  726 

0.  568 

0.  654 

E*/Eo(l  +  (mfi]3/2) 

0 

0.  726 

0.  570 

0.  655 

E/Eo(l  +  [mfrJ  ^/2) 

+  2 

0.  601 

0.  746 

0.  754 

E*/Eo(l  +  [mfr]3/2) 

+  2 

0.  601 

0.  747 

0.  760 

NOTE; 


The  momentum  flux  ratio  for  the  streamline  co-ordinates 
is  not  tabulated  because  there  is  no  reliable  way  to  generate 
a  reference  value. 


has  the  values  1.  10  and  0.  876  for  z/  a  values  of  0  and  +2 
respectively.  These  values  imply  the  existence  of  a  strong 
mfr  effect  on  the  vortex  stretching  since  increasing  the  mfr 
causes  a  shift  of  the  mass  flux  toward  the  walls.  The  same 
ratio  formed  for  mfr  values  of  0.  1  and  0.  2,  has  values  of 
0.  90  and  0.  935  respectively  for  z/a  =  0  and  +2.  This  is  not 
believed  to  be  a  rejection  of  the  conclusion  that  the  mfr  level 
has  an  effect  on  the  vortex  stretching;  rather,  it  is  assumed 
to  be  indicative  of  an  increase  in  the  infiltration  rate  with 
increasing  mfr.  Two  factors  support  this:  i)  the  mass  flux 
ratio  increases  with  increasing  mfr  at  z/ a  =  +2;  ii)  the 
velocity  profiles  in  Figures  62  and  63  show  that  the  profile  is 
more  symmetric  (i.  e.  greater  infiltration)  for  z/  a  =  0  than 
for  z/a  *  +2. 

It  is  concluded  that  the  mfr  does  exert  a  marked 
effect  on  the  resultant  jet;  however,  the  character  of  this 
effect  is  difficult  to  specify  from  the  available  data  because  there 
is  an  interplay  between  the  vortex  stretching  and  the  infil¬ 
tration  rate  effects. 
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CONCLUSIONS  FROM  THE  STUDY  OF  INTERSECTING  JETS 

The  following  conclusions  are  supported  by  the  results 
of  the  study  on  intersecting  jets  for  an  aspect  ratio  of  six  and 
three  nozzle  geometries. 

1)  For  the  range  of  setback  to  standoff  wall  ratios 
of  0.  5  to  1.  5,  the  smallest  ratio  provides  the  conditions  for 
the  highest  gain,  viz. ,  the  largest  deflection  angle  and  the 
farthest  upstream  location  of  the  apparent  pivot  point. 

2)  The  analytical  model  proposed  provides  a  satis¬ 
factory  prediction  method  for  the  deflection  angle  for  the 
geometry  of  the  smallest  ratio  of  setback  to  standoff  wall 
lengths. 

3)  The  analysis  of  Manion  and  Goto  (17)  provides 
a  satisfactory  prediction  of  the  deflection  angle  for  the  geo¬ 
metry  of  the  largest  ratio  of  setback  to  standoff  wall  lengths. 

It  should  oe  noted  that  this  analysis  requires  an  apriori 
knowledge  of  the  far  side  pocket  pressure. 

4)  Detailed  mean  flow  surveys  of  the  resultant  jet 
for  one  momentum  flux  ratio  and  one  nozzle  geometry  indicate 
that  the  effects  of  the  vortex  stretching  are  increased  following 
the  interaction  of  the  jets  as  compared  with  the  single  bounded 
jet. 

5)  The  nozzle  geometry  and  the  momentum  flux  ratio 
have  marked  and  complex  effects  on  the  character  of  the  vor¬ 
tex  stretching  and  on  the  infiltration  rate  of  the  control  jet  into 
the  power  jet. 
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APPENDIX  A  -- 

NOZZLE  AND  CONTRACTION  DESIGN 


In  order  to  provide  for  a  smooth  acceleration  of  the 
fluid,  a  beam  deflection  curve,  pictorially  defined  in  Figure 
A  1,  was  employed  for  all  nozzles  and  contractions  used  in 
the  present  study,  A  smooth  acceleration  curve  is  a  necessary 
condition  to  produce  a  uniform  efflux  velocity  from  the  nozzle. 
An  analysis  is  presented  below  which  esUblishes  the  geometric 
conditions  that  are  necessary  for  a  smooth  acceleration. 

Consider  the  one-dimensional,  irrotational,  and  in¬ 
compressible  flow  of  an  inviscid  fluid  through  the  nozzle  of 
Figure  A  2  whose  cross-sectional  area  (A)  is  given  by  the 
expression 


A  =  W  -  2g(x) 

For  steady  flow,  the  convective  acceleration  is  V{dV/dx). 
The  continuity  of  d*  V/  dx*  will  guarantee  a  smooth  accel¬ 
eration  of  the  fluid  through  the  nozzle.  After  rewriting  the 
continuity  equation,  p  AV  ■  constant,  in  the  form 

dv 


1 

A 


dA 

dx 


V  dx 


d^  V/  dx^  can  be  expressed  as 
d^V  V/  ^ 

jTT  *  -  a 

dx‘ 


d;^^  ^  'a^'  'dT' 


Therefore,  the  requirement  that  d^V/dx*  be  con¬ 
tinuous  can  be  replaced  by  the  same  requirement  for  d^  A/  dx^ 


DEFLECTION  CURVE 


FIG.  A2 


DEFINITION  OF  TERMS  USED 
IN  APPENDIX  A 
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Since  A  is  proportional  to  g(x),  the  continuity  of  d^V/dx* 
is  also  satisfied  by  the  continuity  of  d*  g(x)  /  dx*  .  But,  this 
condition  is  guaranteed  for  the  beam  deflection  curve  since 
the  moment,  which  is  related  to  the  term  d*  g{x)  /  dx^  ,  is 
continuous. 

Because  of  construction  simplifications,  a  circular 
arc  IS  commonly  used  for  the  curve  g(x) ,  i.  e.  g{x)  = 

R  sin  (cos  ^  x/R).  However,  since  the  radius  of  curva¬ 
ture,  which  is  related  to  d^g(x)/dx*,  is  discontinuous  at 
X  «  0,  the  fluid  (under  the  restrictions  of  the  analysis)  will 
undergo  an  instantaneous  change  in  its  rate  of  acceleration. 

Since  this  is  physically  impossible,  the  fluid  will  adjust  itself 
in  such  a  manner  that  the  flow  streamlines  will  not  follow 
similar  paths  over  the  cross-sectional  area,  and  this  ad¬ 
justment  will  produce  a  non-uniform  flow  at  the  nozzle  exit. 

The  nature  of  this  non-uniformity  cannot  be  predicted  by  this 
analysis,  since  the  one-dimensional  assumption  and  a  non- 
uniform  profile  are  incompatible. 

The  magnitude  of  this  phenomenon  is,  of  course,  de¬ 
pendent  on  the  magnitude  of  R  with  respect  to  the  nozzle  exit 
width;  the  non-uniformity  becomes  less  pronounced  with  in¬ 
creasing  values  of  the  ratio  R/(W  -  2R),  A  pronounced  "dog 
eared"  velocity  profile  results  when  this  ratio,  R/ (  W  -  2R), 
is  too  small.  It  is  important  to  note  that  the  above  presentation 
is  quantitatively  true  only  for  the  assumed  inlet  flow  conditions 
and  constitues  no  more  than  a  statement  of  necessary  (but  not 
sufficient)  conditions.  If,  for  example,  a  nozzle  were  designed 
which  satisfied  the  above  requirements  but  incorporated  too  small 
a  radius  of  curvature;  then,  for  a  real  fluid,  the  flow  pattern  could 
also  lead  to  a  "dog  eared"  profile  because  of  high  local  accelera¬ 
tions  near  the  walls. 
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APPENDIX  B  -- 

PROBE  CALIBRATIONS 

Three  probes,  the  forward  facing  12  hole  static,  the 
tapered  toUl,  and  the  three  tube  trident,  were  calibrated  for 
the  effects  of  yaw.  Since  the  same  basic  calibration  technique 
was  employed  for  the  three  probes,  the  general  method  of  test¬ 
ing  will  be  described  and  then  the  data  and  analysis  for  each 
probe  will  be  given  separately. 

The  probe  to  be  calibrated  was  mounted  on  a  vertical 
support  rod  attached  to  the  x  -  y  -  0  device.  The  device  was 
located  approximately  2  diameters  below  the  exit  of  an  8  inch 
O.  D.  pipe;  the  probe  tip  was  on  the  center  line,  one  diameter 
from  the  exit  of  the  pipe.  A  calibration  run  consisted  of 
measuring  the  dynamic  head  (q)  and  the  probe  response  as  a 
function  of  its  angular  position  with  respect  to  the  flow. 

Static  Pressure  Probe 

Results  of  the  calibration  runs  for  this  probe  are 
shown  in  Figure  B  1.  In  order  to  develop  an  empirical  for¬ 
mula  for  Ah^  as  a  function  of  6,  three  regions  were 
specified  as: 

Region  I  -  ©  >  15  degrees,  all  values  of  q 

Reg  ion  11-6  <  15  degrees,  q  >0,4 

Region  III  -  0  <  15  degrees,  q  <0.4 

Since  all  the  Ah^  values  of  region  III  were  less  than 

0.  02  inches  of  water,  was  assigned  the  arbitrary  value 

of  0. 

A  second  order  polynominal  was  assumed  for  Ah  , 

s 
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such  that: 

Ah^  r  A(q)  e  *  +  B{q)  0 

From  the  three  data  runs  in  region  II,  A(q)  = 

(0.  16  q-  0.07)  and  3  a  -0.44. 

For  I,  the  relationship  between  Ah^  and  0  is 
linear,  and  therefore  can  be  described  as: 

Ah^  a  M(q)  0  +  N(q) 

Figure  B  2  shows  the  values  of  M(q)  and  N(q) 
plotted  for  the  values  of  q  tested.  The  formulas  for  these 
functions  are: 

Mr  [9q  -  0.3  5]  X  10 
N  a  [  5q  -  0.  18]  X  10 

Tapered  Total  Pressure  Probe 

Unlike  the  static  pressure  probe,  the  results  for  this 
calibration  are  accurately  expressible  as  a  single  curve  for  all 
q  values  tested.  Figure  B  3  shows  the  experimentally  deter¬ 
mined  points  and  the  assigned  calibration  curve  which  follows 
the  equation 

^  -  0  '  -  O  =  C(0-6.7)  +D(0-6.7) 

^ 

for  0  >  10.7 

and 

Ah 

_ T  r  0  for  0  <  10.7 

q 

The  values  of  C  and  D  were  evaluated  as  0.4  x  10 
and  -  1.  55  X  10  respectively.  This  formulation  yields  a  value 
of  Ah^  to  within  1  per  cent  of  the  measured  value  for  all  values 
of  0  and  q  investigated. 


Trident  Probe 


The  purpose  of  this  calibration  was  to  allow  a  dif¬ 
ferential  reading  of  the  trident  probe  to  be  used  as  a  measure 
of  the  yaw  angle  with  respect  to  the  probe  axis.  Since  the 
differential  reading  was  a  function  of  the  dynamic  head  as  well 
as  the  angle,  the  calibration  data  are  presented  as  plots  of  © 
vs.  Ah,j.  /q  with  q  serving  as  a  parameter.  These  data, 
which  are  presented  in  Figures  B  4  and  B  5  are  seen  to  follow 
straight  lines  passing  through  the  origin.  These  straight  lines 
may  be  represented  by  the  relation 

0  =  m(q)  [  Ah  /  q] 

45 

Figure  B6  indicates  the  variation  of  m{q)  with 
q.  Since  an  empirical  formula  that  could  be  adapted  to  a  com¬ 
puter  program  was  desired,  the  following  relationship  was 
chosen. 

m(q)  =  31  -  30  q  for  9<q  <0.  1 
m(q)  =  28  for  q  >0.  1 

The  accuracy  of  this  formulation  is  dependent  upon  the 
magnitude  of  Gas  indicated  by  Figures  B  4  and  B  5.  However, 
the  maximum  error  of  the  computed  angle  given  by  the  above 
formulation  will  be  of  the  order  of  5.  5  per  cent.  Because  of 
the  complicated  nature  of  the  computation  procedure,  no  genera 
statement  can  be  made  as  to  the  effect  of  this  error  on  the 
velocity  magnitudes  except  that  the  errors  will  not  be  large 
enough  to  become  unacceptable  (viz.  greater  tlian  about  5  per 
cent) . 


iillD4  «  C.I2 


FIG*  B2  EVALUATION  OF  CALIdRATlCN  PARAMETERS  *  STATIC  PROBE 
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FrG.  Q9  CALlBftATtOM  CUffVE  FOR  THE  TRIDENT  PROBE 


SYMBOL 
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APPENDIX  C  -- 

SOLUTION  OF  EQUATIONS 
FROM  THE  JET  DEFLECTION  ANALYSIS 


Eight  equations  involving  eight  unknowns  were  developed 
in  section  VII ;  this  appendix  presents  the  reduction  of  the  eight 
equations  to  two  computation  equation'^.  In  the  presentation 
below,  the  steps  will  be  referred  to  by  reference  to  the  equations 
used;  the  detailed  algebraic  steps  will  not  be  shown.  The  equa¬ 
tions  from  the  analysis  in  the  body  of  the  text  are  reproduced 
here  for  reference  with  the  exception  that  equations  2  and  3 
have  already  been  used  in  equation  4  in  order  to  eliminate 


the  variables  P,  and  P, 


Equation  Number 
4 


6 

7 

8 


Equation 


*  tqi/Rilx' 

^qi/Rz  -  qi/Ri 


]Y'  +q 


tan  a  s  qi  sin  y  q?  cos  g 
qz  cos  y  +  qi  ein  /3 

1  +  Y/Rl  «  cos  y  +  (Ri/  Ri)  sin  fi 

l+X'/Ris  (Ri/Rj)  siny  4  cos /} 

tana  =  Ri(l-cosy)/[Ri  siny-x  1 

app 

X  /X'  =  0.25  k  (Y)  +  0.875(X'X^ 
app  v' 


A  non-dimensionalizing  procedure  will  be  employed 

for  all  factors  as  given  by  the  relationships: 

qz/  qi  =  mfr,  Y/  X  =  k,  x  /  X*  =  d, 

app 

Rl/X'  =  r'i,  Rz/X  *  b!z 
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Equations  4  and  9  can  be  combined  to  solve  for  Rj 
in  the  form  of  R'a  =  R'j  (R'l ,  y) , 

mfr  ^  [1  +  R'l  mfr][Ri  sin  y  -  d]  -  R‘i^  (1  -  cos  y)  +  kR'i(l  -  cos  y) 
R't  2  R'l'  { 1  -  cos  y)  k  ' 

(10) 

Solving  equations  6  and  7  for  sin  fi  and  cos  p  ,  then 
squaring  and  adding  the  resulting  equations,  results  in  another 
formulation  for  Rz  =  R'2(R'i  r  y) 

nifr  2  (R'l  sin  y  -  1)  mfr 

*  [Ri  (1  -  cos  y)  +  k]*  +  [1  -  R'l  sin  y]* 

(11) 

Equations  10  and  11  can  be  equated,  eliminating  R^, 
and  the  resulting  expression  can  be  reduced  to 

2{  1  -  cos  y)  [mfr  (sin  y)  -  (1  -  cos  y)] 

+R'l*{  4  sin  y  (1  -  cos  y)  -2  mfr  6(1-  cos  y) 

-  2  mfr  (sin  y  -  k(l  -  cos  y)  sin  y  } 

+R  l*  {2  6  [mfr  sin  y  -  k  ( 1  -  cos  y)  mfr  -  ( 1  -  cos  y)  } 

+R'i{[2k  ( 1  -  cos  y)  -  2  sin  y]  6 

+  [sin  y  -  mfr  6  +  k  ( 1  -  cos  y)  ]  (k*  +  1)  } 

+  (k*  +  1)  (-8)  «  0  (12) 

If  equations  6  and  7  are  solved  for  sin  p  and  cos  p  and 
if  these  factors  are  used  in  the  combined  equations  4  and  9,  a 
quadratic  expression  for  R'l  results. 

«  -B  +  -/[B*  -  4AC] 

2A 
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where 

A  *  2k  ( 1  -  R*|  Bin  y) 

B  =  {2k  (mfr  +  Bin  y)  +  (1  -  Bin  y)  (1  -  1/  R'|) 

-  [  1/  Ri  +  mfr  ]  [  k  +  R'l  ( 1  -  cob  y)  ]  } 

C  ■  (mfr  +  Biny)  (1  -  k/R'i)  -  cob  y(l/R'|  +  mfr) 

EquationB  11  and  13  can  be  combined  to  eliminate 
Ra,  thus  yielding  a  Becond  equation  for  R'l  and  y. 

(mfr)  *  k  8  ( 1  -  R'l  sin  y)* 

+  {(mfr  +  sin  y)  (1  -  k/R'i)  -  cos  y(l/R'i  +  mfr)} 

{[R'l(  1  -  COB  y)  +  k]*  -f  ( 1  -  R'l  sin  y)  * }  * 

-2mfr(l  -  R'l  sin  y)  {Zk  (mfr  +  sin  y)  'f  ( 1  -  R'l  sin  y)  (1  -  l/R'i) 
-  (  1/R'i  -  mfr)  [k  +  R'l  ( 1  -  cos  y)  ]} 

*{[R'l  ( 1  -  COB  y)  +  k]*  +  [  1  -  R'l  sin  y]*  }  ■  0  (14) 

The  computational  procedure  which  was  chosen  from 
this  point  was  dictated  by  the  complicated  algebraic  equations. 
Using  an  IBM  7090  for  computation,  equation  12  was  solved 
for  a  given  value  of  y.  With  the  R'l  from  this  computation  the 
left  side  of  14  was  evaluated  and  its  value,  as  well  as  R'f,  was 
printed  as  output.  The  value  of  R't  and  tana  were  then  compu¬ 
ted  from  equations  11  and  4  and  the  resulting  value  printed  ar 
output.  Finally,  two  non-dimensional  representations  of  the 
pressure  were  formed  as  (using  equations  2  and  3): 


»  mfr 

(2a) 

qi 

^Pr-  ^'Pr 
qi 

■  mfr 
R‘l 

-  1 

R'l 

(3a) 

These  values  were  also  computed  and  printed  as  the  final 
output.  This  procedure  was  repeated  until  all  the  values  of 
y(l  <y  <42)  were  employed,  then  the  entire  procedure  was 


repeated  until  all  of  the  valuea  of  k  and  the  mfr  were  used 
(viz,,  k»0.5,  1,0,  1, 5  and  mf r  0, 05,  0,  10,  0,  15,  0,20, 
0,25), 

Since  equation  12  ia  a  fourth  order  equation,  there 
are  four  roots;  however,  only  one  of  these  is  a  physically 
possible  solution.  Consequently,  the  criterion  to  select  the 
correct  value  of  <y  was  established  as  that  y  which  led  to  a 
zero  value  of  equation  14  and  positive  values  of  R't  and  R't. 
Linear  interpolations  were  employed  where  necessary  in 
the  evaluation  of  the  computer  print  out. 


APPENDIX  D  — 

PROGRAM  FOR  y*  TRAVERSES 
The  Fortran  IV  program  which  was  used  for  the  data 
reduction  on  an  IBM  7090  is  presented  below. 

Method  1  (direct  measurement  of  6)  and  Method  2, 
as  described  introduced  in  the  "Detailed  Mean  Flow  Survey" 
section,  VIII,  required  a  different  format  for  the  data  cards 
and  a  change  of  four  Fortran  statements.  The  four  chiinges 
for  Method  2  are  listed  on  the  program  statement  to  the  right 
of  the  regular  Method  1  statements.  Static  pressure  and  total 
pressure  corrections  were  made  in  accordance  with  the  for¬ 
mulas  determined  in  Appendix  B;  note  that  these  corrections 
are  not  iterative.  The  calculation  for  the  angle  from  the 
A  h^  data  (in  Method  2)  was  applied  twice.  This  calculation 
actually  should  be  an  iterative  process,  but  there  was  no  way 
to  generate  a  goodness  criterion.  Repeated  application  of  the 
computation  would  lead  to  an  unbounded  6  value;  consequently, 
the  arbitrary  choice  of  two  repetitions  was  made.  The  pro¬ 
gram  was  also  written  such  that  the  maximum  (or  minimum} 
angle  could  be  +  80  degrees.  This  was  necessary  because 
certain  situations  arose  involving  a  very  small  q  while 

A  h  was  large  enough  to  yield  unrealistically  large  angles 
^  45 

(e.g.  a  q  value  of  0.001  inch  of  water  and  a  Ah  value  of 

^45 

0.  004  inch  of  water  would  lead  to  a  6  of  12  1  degrees) . 
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A  complete  traverse  was  taken  at  z/a  e  0,  at  the 

three  stations  x/ a  =  5,  15,  30,  and  for  a  mfr  =  0.  1  using 

Method  1  and  Method  2.  A  comparison  of  the  results  from 

these  six  traverses  gives  an  indication  of  the  accuracy  of 

Method  2.  Figures  D- 1  to  D-3  show  the  velocity  profiles 

(u/  u  vs  y  )  and  Figures  D-4  to  D-6  show  the  values  of 
♦  * 

0  vs  y  .  A  plot  of  u/  u  (from  Method  1)  is  included 

c 

on  the  figures  to  provide  a  reference  for  the  dimensions  of 

the  jet.  The  agreement  shown  by  the  velocity  profiles  is 

quite  satisfactory.  The  agreement  shown  by  the  two  methods 
* 

for  ©  is,  in  general,  acceptable  except  for  x/  a  e  5  and 
y/  a  less  than  about  -0.  75.  Since  the  principal  purpose  of 
the  data  of  Method  2  was  the  establishment  of  mean  velocity 
profiles,  the  method  was  satisfactory. 

The  variables  used  in  the  text  of  this  writing  were 
necessarily  changed  into  Fortran  Language  variables;  the 
text  and  the  equivalent  program  variables  are: 


Program 

Text 

Program 

Text 

Program 

Text 

YS 

y* 

HT 

T 

HS 

h 

'  s 

DIAL 

t 

DELHT 

T 

46 

XOA 

x/a 

ZOA 

z/  a 

DC 

X 

lY 

Y 

FMR 

mfr 

ZETA 

t 

UO 

“o 

TEMP 

T 

PRESS 

B 

TC 

e 

C 

QO 

*10 

HTC 

•'t 

HSC 

h 

8 

C 

DIAL  0 

^0 

CP 

Cp 

Q 

q 

VE 

V 

U 

u 

V 

V 
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Program 

Text 

Program 

Text 

Program 

Text 

Uouc 

u/  u 

c 

UOUCSQ 

(u/u  )* 
c 

VOUC 

(v/  u^i 

US 

U* 

T 

e 

XI 

i 

USOUC 

USOUCS 

YOB 

y/  b 

XIOBS 

4/b, 

AK 

K 

EM 

m 

B 

b 

E 

E 

EMOMO 

M/  Mq 

RATIO 

2bqc/qo 

EOEO 

E/Ejj 

EMS 

BS 

ES 

DUSDZ 

eu,(x, 

EMSOMO 

RATIOS 

2b^q^/aqo 

ESOEO 

^,'^0 

Note  TT  is  used  in  the  connputation  for  6 

HSP  is  used  in  the  computation  for  h^. 


lyVa)  FOR  x/o=5,  z/a  =0  FROM  METHOD  I  AND  2  TO  SHOW 

the  effect  of  the  measurement  technique 


THE  EFFECT  OF  THE  MEASUREMENT  TECHNIQUE 


I 


1 

I 


FIG.  0-4  y  VERSUS  yVo  FOR  x/o=5,z/a=0,  METHODS  182 


Smaot.  bflTHOD 


FIG.  D5  U  VERSUS  yVa  FOR  x/o  =  15  ,  z /o=  0  ,  METHODS  182 


FIG.  06  0  VERSUS  yVa  FOR  x/o  =  30,  z/a=0  ,  METHODS 
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DATA  PACCESSING  FOR  THE  PRESENT  STUCV  V  STAR  TRAVERSES 
OtPEKSICN  TSISOli  CPiSOli  CISOli  VEiSOIt  UtSCIi  VI50I«  UOUCTSOIt 
lUClCSQISOt  VCUCTSOli  USCSOti  T(SOI»  V  ISCii  XICOTi  USOUC  (SOI  lUSO 
2t'C  S  I  SO  I .  TT  (  SC  I  t  YOB  I  SC  I .  X  100  S  (  SO  I 

120C  READ  IS,  2001  IXOA,  20A,  IX,  lY,  FHR,  IDAV,  INCNTH,  IVEAR,  NO 


2CC  FCPMAT  (1S,FS.2,  2IS,FS.2,  AI5I  .  . 

READ  IS,  2021  N,  {ETA,  UO.  lEHP,  PRESS 
202  FORNAT  1112,  4F12.S  I 

READ  IS,  20AI  TC,  00,  HTC,  FSC,  OIALO  READ  ISt  20AI  TC,  00,  HTC,  HSC 
2CA  format  ISF12.SI  204  FORMAT  I4F12.SI 

TC--TC  ^200  READ  IS.  206IVItl,HT,  HS.  DIAL,  OELHT, OIALO 

00  48  I ‘UN  2®*  format  (AFIO.OI  _ 

22CC  REAO  IS,  2C6IYIII,FT,  PS,  DIAL,  CELH1 
2C6  FCPMAT  ISFIO.OI 

IF  lOIALO)  SOU,  2,  SCO  .  . 

SOC  TSIl  >•ICIAL-CIAI.0I/C.164 
IS  01  I  I*  hT-MS 
6SC  AK«lS.e9*S0RTITEMP/PPESSI 


OC-MTC-HSC 

6S1  UC>AK*SCRTIOC>*CaSI  TC/S7.20C I 

ucs«ak,sgrticci 
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IF  (hSP-MSI  6,6,301 
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IF  lOlAlOl  620,800,620 
620  HS<HSP 
GC  TC  6 

4  IF  IQIII-0.41  6,6,302 

SC2  HSF>hS*  I0.16,cill>a.0ri*  IABSITSI  IMTCI  ••2.0-C.44*ABS(1SI||a 
1  Tcn*io.ooii 

IF  IHSP-HSI6,  6,  303 
3C3  OII)*MT-MSP 

IF  I0IA10I304,  800,  304 
304  HS'HSP 

6  IF  ICIALOl  3CS,  14,  SOS 
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IF  IQlin  47,76.76 
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UOLCIIIo  UIII/UC 

ucccsoiii*  ucuci 1 1*00001 1 1 
vccoiii*  VII  i/uc 

USIl  l•AK*SORT(CI II l•COS(TSI II/S7.2S6I 
Tl  I  l•TSln*TO 
GC  TC  f? 


47  VE(II«C. 
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OCLCI  11*0. 
OOLCSUI 1 1*0. 
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_USCUC(ll*US(l)/UCS 

uscucs(n>uscuc(  1  l••2.o 

00  36  J«2,N 

XI(JI»Y(Jl*CCS(rC/67.296l 

0US0Z«(UCS«0.0S67fl.(ABS(Y(JI/Bn»#2.0/B-1.14«UCi/7ETA»*l.5l« 
1  EXP(-0.392  7*( (A6S(Y(J|/B»  »**2.C) I 
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109  NKI Tr 16, 101 1  7  1 1  1 ,TS I  I  I >  T I  I  1 ,CI I l.CPlIl, VEI 1 1 ,UOUC  I  I  l.VOUCI I  I 

101  FUBPAT(6X,F9.3,3X,FS.2,3X,F  ...2,1PE10.2,1X,1PE1C.2,  1PE1C.2,9X, 

21Pt 10.2  ,9x,lPElo.2) 
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1  (M,  E,,l).2, 

2  ICH,  EMOP  *,  E9.2,  1  RATIO  *,  E9.2,  91-,  EOlO  *, 

3  E9.2  I 
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108  FORMAT  (///6P  UCS  -,  IPIS.2,  7E,  EPS  *,  E9.2,  EM,  BS  *,  E9.2. 

1  EH,  ES,,C9.2, 
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APPENDIX  E  -- 

y"'  TRAVERSE  RESULTS 

This  appendix  presents  the  computer  results  from 
the  37  data  runs  that  comprised  the  y*  traverses.  The 
other  traverse  data  {u(z)  )  is  implicitly  presented  in  Figures 
9  and  10. 

The  presentation  of  the  y*  traverses  involves  repro¬ 
ductions  of  the  original  results  sheets.  They  are  presented 
in  the  following  order; 

i)  "Single  Bounded  Jet" 

ii)  "Detailed  Mean  Flow  Survey"  x/  a  *  S,  15,  30 

iii)  "Nozzle  Geometry  Effect  on  the  Resultant  Jet" 

iv)  "mfr  Effect  on  the  ResulUnt  Jet” 

v)  Comparison  of  Method  1  and  2. 

For  each  data  set  in  which  the  mfr  =  0,  the  complete 
results  are  presented  on  one  page;  for  mfr  not  equal  to  zero 
the  second  page  of  the  set  presents  the  results  for  the  natural 
co-ordinates  ( t  ,  4  ,  z)  . 

Each  set  of  data  is  specified  by  9  numbers.  These 
are:  x/ a ,  z/a,  X,  Y,  mfr,  day,  month  year,  program  path 
number.  The  9  numbers  are  listed  at  the  beginning  of  each 
set  of  results. 

An  interesting  fact,  that  is  not  shown  in  these  results, 
is  the  existence  of  negative  gage  total  pressure  readings  for  the 
extremities  of  the  jet.  Since  these  measurements  were  made 
with  the  total  pressure  probe  facing  the  flow  (x  -  y  -6  traverse 
device)  the  negative  gage  pressure  is  not  a  result  of  yaw  induced 
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